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VORWORT 


Dieser Band enthält 40 mathematische Programme aus den Bereichen 


— Mehr-Register-Arithmetik 
— Zahlentheorie 

— Kombinatorik 

— Algebra 

— Geometrie 

— Numerische Mathematik 


Besonders nützlich werden Praktiker und Computerfans die hier neu vorgelegte 
Langzahl-Arithmetik finden, die die Grundrechenarten für Zahlen bis 255 
Stellen gestattet! 


Zahlreiche Anwendungen finden auch die hier angegebenen kombinatorischen 
Prozeduren: Wer wollte nicht schon mal alle Permutationen oder Kombina- 
tionen eines Problems durchspielen? Vielfältige Unterstützung findet auch das 
Rechnen mit Polynomen, Matrizen und komplexen Zahlen: Neben dem kom- 
plexen Hornerschema werden insbesondere Algorithmen zur Polynomdivision 
und Matrizeninversion gegeben. 


Anwendungsbezogen sind die Programme der Numerischen Mathematik. Hier 
werden Verfahren zur Lösung von linearen und nichtlinearen Gleichungen ge- 
geben, zum Eigenwertproblem von Matrizen und zur numerischen Integration 
und Differentiation. Alle Programme werden durch vollständige Beispiele und 
komplette Programmausdrucke erläutert. 


Dem Verlag danken wir für die Herausgabe des Bandes und für die stets freund- 
liche Zusammenarbeit. 


Anzing, Burscheid, im Juni 1984 
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La machine arithmetique fait des effects 
qui approchent plus de la pensee 
que tout ce que font des animaux 

Pascal 


Daß die niedrigste aller Geistestätigkeiten 
die arithmetische ist, wird dadurch belegt, 
daß sie die einzige ist, welche auch durch 
eine Maschine ausgeführt werden kann 
Schopenhauer 


EINLEITUNG 


Mathematik und Computer 
Schon zu allen Zeiten hat es bemerkenswerte numerische Leistungen gegeben. 


Die berühmte Stieraufgabe, mit der Archimedes die Mathematiker von Alexan- 
dria zur Verzweiflung trieb, hatte keine kleinere Lösung als 


10366482, 7460514, 7358060 und 4149387. 
Ähnlich schwierig war die Lösung x = 1.3688081875 der Gleichung 
x3 +2 x2 + 10x = 20 


die Leonardo von Pisa, Fibonacci genannt, um 1220 auf 10 Dezimalen genau 
berechnete. 

Van Ceulen (1539-1610) benötigte mehr als die Hälfte seines Lebens, um die 
Zahl pi auf 35 Dezimalen zu berechnen. Als Dank dafür hat man sie in seinen 
Grabstein in Leiden gemeißelt. Fermat (1601-1655) stellt in Briefen an seine 
armen Zeitgenossen mehrfach die Frage nach den Lösungen von 


x2 — 149y2 = 1 
deren kleinste Werte x = 25801741449 und y = 2113761020 sind. 


In mühevoller Arbeit berechnete W. Shanks 1873 die ersten 707 Dezimalstel- 
len von pi. Ein Großteil dieser Mühe war jedoch vergeblich, da, wie sich später 
herausstellte, die 528ste Stelle bereits falsch war. 

An dieser Stelle muß auch der immense Rechenaufwand bedacht werden, den 
Ptolemaeus für die Erstellung seiner Sehnentafel Buergi und Brigg mit ihren 
Logarithmentafeln und Kepler mit den Rudolffinischen Tafeln hatte. 

Alle diese Probleme hätten in kurzer Zeit mit Hilfe eines Rechners gelöst wer- 
den können. Es erhebt sich die Frage, ob der Computer auch prinzipiell neue 
Fragestellungen in Angriff nehmen kann. 
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1976 kam es zu dem seltenen Ereignis, daß ein mathematischer Lehrsatz 
Schlagzeilen in der berühmten New York Times machte. Das sensationelle Er- 
eignis war, daß der Beweis des Vierfarbensatzes mit Hilfe von Computern ge- 
lungen war. 

K. Appel und W. Haken ernteten jedoch mit ihrem Beweis bei den meisten 
Mathematikern nur Kopfschütteln. Ein mathematischer Beweis war traditio- 
nell etwas, was man mit Papier und Bleistift nachvollziehen konnte. 


Wer aber sollte das 460-seitige Computerprogramm von Appel und Haken 
nachprüfen? Dazu kamen weitere 133 Seiten des Illionis Journal of Mathe- 
matics, auf denen die Autoren ihr Computerprogramm erläuterten. 


Das Vier-Farben-Problem, 1852 von Francis Guthrie aufgeworfen, besagt, 
daß jede ebene Landkarte mit 4 Farben so gefärbt werden kann, daß je 2 
benachbarte Länder verschiedenfarbig sind. Trotz erheblicher Anstrengungen 
widerstand die Vermutung über 120 Jahre lang den Beweisbemühungen der 
Mathematiker. 

Die Schwierigkeit war die, alle möglichen Lagemöglichkeiten der Länder in 
den Griff zu bekommen. Auf Grund der Vorarbeiten von Heesch (Hannover) 
gelang es Appel und Haken, alle möglichen Konfigurationen in einigen tau- 
send Fallunterscheidungen zu erfassen und vom Computer überprüfen zu 
lassen. 


Großen Aufschwung erfuhren die Teile der Mathematik, die sich mit endli- 
chen Mengen befassen, wie die Kombinatorik, Endliche Geometrie, Graphen- 
und Codierungstheorie, durch die Anwendung von Computern. Auch die 
Zahlentheorie erfuhr wesentliche Anregungen, so wurden z.B. sehr effektive 
Primzahltestverfahren entwickelt. Mit ihrer Hilfe konnte im September 1983 
die z.Z. größte bekannte Primzahl 

2132 049 _1 


mit 39751 Stellen gefunden werden. Die Rechenzeit beträgt bei den hier be- 
nützen Lucas-Lehmer-Tests nur wenige Stunden. 

Konträr ist auch die Meinung der Mathematiker zu den stochastischen Prim- 
zahltests, die Primzahlen nur einer gewissen Wahrscheinlichkeit nachweisen. 
Beim Rabintest läßt sich die Irrtumswahrscheinlichkeit bei einigem Aufwand 


auf 
10-60 


senken. Kann man dies als Beweis der Primzahleigenschaft ansehen? Oder ist 
nicht schon die Irrtumswahrscheinlichkeit eines Menschen größer? 


Bezeichnend ist die Geschichte von Rosser, Schönfeld und Yohe (Wisconsin), 
die bei ihrer Arbeit über die noch unbewiesene Riemannsche Hypothese einen 
Rechner benützten. Als ”reine’’ Mathematiker mißtrauten sie dem Computer 
und überprüften in mühseliger Kleinarbeit dessen Betriebssystem. Tatsächlich 
fanden sie mehrere Fehler in der internen Rechnerlogik. Diese Fehler hatte 
bisher noch niemand bemerkt, obwohl der Computer bereits seit mehreren 
Jahren in Betrieb war. 


10 


Im Gegensatz zu den Primzahltests fehlt es bei der Primzahlfaktorisierung 
noch an effektiven Algorithmen. So weiß man seit langer Zeit, daß eine Zahl 


wıe 
2256 + 1 


zerlegbar ist, man kennt jedoch bis heute noch nicht ihre Primzahlzerlegung. 


Die Primfaktorisierung hat große Bedeutung in der Kryptologie gewonnen. 
Man kann nämlich damit Verschlüsselungsverfahren konstruieren und ihr Prin- 
zip offenlegen, ohne befürchten zu müssen, daß es "‘geknackt”' wird (Verfah- 
ren von Rivest, Shamir und Adleman). Die Rechenzeit zur Faktorisierung 
einer 1000stelligen Zahl beträgt z.Z. noch 4 Mrd. Jahre. Dies hat große poli- 
tische Bedeutung, falls es einmal zu einem generellen Atomwaffen-Sperrver- 
trag kommen sollte und die Codierung der Beobachtungssonden bekanntgege- 
ben werden müßte. 
Ein weiteres bekanntes bisher ungelöstes Problem der Zahlentheorie ist der 
sog. große Fermat-Satz. Er besagt, daß es keine ganzzahlige Lösung der Glei- 
chung #0 

xn + yn = zn 
für n>2 gibt. Gleichungen, von denen nur ganzzahlige Lösungen gesucht wer- 
den, nennt man Diophantisch. Vor kurzem (1983) konnte Faltings (Wupper- 
tal) zeigen, daß Diophantische Gleichungen vom Grad >3 höchstens endlich 
viele Lösungen haben können (Mordellsche Vermutung). Durch die Reduzie- 
rung auf endlich viele Fälle scheint das Problem nun mit dem Computer greif- 
bar. 
Im Falle der sog. Catalanschen Vermutung ist dies tatsächlich der Fall. Diese 
besagt, daß es keine zwei Potenzen >2 von natürlichen Zahlen — außer 8 und 
9 — gibt, deren Differenz 1 ist. Tijdemann konnte ebenfalls zeigen, daß das 
Problem höchstens endlich viele Lösungen hat. Durch umfangreiche Compu- 
terberechnungen glaubt er die Catalansche Vermutung mit an Sicherheit gren- 
zender Wahrscheinlichkeit bewiesen zu haben. Was ist von einem solchen Be- 
weis zu halten? 
Auch viele Fragestellungen der Kombinatorik, die infolge des großen Rechen- 
aufwandes bisher ungelöst waren, konnten seitdem beantwortet werden. Be- 
kannt ist das Eulerische Offiziersproblem (1782): Ist es möglich, 6 verschiede- 
ne Offiziere aus 6 verschiedenen Regimentern so in einem Quadrat aufzustel- 
len, daß in jeder Reihe und Spalte jeder Offiziersrang und jedes Regiment ver- 
treten ist? 


Eine solche Anordnung nennt man ein Paar von orthogonalen Lateinischen 
Quadraten. Euler vermutete, daß es keine orthogonalen Quadrate der Ord- 
nung 6, 10, 14, usw. gäbe. Erst in den Jahren 1959/60 konnten Bose und 
Shrikhande nach umfangreichen Rechnungen orthogonale Lateinische Qua- 
drate der Ordnung 10 und 14 angeben. Sie erhielten dafür den Spitznamen 
“the Euler Spoilers’’. 

Verwandt mit den orthogonalen Lateinischen Quadraten sind die Blockpläne 
der Endlichen Geometrie. Ausgangspunkt war das Schulmädchenproblem des 
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Geistlichen Kirkman (1851): Kann man 15 Schulmädchen an den 7 Wochen- 
tagen so in Dreierreihen spazieren führen, daß jedes Mädchenpaar an genau 
einem Tag zusammentrifft? 

Spezielle Blockpläne sind die sog. Steinersysteme S{t,k,v). Dabei müssen 
v Dinge zu je t Stück in k Reihen angeordnet werden können. Die Lösung des 
Schulmädchen-Problems stellt somit ein (2, 3, 15)-Steinersystem dar. Es ist 
im Rahmen des Buches nicht möglich, die Vielzahl der ungelösten Probleme 
im Zusammenhang mit den Steinersystemen aufzuzählen. Erwähnt sei, daß 
die kleinsten Steinersysteme, deren Existenz bisher unbewiesen sind, folgende 


sind: 


S (4,5, 17), S (4,5, 21),S (4,5, 27) 


Durch umfangreiche Rechnung konnten Mendelsohn und Hung vor wenigen 
Jahren z.B. zeigen, daß das Steinersystem S (4, 5, 15) nicht existiert. Hier 
bleibt noch viel Forschungstätigkeit für künftige Mathematiker und viel 
Rechenarbeit für die nächste, noch schnellere Computergeration zu tun. 


N NEIN LAN 


I 
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"Aus Margarita Philosophica des Gregor Reisch, Freiburg 1503. 
Boethius mit indischen Zahlen und Pythagoras als Linienrechner 
vor der Arithmetica sitzend” 
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Es macht mich traurig, daß gebildete Leute 
nicht einmal von der Existenz 
meines Fachgebietes wissen 
P. R. Halmos 


Die Furcht vor der Mathematik .. 
steht der Angst erheblich näher 
als die Ehrfurcht 

F. Auerbach 


Mathematik — die ungeliebte Wissenschaft 
Kein Zweifel, die Mathematik war schon immer eine unpopuläre Wissenschaft. 


Besonders dem Christentum waren die heidnischen Schriften eines Euklid- 
Diophant und Archimedes als heidnisch suspekt. Tertullian schreibt: 
“Nach Christus brauchen wir keinerlei Wissbegier mehr, nach den 
Evangelien sind keinerlei Forschungen mehr nötig”. 


Als Folge der Verketzerung der Heiden, wurde die berühmte Mathematikerin 
Hypathia, Tochter des Akademievorstehers Theon, 415 von einer fanatischen 
Menge zu Tode gemartert. Dies, obwohl der Kirchenlehrer Synesios von 
Kyrene, Bischof von Ptolemais, ein Hörer ihrer Vorlesung war. Sogar die 
Mathematik hatte den Beigeschmack des Heidnischen. So lautete eines der 
Gesetze des Justinian-Kodex: “de Maleficis, mathematicis et ceteris similis’’ 
(Über Übeltätern, Mathematikern und dgl.): 
“Vollständig verboten ist die verdammenswerte Kunst der Mathe- 
matik’’ 
Die Mathematiker werden dort mit Giftmischern, Zauberern und Sterndeutern 
{Veneficis, Magis, Chaldeis) in einen Topf geworfen. Noch 1614 bezeichnete 
der Pater Caccini in einer Predigt, die gegen Galilei gerichtet war, die Mathe- 
matik als Teufelskunst und die Mathematiker als Urheber aller Ketzereien, die 
man aus allen Staaten vertreiben müsse. 
Er stützte sich dabei auf den heiligen Augustinus, der in De genesi ad literam 
geschrieben hatte: 
“Es besteht die Gefahr, daß die Mathematiker mit dem Teufel im 
Bunde, den Geist trüben und in die Bande der Hölle verstricken’’, 


Als 529 Kaiser Justinian die Athener Philosophenschule schloß, erlosch eine 
lange Tradition mathematischer Lehre. Noch bis 485 hatte an dieser Schule 
Proklos gelehrt, dem wir einen grundlegenden Euklid-Kommentar verdanken. 
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Es ist ein ganz wesentlicher Verdienst der arabischen Kultur ab 800 die wich- 
tigsten Werke des Apollonius, Archimedes, Heron, Euklid, Diophant, Ptole- 
maeus usw. in ihre Sprache übersetzt und somit der Nachwelt erhalten zu 
haben. Einige Übersetzungen aus dem Griechischen ins Lateinische führte 
auch Boethius (480 - 524) durch. Von ihm stammt auch die im Mittelalter 
übliche Einteilung der Künste in das Trivium (Grammatik, Rhetorik, Logik) 
und das Quadrivium (Arithmetik, Geometrie, Astronomie und Musik). 


Um eine mathematische Allgemeinbildung in Gang zu bringen, führte der 
Schotte Alkuin von York, im Auftrag Karls des Großen, 789 eine Schul- 
reform durch. Sie forderte, daß neben Psalmen, Noten und Grammatik ins- 
besondere auch der Computus gelehrt werden solle. Unter dem Computus 
Ecclesiaticus verstand man die Berechnung der beweglichen Kirchenfeiertage. 
Es war nämlich der Kirche ein Dorn im Auge, daß wegen der mangelhaften 
Rechenfertigkeit der Geistlichen das Osterfest zu verschiedenen Zeitpunkten 
gefeiert wurde. Alkuin selbst schrieb ein Rechenbuch zu Unterrichtszwecken. 


Der Rechenunterricht scheint — wie heute — nicht die besondere Freude der 
Schüler gewesen zu sein. Ein Mönch namens Strabo der Klosterschule Reiche- 
nau schreibt 822 in sein Tagebuch: 


‘Zur Abwechslung und Unterhaltung lösten wir mathematische Rät- 
sel, welche Alkuin für den großen Karl gefertigt hatte. Viele vermoch- 
ten nicht allen diesen Rechnungen zu folgen, und bevor wir zur Geo- 
metrie übergingen, traten diejenigen aus, welche sich fortan den 
Studien der Medizin, Rechtswissenschaft und den Künsten der Male- 
rei und der Bildhauerei widmen wollten’, 


Ein Mathematiker, Gerbert von Aurillac (940 - 1003) bestieg als Sylvester II 
999 den Stuhl Petri. Gerbert hatte einige Jugendjahre in Spanien verbracht 
und dort von den Arabern den Abakus (Rechenbrett) kennengelernt. Als Leh- 
rer der Domschule Reims verfaßte er ein Lehrbuch zum Rechnen am Abakus. 
Auch an den aus den Domschulen hervorgegangenen Universitäten Bologna 
(1088), Paris (1150), Salerno (1173) und Montpellier (1180) spielte die 
Mathematik keinerlei Rolle. 


Von der Universität Erfurt weiß man aus Quellen, daß die dort wirkenden 
Humanisten die Anstellung eines Mathematikers strikt ablehnten. Noch An- 
fang des 16. Jahrhunderts legten an der Universität Paris die Kandidaten für 
den Grad des Magisters der Künste kein Examen im Fach Geometrie ab. Sie 
mußten lediglich beeiden (l), Vorlesungen über die ersten 6 Bücher des Euklid 
gehört zu haben. 

Das Risiko einer wirklichen Prüfung wollte man nicht eingehen. Bezeichnend 
ıst, daß der Lehrsatz des Pythagoras am Ende des 1. Kapitels den Namen 
”"Magister Matheseus’’ trug und somit als Weisheit letzter Schluß galt. 


Welche Schwierigkeiten das Rechnen damals machte, wird aus der Rede 
Philipp Melanchtons vor der Universität Wittenberg (1517) deutlich. 
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".... Deshalb können ihre Anfangsgründe (der Rechenkünste) gar 
nicht dunkel und schwer sein; sie sind im Gegenteil so durchsichtig, 
daß Kinder sie begreifen können, weil ja alles so natürlich vor sich 
geht. Die Regeln des Vielfachen und Teilens allerdings erfordern viel 
mehr Fleiß, aber ihr Sinn wird schon bald von den aufmerksameren 
eingesehen werden. Übung und Anwendung erfordert diese Fertig- 
keit wie alle anderen.’’ 
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“Ein Vater meldet seinen Sohn bei einem Rechenmeister an. 
Holzschnitt von Hans Weiditz, 1535‘ 


Die Mathematik-Kenntnisse dieser Studenten aus den Lateinschulen sind ent- 
sprechend schlecht gewesen. Der Rechenunterricht wurde in der Regel von 
einem externen Rechenmeister übernommen, die dazu meist eigene Rechen- 
bücher verfaßten. Der bekannteste Rechenmeister war Adam Riesle) in Erfurt 
und Annaberg. 

Es dauerte noch sehr lange, bis Mathematik ein reguläres Unterrichtsfach war. 
So wurde noch 1805 (!) im Jahresbericht eines lutherischen Gymnasiums in 
Essen vermerkt, daß an der Schule kein Rechenunterricht erteilt werde. Dies 
sei der Fall, weil ältere Schüler den Lehrern bekannt gemacht hätten, ’’daß in 
der Stadt jetzt ein fertiger Rechenmeister sei, der eine kürzere Methode habe’. 


Oft wurde der Mathematik-Unterricht von nicht ausgebildeten Lehrkräften, 
z.B. von Theologen übernommen, die in Wartestellung auf eigene Pfarrpfründe 
waren. Ob es wohl diese Theologen waren, über die es in der Prüfungsordnung 
von 1820 für Preussen heißt: Um das Eindringen untüchtiger Objekte in den 
Höheren Schulen Einhalt zu gebieten... .'’. Auch in Bayern waren viele Mathe- 
matik-Lehrstühle im 19. Jahrhundert noch mit Theologen besetzt. 
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Daneben litt die mathematische Lehre noch unter anderen Schwierigkeiten. 
Mollweide in Leipzig erklärte es für unmöglich, neben der, für alle Fakultäten 
verbindlichen Mathesis pura, auch höhere Mathematik zu bringen, weil es da- 
bei zuviel Schreibens an der Tafel gäbe. 

Bezeichnend ist auch der Bericht von Neumann in Berlin von 1818: "Als ich 
mich beim Professor für Mathematik meldete, sagte dieser, ja ich habe die 
Vorlesung angezeigt, aber sie pflegt niemals zustande zu kommen. Ich verab- 
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“Darstellung der Arithmetica und Geometria aus dem Quadrivium, 
Holzschnitte um 1500 
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redete mich mit 5 anderen, zu ihm zu gehen. Der Professor kam ins Audito- 
rium und schrieb ununterbrochen Formeln an die Tafel und sprach dabei kein 
Wort, bis die Zeit um war. Am 2. Tag kamen nur noch 2 Zuhörer. Der Profes- 
sor stellt sich wieder an die Tafel und zeichnete wieder ununterbrochen 
mathematische Formeln an diese. Am dritten Tag kam außer mir nur noch 
ein Zuhörer. Der Professor erschien, ging ans Katheder, wandte sich an uns 
und sagte: Sie sehen, meine Herren, es kommt kein Kolleg zustande’’. 


"Darstellung der Musica und Astronomia, Holzschnitte um 1500’ 
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1. Addition langer Zahlen 


Jeder, der schon länger mit dem Computer rechnet, ist schon einmal an die 
Grenze der Rechengenauigkeit gestoßen. Während bei Großrechenanlagen 
meist eine Langzahl-Arithmetik hardware-mäßig implementiert ist, ist dies bei 
Mikrocomputern nicht der Fall. 

Daher sollen im folgenden, Programme zu den Grundrechenarten für lange 
Zahlen angegeben werden. 


Diese Programme setzen voraus, daß die eingegebenen Zahlen kommafrei und 
positiv sind. Dies ist keine große Einschränkung, da die Vorzeichenregeln 
leicht von Hand ausgeführt werden können. Kommafreie Zahlen können 
durch entsprechendes Ausklammern von Zehnerpotenzen erhalten werden, 
z.B. 

12.3456789 = 123456789 + 10°? 

0.008987654 = 987654 * 10° 


Alle Zahlen werden als Stringvariablen eingegeben und können daher bis zu 
255 Stellen haben, sie werden innerhalb des Programms in Felder zerlegt. Da- 
bei kann jedes Feld 1 bis 8 Ziffern umfassen. 

Umfaßt ein Feld mehrere Ziffern, so ist zu beachten, daß der Rechner führen- 
de Nullen nicht ausdruckt. Mit Hilfe der Stringfunktion kann man auch füh- 
rende Nullen ausdrucken lassen, wie es z.B. im Programm 7 geschieht. 


Wird der 1. Summand durch das Feld U(l), der zweite durch V(I) dargestellt, 
so kannbei ziffernweiser Addition das Addieren wie folgt beschrieben werden: 


98 K=0 

108 FOR I=NTO1STEP -1 
119 w= UN +VÜI+K 

128 W(l) =W —- INT(W/10) * 10 
138 K = INT(W/19) 

1498 NEXTI 

158 w{@)=K 


Dabei ist N die Länge des größten Summanden, K der Übertrag. Die Summe 
ist im Feld W(I) gespeichert. Bei sehr großen Rechnungen, wie z.B. bei der 
Berechnung von 1080 Dezimalzahlen, zerlegt man die Summanden z.B. in 
Fünferblöcke U(l) bzw. V(l). 
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Die Fünferblock-Addition läuft wie folgt ab: 


99 K=0 

188 FOR I=N TO 1STEP —1 
198 w=U() +V() +K 

1108 K = INT (W/1015) 

129 w(l)=W-K + 1015 
136 NEXTI 

1498 w(0)=K 


Dabei ist N natürlich die Anzahl der Fünferblöcke. 


Zum folgenden Programm 
Als Beispiel werden im Programm die Zahlen: 


7839412596412348 
und 
914935812769267 
addiert. Hier ergibt sich die Summe 
8744448439181615 


Das hier gezeigte Programm arbeitet mit Einerblöcken, so daß auch Nullen 
richtig ausgegeben werden. 
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REM ADDITION LANGER ZAHLEN 


HOME : INVERSE 
PRINT " n 
PRINT " ADDITION LANGER ZAHLEN " 
PRINT " " 
NORMAL : PRINT : INVERSE 
PRINT "1. SUMMAND "5: NORMAL 
PRINT : PRINT 
INPUT A$: PRINT : INVERSE 
PRINT "2. SUMMAND "3: NORMAL 
PRINT : PRINT 
INPUT B$ 
PRINT 
N= LEN (A$):M = LEN (B$) 
IFEN> = M THEN 198 
H=N:IN = MM =H 
HS = A$:A$ = B$:B$ = Hs 
DIM UN) ,ViN) WIN) 
IF N = M THEN 269 
FORI=S1TON-M 
B$ = "g" + B$ 
NEXT I 
REM UMWANDLUNG IN ZAHLEN 
FOR I=1TON 
U$ = MID® (A$,1,1) 
IF Us < "O" OR Us > "9" THEN 538 
UT) = VAL (U$) 
V$ = MID$ (B$,1,1) 
IF V$ < "9" OR Vs > "9" THEN 539 
VD) = VAL (v8) 
NEXT I 
REM STELLENWEISE ADDITION 
J=NK=9 
W=UL) +VID) +K 
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392 WiJ) = W- INT (W / 19 
488 K= INT (W / 19 


410 J 
420 


=J-14 
IF J > 8 THEN 388 


4358 W(iO) = K 


PRINT 


INVERSE : PRINT "SUMME = 


PRINT : PRINT 
PRINT " "; 


489 K=1-K 


FOR I=KTON 

PRINT W(DD5 

NEXT I 

GOTO 578 

PRINT © PRINT CHR$ (7) 
INVERSE 


PRINT " FEHLERHAFTE EINGABE ! 


NORMAL 
PRINT © PRINT 
END 


ADDITION LANGER ZAHLEN 


1. SUMMAND 


7830412596412348 


2. SUMMAND 


9149835842769267 


SUMME = 


8744448459181615 
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* 18 


NORMAL 


2. Subtraktion langer Zahlen 


Die Subtraktion langer Zahlen läuft als inverse Rechenoperation ähnlich wie 
die Addition ab. Da, wie anfangs erwähnt, auf das Rechnen mit Vorzeichen 
verzichtet wird, muß vorausgesetzt werden, daß der Minuend größer ist als der 
Subtrahend. Ist dies nicht der Fall, so vertauscht man Minuend mit Subtra- 
hend und versieht die Differenz mit einem negativen Vorzeichen. 


Wird der Minuend durch das Feld U(l), der Subtrahend durch V(I) dargestellt, 
so kann bei ziffernweiser Subtraktion das Subtrahieren wie folgt beschrieben 
werden: 


989 K=0 

188 FOR I=N TO 1STEP —1 
119 w= U) - VI +K 
128 K = INT (W/10) 

138 wi)=W-K * 10 

1408 NEXTI 


Dabei ist N die Länge des Minuenden, K der Übertrag. Die Differenz ist im 
Feld W(I) gespeichert. Übersteigt der Subtrahend den Minuend, so wird eine 
entsprechende Fehlermeldung ausgegeben. 

Bei sehr großen Rechnungen, wie z.B. bei der Berechnung von Fünferblöcke 
UN) bzw. V(l). 

Die Fünferblock-Subtraktion läuft wie folgt ab: 


99 K=0 

188 FOR I=N TO 1STEP —1 

1898 w= UN) - VI —K 

118 IFW>=BTHENK = dB: GOTO 130 
128 K=1:W=W + 1015 

138 w(il)=W 

149 NEXTI 


Dabei ist N natürlich die Anzahl der Fünferblöcke. 


Zum folgenden Programm 


Als Beispiel werden im Programm die Zahlen: 


7841035 16910456 

und 
9541217036428 

subtrahiert. Es ergibt sich die Differenz 
774562299874028 
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1988 REM SUBTRAKTION LANGER ZAHLEN 
110 : 

128 HOME 

139 INVERSE 

148 PRINT " " 
158 PRINT " SUBTRAKTION LANGER ZAHLEN " 
168 PRINT " n 
178 NORMAL : PRINT 

188 INVERSE : PRINT " MINUEND: " 

198 NORMAL 

280 PRINT 

218 INPUT As 

228 PRINT : PRINT : INVERSE 

238 PRINT " SUBTRAHEND: " 

248 NORMAL : PRINT 

258 INPUT B$ 

268 PRINT 

270 : 

280 N= LEN (A$):M = LEN (B$) 

298 IF N< M THEN 568 

388 DIM U(N),V(N),WiN) 

3190 : 

328 IF N = M THEN 378 

338 FORI=1T0ON-M 

340 B$ = "9" + B$ 

350 NEXT I 

368 : 

378 REM UMWANDLUNG IN ZAHLEN 

388 FOR I=1 TON 

398 U$ = MID$ (A$,1,1) 

480 U(I) = VAL (U$) 

419 IF U$ < "@" OR U$ > "9" THEN 498 
420 V$ = MID$ (B$,1,1) 

438 IF Vs < "@" OR V$ > "9" THEN 498 
448 V(l) = VAL (V$) 

450 NEXT I 

460 : 

478 REM STELLENWEISE SUBTRAKTION 
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NK=9 

U) - VI) +K 
INT (W / 19) 

‚)=W-K+ 19 

J 

J 


Lo) 
nıcuı mn u 


BEZEARZCO 


-1 
IF > 8 THEN 498 
WB) = K 
IFK< > - 1 THEN 69298 
INVERSE 
PRINT CHR$ (7)5" MINUEND KLEINER ALS SUBTRAHEND ! 
NORMAL : END 
INVERSE 
PRINT " DIFFERENZ = ": NORMAL 
PRINT : PRINT " "5; 
K=9 
IF W(K) = 8 THENK = K + 1: GOTO 848 
FOR I =K TON 
PRINT W(DD5 
NEXT I 
GOTO 728 
PRINT : INVERSE 
PRINT CHR$ (7)5" UNERLAUBTE EINGABE ! " 
NORMAL 
PRINT : PRINT 
END 
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SUBTRAKTION LANGER ZAHLEN 


MINUEND: 


784183516912456 


SUBTRAHEND: 
9541217056428 
DIFFERENZ = 


774562299874928 
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3. Multiplikation langer Zahlen 


Nicht ganz so einfach wie die Addition und Subtraktion ist die Multiplikation 
langer Zahlen. 

Wird der größere Faktor stellenweise in das Feld U(l) und der zweite in V(l) 
zerlegt, so läßt sich die stellenweise Multiplikation wie folgt programmieren: 


198 FORJ=MTO1STEP —1 

118 IF V(J) =8THEN W(J) = GB: GOTO 190 
120 K=0 

138 FOR I=N TO 1 STEP —1 

148 T=U() * VII) +WÜHI +K 

158 W(l+J) =T — INT(T/10) * 10 

168 K = INT(T/10) 

178 NEXTI 

188 wi) =K 

198 NEXT J 


Dabei ist M die Stellenzahl des größeren Faktors, N die des kleineren Faktors. 
Im Feld W(I) wird die Ziffernfolge des Produkts gespeichert. 

Ähnlich, wie bei der Addition und Subtraktion, könnte man auch die Multi- 
plikation in Dreier- bzw. Viererblöcken durchführen. 


Zum folgenden Programm 


Das Programm berechnet nach dem oben angegebenen Verfahren das Produkt 
zweier langer Zahlen. 


Als 1. Beispiel wird die Zahl 
999999999999999 

quadriert. Es ergibt sich 
999999999999999809 00 090900000 1 


Im 2. Beispiel werden die beiden Primfaktoren 
59649589127497217 

und 
5704689200685129054721 
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der Zahl 
2128 +1 
miteinander multipliziert. Diese Primzahlfaktorisierung gelang erst 1970 den 


Amerikanern Morrison und Brillhart, obwohl man schon sehr lange wußte, 
daß diese Zahl keine Primzahl ist. 


Das Programm liefert 


340 82366 92093 84634 63374 60743 17682 11457 
(vgl. Programmausdruck). 


olıg L_ 
Abc- mir leiß das cinmalein Sownet 
be 
N 


ala 411618 
9 I12lı slı 82 1 247 
sB I1zlieze 283 215 


“"Multiplikationstafel aus der Arithmetik von Johannes 
Widmann, Leipzig (1489) 
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REM MULTIPLIKATION LANGER ZAHLEN 


HOME : INVERSE 

PRINT " " 
PRINT " MULTIPLIKATION LANGER ZAHLEN " 
PRINT " " 
NORMAL =: PRINT : INVERSE 

PRINT " 1.FAKTOR: "= NORMAL 

PRINT 

INPUT A$ 

PRINT : INVERSE 

PRINT " 2.FAKTOR: "2 NORMAL 

PRINT 

INPUT B$ 

PRINT 

REM GROESSEREN FAKTOR BESTIMMEN 
N= LEN (A$):M = LEN (B$) 

IFN> = M THEN 312 
H$ = A$:A$ = BS:B$ = H$ 
H=NIN = MM =H 

DIM U(N),V(M,W(M + N 


REM UMWANDLUNG IN ZAHLEN 

FOR I = 1 TON 
U$ = MID$ (A$,1,1) 

IF US < "9" OR Us > "9" THEN 708 
UT) = VAL (U$) 

NEXT I 

FOR I=1 TOM 
V$ = MID$ (B$,1,1) 

IF v8 < "Q" OR v& > "9" THEN 708 
V(TD) = VAL (v8) 

NEXT I 

FORI=M+1TOM+N 
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REM STELLENWEISE MULTIPLIKATION 

IF V(J) = 8 THEN W(J) = 8: GOTO 588 
I=N:K=09 
T=UMD *« VO) + WI + JI +K 
WII +J) =T- INT (T / 18 + .DOBBBl) * 10 
K= INT (T / 19 


I=1-1 

IF I > 28 THEN 528 
Wi) = K 
J=J-1 


IF J > 8 THEN 598 

INVERSE 

PRINT " PRODUKT = " 

NORMAL =: PRINT ©: PRINT " "; 
K=1 

IF WiK) = 8 THENK = K + 1: GOTO 559 
FORI=SKTON+M 

PRINT W(DDJ 

NEXT I 

GOTO 739 

INVERSE 

PRINT CHR$ (7)5" UNZULAESSIGE EINGABE ! 
NORMAL 

PRINT : PRINT 

END 
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MULTIPLIKATION LANGER ZAHLEN 


1. FAKTOR: 
?IIIIIIIIIIIII79 

2. FAKTOR: 
PIIITIIIIIIIIITI7 

PRODUKT = 


9999997979979797977782902202202209091 


1. FAKTOR: 
75704689200685129054721 
2. FAKTOR: 
759649589127497217 
PRODUKT = 


3492823646920938453463374607431768211457 


1. FAKTOR: 
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4. Division langer Zahlen 


Wesentlich komplizierter als Addition und Multiplikation ist die Division lan- 
ger Zahlen. 

Das Programm führt eine stellenweise Division durch, entsprechend dem 
schriftlichen Rechnen. Zuerst wird geprüft, wie oft der Divisor in den Dividen- 
den hineingeht. Sodann wird der Divisor mit dem entsprechenden Vielfachen 
multipliziert. Diese Multiplikation einer langen Zahl mit einem einstelligen 
Faktor wird bei Programm 5 erklärt. Das sich ergebende Produkt muß noch 
vom Dividenden subtrahiert werden. Der genaue Algorithmus kann dem fol- 
genden Programm entnommen werden. 


Zum folgenden Programm 

Als Beispiel wird die schon aus dem Programm 3 bekannte Zahl 
2128 +1 

durch den Primfaktor 
59649589127497217 

dividiert. Das Programm liefert den ganzzahligen Quotienten 
5704689200685129054721 


Das Programm ist so aufgebaut, daß Zahlenwerte sowohl aus DATA-Zeilen, 
als auch über INPUT-Anweisungen übernommen werden können. Sollen DATA- 
Zeilen gelesen werden, dann müssen die Programmzeilen 199 und 148 entfallen 
bzw. durch eine vorgestellte REM-Anweisung unwirksam gemacht werden. 
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REM DIVISION LANGER ZAHLEN 


HOME 
INVERSE 
PRINT " 


PRINT " DIVISION LANGER ZAHLEN 


PRINT " 
NORMAL : PRINT «: 


INVERSE : PRINT " DIVIDEND " 
INPUT A$: GOTO 128 
READ A$: PRINT A$ 


PRINT : PRINT 


INVERSE : PRINT " DIVISOR 
INPUT B$: GOTO 169 
READ B$: PRINT B$ 


PRINT : PRINT 


LEN (AS):N = 
-N 


FOR I=-1T0OM+N 
VAL ( MID$ (A$,I,1)) 


UI) = 
NEXT I 
FOR I=1 TON 

von = 
NEXT I 


D 


K = 


FOR I=N +M TOD 1 STEP 


U=SUM «D+rK 
K= INT (U7/ 19 
UD =U-K%* 19 
NEXT I 

U) =K 
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LEN (B#) 


= M 
FMX< 8 THEN PRINT CHRS 
IM UN + M),VIN),WIN),QM 


VAL ( MID$ (B$,1,1)) 


(Ms"EINGABEFEHLER !": 


= INT (19 7 (VD + DISUID 
IFD= 1 THEN 478 
=® 


1 


398 : 

409 K = 9 

418 FOR I =N TO 1 STEP - 1 

420 VeV(ID) «D+rk 

438 K= INT (V / 1 

449 V(I) =V-K #19 

4598 NEXT I 

468 VID) = K 

478. 

A809 J= 98 

498 IF UI) = VYt1) THEN Q = 9: GOTO 529 

s80 IF Vti) = 8 THEN PRINT CHR$ (7)5"FUEHRENDE NULL !": END 
siea= INT (UM *12 +UJ + 1)) ZVAD 

328 @1 = (UV) #12 + U + DD - AH VI) * 10 + U(J + 2) 
538 IF v(2) *#@< = Qi THEN 569 

349 Q = Q - 1: GOTO 529 


Sad K= 9 

578 FOR I=N TO 1 STEP -1 

SB W= VD *«Q+rK 

5978 K= INT (W / 19 

688 WID =W-K#* 19 

619 NEXT I 

628 WO) = K 

658 : 

649 K= BF =9 

658 FOR I = J +NTO JSTEP -1 
660 U = UI) - WII - J) HK 

678 K= INT (U / 19 

688 UIID = U-K * 19 

&978 NEXT I 

#8 IFK= -1 THENF = 1 

7118: 

728 @(J) = Q 

738 IF F = 8 THEN 858 

732 Q(J) = QA(J) - I 

759 K=09 

768 FOR I=J+NTOJSTEPr -1 
778 U =UNM +V(I-J +K 

788 K= INT (U 7/7 19 


35 


79 UM) = U-K #19 


NEXT I 


U - DD =UlJ-UVD)+K 


= ) + 
IF J< = M THEN 498 


7) 
I=M+1i1T7T0OM+N 
= U(I) +K # 18 
(I = INT(U/D 
= UÜ- INTU/D *«D 


= 0%: REM AUSGABE 

INVERSE : PRINT " QUOTIENT ";2 NORMAL 
PRINT ” "5 

IF Q(J) = 8 THEN J = J + 1: GOTO 988 
FOR I=J TOM 
PRINT Q(TI}$ 
NEXT I: PRINT "."; 
FORI=M+1T0ON+M 


PRINT U(D); 

NEXT I 

PRINT : PRINT 

END 

DATA "3492823689209584534533574609743178688211457" 
DATA "59649589127497217" 
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DIVISION LANGER ZAHLEN 


DIVIDEND 34028236492089738463463374697431768211457 


DIVISOR 59649589127497217 


QUOTIENT 57846889208685127054721.2000800090000290090 
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5. Potenzieren 


Eine weitere Anwendung der Mehr-Register-Arithmetik liefert das Potenzie- 
ren, da hier sehr schnell große Zahlen auftreten. 

Das Potenzieren wird im Programm als wiederholte Multiplikation durchge- 
führt. Eine andere Möglichkeit wäre das logarithmische Rechnen, jedoch ist 
die mehrfachgenaue Berechnung von Logarithmen noch aufwendiger. 


Teilt man die mit B zu multiplizierende Zahl in Achterblöcke R{l), so kann 
die Multiplikation wie folgt durchgeführt werden: 

1098 U=® 

118 FOR I=S1TOR 

128 H=R({ *B+U 

138 IFH <= -1E8 THEN U = 8: GOTO 169 

148 U= INT (H/1E8) 

158 H=H-U* 1E8 

168 R{l) =H 

178 NEXTI 


Dabei ist R die Anzahl der Achterblöcke und U der jeweilige Überlauf. 


Soll die Zahl B zur N-ten Potenz erhoben werden, so muß das angegebene 
Programmstück noch in die Schleife 


FORK=-1TON..... NEXTK 
eingebaut werden. 
Die benötigte Anzahl der Achterblöcke wird logarithmisch berechnet: Der 
Zehnerlogarithmus der gesuchten Potenz BN ist 

N * LOG (B)/LOG (19) 


Da der aufgerundete Zehnerlogarithmus einer Zahl die Stellenzahl angibt, 
kann daraus durch Division mit 8 die notwendige Registerzahl R berechnet 
werden. 
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Zum folgenden Programm 


Das Programm berechnet nach dem oben beschriebenen Verfahren beliebige 


Potenzen mit einstelliger Basis. 
Soll die Basis mehrstellig sein, so wird empfohlen die Registerzahl entspre- 
chend zu vermindern. Dazu müssen alle Achter im Programm durch 7 bzw. 6 


ersetzt werden. 
Als Beispiel wird die hundertste Potenz von 2 berechnet. Es ergibt sich: 


1267650 60822822 94914967 032053756 


(vgl. Programmausdruck). 
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REM BERECHNUNG GROSSER POTENZEN 


HOME =: INVERSE 


PRINT " " 
PRINT " BERECHNUNG GROSSER POTENZEN " 
PRINT " " 
NORMAL : PRINT : PRINT 

INVERSE : PRINT " BASIS "j: NORMAL 


INPUT B: INVERSE : PRINT : PRINT 
PRINT " EXPONENT "5: NORMAL 
INPUT N 


PRINT : PRINT : INVERSE 
PRINT ss "B; 1 HOCH "N; N m 1 
NORMAL : PRINT : PRINT 

DIM R(199) 

REM ZAHL DER 8-ER BLOECKE 
R=N* LOB (B) / LOB (19) 
R= INT(R/B +1 

REM MULTIPLIKATIONSSCHLEIFE 
U=@RU) =1 

FORK = 1 TON 

FOR I = 1 TOR 
H=RUD#B+U 

IFH< - 1ES THEN U = 9: GOTO 388 
U= INT (H / 1EBO 
H=H-U+* 1E8B 
RD =H 

NEXT 1 

NEXT K 

REM AUSGABE 
Fg nz " 

FOR I=RTOALSTEP -1 


R$ = STR$ (R(D)) 


IF LEN (R$) < 8 THEN R$ = F$ + R$: GOTOD 469 


PRINT R$5" "3 
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4898 IF(iR-I+1) /4= INT(CR-TI+1) 74) THEN 


490 FR = "9" 

590 NEXT I 

518 PRINT : PRINT 
520 END 


BERECHNUNG GROSSER POTENZEN 


BASIS ?2 


EXPONENT 71929 


2 HOCH 198 = 


1267659 60022822 94014967 83295376 
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PRINT 


6. FAKULTATS-BERECHNUNG 


Die Berechnung der Fakultäts-Funktion liefert ebenfalls eine Anwendung der 
Mehr-Register-Arithmetik. Die Fakultätsfunktion n!, definiert als das Produkt 
aller natürlichen Zahlen von 1 bis n, wächst nämlich noch schneller als eine 
Potenzfunktion. 

Die Fakultätsberechnung wird im Programm als wiederholte Multiplikation 
durchgeführt. Eine andere Möglichkeit wäre auch hier das logarithmische 
Rechnen, jedoch ist die mehrfachgenaue Berechnung von Logarithmen noch 
aufwendiger. Die Multiplikation erfolgt nach dem bei Programm 5 angegebe- 
nen Verfahren. 

Komplizierter ist nur die Berechnung der benötigten Register. Sie erfolgt 
ebenfalls logarithmisch. Da der Logarithmus eines Produkts gleich der Sum- 
me der Logarithmen ist, ergibt sich log(n!) aus 


198 FOR I=2 TON 
118 L=L+LOG(I) 
120 NEXT I 


Da die so ermittelte Zahl L, aufgerundet die Stellenzahl angibt, ist L/6 die 
benötigte Anzahl von Sechserblöcken. Da der Computer jedoch mit natürli- 
chen Logarithmen rechnet, muß die Zahl zuvor über den Faktor 

In (10) = .4342945 


auf den Zehnerlogarithmus umgerechnet werden. 


Zum folgenden Programm 


Das Programm berechnet die Fakultätsfunktion für dreistellige Zahlen in 
Sechserblöcken. Ist die Fakultät einer größeren Zahl gesucht, so wird aus 
Genauigkeitsgründen empfohlen mit Fünfer- bzw. Viererblöcken zu rechnen. 
Dazu müssen die im Programm auftretenden Sechsen entsprechend ersetzt 
werden. 


Als Programmbeispiel wird 
190 ! 
berechnet. Das Ergebnis kann dem Programmausdruck entnommen werden. 
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REM FAKULTAET GROSSER ZAHLEN 


HOME : INVERSE 

PRINT " ” 
PRINT " FAKULTAET GROSSER ZAHLEN " 
PRINT " " 
NORMAL © PRINT = PRINT : INVERSE 
PRINT " WELCHE FAKULTAET "5: NORMAL 
PRINT : PRINT 

INPUT N 

FRINT : PRINT 

INVERSE 

PRINT " Da N Fa Bu Bu u 

NORMAL : PRINT 

DIM R(198) 

REM ZAHL DER 4-ER BLOECKE 
L=1l 
Ril) = 
FOR I=2T0ON 
L=L+ LOG (D 

NEXT I 
L=L * „4342945 
R= INT (L/6) +1 


1 


FOR I=NTO2STEP -1i 
L=L+ LOG (I * „A342945 
R= INT(L/d +1 


REM MULTIPLIKATIONSSCHLEIFE 
U=8 


U 
THEN U = 8: GOTO 419 
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sS7 U = INT (H / 1E6) 

6 H=H- U%* I1ES 

418 RU) =H 

420 NEXT J 

450 NEXT I 

449 : 

458 REM AUSGABE 

420 FORI=RTOAISTEP - 1 

478 IF R{I) = 8 THEN PRINT "908008 ":: GOTO 518 
480 R$ = STR$ (R(TD)) 

4798 IF LEN (R$) < 6 THEN R$ = "8" + Rs: GOTO 490 
S20 PRINT Rss" "5 

318 IFÄR-I+1) /5= INTCR-I+1) 75) THEN PRINT 
3208 NEXT I 

358: PRINT : PRINT 

548 END 


FAKULTAET GROSSER ZAHLEN 


WELCHE FAKULTAET 


?198 


198 | = 


289995 3265215 443944 152681 6997238 
856246 708490 715968 264381 521468 
392965 895217 599993 229915 628941 
463976 156518 286253 677920 827223 
758251 185218 916864 00022099 290009 
2OoOODa 200009 
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7. mAUF 1999 STELLEN 


Als weitere Anwendungen der in den Programmen 1 — 4 angegebenen Mehr- 
Register-Arithmetik werden noch 2 mathematische Konstanten auf jeweils 
18090 Stellen berechnet. 

Die Kreiszahl 7 , definiert als 


Kreisumfang/Durchmesser 


ist wohl die berühmteste aller mathematischen Konstanten. Sie ist auch kultur- 
historisch interessant, da die Genauigkeit, mit der ein Volk 7 bestimmt hat, 
Auskunft über das erreichte mathematische Niveau gibt. 

Während in älteren babylonischen Keilschrifttafeln die Kreiszahl zu 3 be- 
stimmt, zeigen die 1936 in Susa ausgegrabenen Tafeln den Wert 


3 7'306 

Rechnet man diesen Wert aus dem 6ßer-System um, so ergibt sich 
3 + 7/60 + 30/3600 

oder 
3 1/8. 


Bemerkenswert genau ist auch der ägyptische Wert: Der Papyrus Rhind gibt 
n als das (5) 2_fache des Durchmessers an. Dies liefert 


3 13/81. 


Die semitischen Völker rechneten mit dem Wert 3, der sich zweimal in der Bibel 
findet: 

“Es maß 16 Ellen von einem Rand zum anderen. Eine Schnur von 30 Ellen 
umspannte es ringsum’ (1. Könige 7,23; 2. Chronik 4,2). Archimedes (287 — 
212 v. Chr.) gibt die Schranken 


223/71 <a <22/7 


an. Ptolemäus rechnet 400 Jahre später mit dem Mittelwert aus diesen Schran- 
ken 
3 17/120. 


Um 500 n. Chr. fanden die Chinesen mit dem sehr genauen Wert 
355/113 = 3.1415929.. 

Die Inder sahen zur gleichen Zeit 
3927/1259 = 3.1416 


als genauen Wert an, rechneten aber meist mit 


v 16 =3.1623.. 


47 


Erst 1610 gelang van Ceulen, wie in der Einleitung erwähnt, mit 35 Dezimal- 
stellen eine wesentliche Genauigkeitssteigerung. 


In der Folgezeit wurde m zunehmend genauer berechnet. Zuletzt gelang es 
E. Salamin 1976 mit Hilfe eines Großcomputers auf 33 Millionen Dezimal- 
stellen zu berechnen. Natürlich hat diese Genauigkeit keinerlei praktische Be- 
deutung, die Dezimalbruchentwicklung wurde jedoch als idealer Zufallszahlen- 
generator erkannt. Eine Wiederholung von Dezimalstellen ist dabei ausge- 
schlossen, da man weiß, daß eine transzendente Zahl (d.h. nicht Nullstelle 
eines Polynoms) ist und somit eine nichtperiodische Entwicklung hat. Der 
Transzendenzbeweis erfolgte 1882 durch Lindemann. 


Eine genauere Approximation, als die oben gegebenen Werte darstellen, durch 
eine Bruchzahl kann man mit Programm 10 erhalten. Gibt man dort auf 
8 Dezimalstellen ein, so erhält man 


251699445/80115557 


Die Computerrechnungen stützen sich alle auf folgende Darstellungen von 7: 
Formel von Machin: 
a = 16 arctan (1/5) — 4 arctan (1/239) 


n = 16 * arctan (1/5) — 4 * arctan (1/70) + 4 arctan (1/99) 
oder die Formel von Gauß 
n = 48 arctan (1/18) + 32 arctan (1/57) — 20 arctan (1/239). 


Die Reihenentwicklung des arctan 


-„_1,3,15_1,7,19_, + 
arctanx=x 3x 5x treten. 
gültig für Ixl < 1 konvergiert umso schneller, je kleiner x ist. Insofern wäre die 
Gaußsche Formel vorzuziehen, jedoch muß die Reihenentwicklung dreimal 
ausgewertet werden. Daher benützt man meist die Formel von Machin. 


Die ungeraden Potenzen der Reihenentwicklung erhält man aus der 1. Potenz 
durch fortgesetztes Multiplizieren mit dem Quadrat von x bzw. durch Divi- 
sion durch das Quadrat des Nenners, 

Diese Division einer langen Zahl mit den Stellen U(I) durch die kurze Zahl D 
kann mit Hilfe von folgendem Programmstück durchgeführt werden: 


108 K=0 

118 FOR I=1 TON 

120 U=U() +10 *K 

139 U(I) = INT(U/D) 

1489 K=U - INT(U/D) *D 
159 NEXT I 


Dabei ist N die Stellenzahl und K der jeweilige Übertrag. Bemerkenswert ist 
daß die Zahl m auch zweimal als Ergebnis eines Zufallsversuchs auftritt. 
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Das eine ist das berühmte Nadelproblem von Buffon. Dabei wird eine Nadel 
der Länge I zufällig auf ein Parallelgitter vom Abstand d geworfen (I<d). 
Die Nadel trifft eine Parallele mit der Wahrscheinlichkeit 


2l/nd. 


Ein zweites Ergebnis stammt von R. Chartres: Er zeigte 1904, daß zwei zufäl- 
lig ausgewählte natürliche Zahlen mit der Wahrscheinlichkeit 


6 


m? 


teilerfremd sind. 


Zum folgenden Programm 


Es wird darin die Zahl 7 nach der Formel von Machin berechnet. Die Division 
durch 25 bzw. durch 2392 werden mit Hilfe des angegebenen Verfahrens 
durchgeführt. Der sich ergebende Reihenterm wird entsprechend seinem Vor- 
zeichen addiert bzw. subtrahiert. 

Die Ausgabe erfolgt in Dreierblöcken. 


Mit Hilfe des hier gegebenen Programms kann 7 auf beliebig viele Stellen be- 
rechnet werden. Dazu muß die Variable N in Zeile 119 auf 1/3 der gewünsch- 
ten Stellenzahl gesetzt werden. Es ist ratsam N noch aufzurunden, damit 
nicht das Ergebnis durch einen Rundungsfehler verfälscht wird. 


Die Laufzeit des Programms beträgt mehrere Stunden, so daß Sie entspre- 
chende Geduld aufbringen sollten. Während des Programmlaufes erscheint auf 
dem Bildschirm eine Countdown-Anzeige, die über den Fortschritt der Rech- 
nung informiert. 


Die vom Programm berechneten Stellen können dem Programmausdruck ent- 
nommen werden. 
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419 
420 


REM PI AUF 1288 STELLEN 


PRINT " " 
PRINT " PI AUF 1088 STELLEN " 
PRINT " " 
NORMAL 
PRINT : PRINT 
PRINT "AUSGABE AUF <B>ILDSCHIRM" 
INPUT " ODER <D>RUCKER 7"; P$ 
PRINT : PRINT 
PRINT "BERECHNUNG LAEUFT" 
N = 334 
IZ=eN#k3 
T = 1806: DIM RY(N,2),K(2) 
Kil) = 25:K(2) = 239% 2 
FOR P = 1 TO 2: GOSUB 489 
GOSUB 770:A = 1:D1 = E: GOSUB 428 
IF V > 8 THEN GOSUB 528: G0TO 219 
GOSUB 40R 
E=E+2:V= -V 
A= 8:Di = KiP): GOSUB 429 
HTAB 28: VTAB 18: PRINT zZ" 9; 
27 = 212-1 
IF Z THEN 188 
NEXT P 
IF P$ = "D" THEN PR# I: GOTO 329 
PRINT : PRINT CHR$ (75 CHR$ (7) 
PRINT "RETURN FUER AUSGABE"; 
INPUT W$ 
HOME 
PRINT : PRINT "PI = "3 
PRINT STR$ (RX(B, 2)", " 
FOR I = 1 TON 
IF(I-D /9= INTÜI-1 7/9 THEN 
X$ = STR$ (RX(1,2)) 
PRINT RIGHTS ("@0" + X$,3)" "3 
NEXT I 
PRINT 
PRR 8 
END 
R=0:1=9 


50 


PRINT 


FORI=6TON 
=RX*T + RA(I,A) 
D/D 
D-Di* INT (0) 
ZORQ 
RY(I,A) = Q 
NEXT I 
REIN,A) =Q= „5 
RETURN 
C=9 
FORI=NTO@STEP -1 
S = RAtI,2 + RA(L,U) + C 
C=9 
IFES> =TTHENS=S-Til= 
RY1,D) =8 
NEXT I 
RETURN 
L=9 
FORI=NTO@STEP - 1 
D2 = RY(1,2) - RAL,D) -L 
L=9 
IF D2 < @ THEN D2 = D2 + Til 
R4(1,2) = D2 
NEXT I 
RETURN 
FOR I=8T0ON 
R%(I,G) = 8: IFP = 1 THEN RAX(L,2 = 9 
NEXT I 
RB, = 16 /P* 2 
DI = 234 * P - 29 
A=8 
GOSUB 429 
E=1:V=3-2#P 
RETURN 
FOR I=8T0ON 
RY{I,1) = RA(1,9) 
NEXT I 
RETURN 


N DO 


) 
>= 


u 
_. 
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8. e AUF 1988 STELLEN 


Mit der Zahl ist die Eulersche Zahl 
e=2.7182818... 

über die Gleichung 
eTI= 1 

verbunden, dabei ist i= V—1. 

e ist insbesondere als Basis der natürlichen Logarithmen von Bedeutung 
y-Inx> x>=eY 

Ebenso stellt e den Grenzwert 


e= lim (1+1)x 
X > 0 X 


dar. Mit diesem Grenzwert ist auch die Zinzeszinsformel verknüpft: Wird 
1 DM nicht jährlich oder monatlich, sondern kontinuierlich verzinst, so 
wächst der Zinseszins bei 100% nicht ins Unbeschränkte, sondern nähert sich 
dem Wert e DM. Die Reihenentwicklung der Exponentialfunktion liefert eine 
für numerische Zwecke brauchbare Formel für e 


1,1,1,1,1,1, 


e=-—- + WT 51 61 .... 


++ 


11 21 31 
dabei stellt n! die Fakultätsfunktion dar. 


Erwähnenswert ist, daß auch e als Ergebnis eines Zufallsexperiments auftritt. 
Werden n Dinge zufällig angeordnet, so ist mit der Wahrscheinlichkeit 
1/e = 0.3678794. 


kein Gegenstand an seinem Ort. Solche Zufallsexperimente können mit Hilfe 
von Programm 12 erhalten werden. Der Versuch ist besser bekannt als das 
Problem der ""vertauschten Briefe’’. Legt man n Briefe zufällig in n vorbereite- 
te Umschläge, so liegt mit der Wahrscheinlichkeit 


1 — 1/e = 0.6321205 


mindestens ein Brief im richtigen Umschlag. Im Band 1 der vorliegenden Pro- 
grammsammlung findet sich für n = 5 ein Grafikprogramm. 
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Zum folgenden Programm 


Das Programm berechnet die Dezimalbruchentwicklung von e® mit Hilfe der 
angegebenen Reihenentwicklung. Die Kehrwerte der benötigten Fakultäten 
werden durch wiederholte Divisionen berechnet und fortlaufend aufsummaiert. 


Das hier angegebene Programm läßt sich auf beliebige Stellenzahl ausdehnen. 
Die Variable N in Zeile 180 muß dazu so geändert werden, daß die Dezimal- 
stellen von 1/NI! bei der gewünschten Genauigkeit vernachlässigt werden kön- 
nen. 
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REM E AUF 1898 STELLEN 


HOME =: INVERSE 
PRINT " 


PRINT " E AUF 1898 STELLEN " 


PRINT " 
NORMAL 
FRINT 


PRINT "DRUCKERAUSGABE <J/N> " 


INPUT PR$ 


PRINT : PRINT "FELD ANLEGEN" 


DIM A(202) 

FOR I = 1 TO 282 
aD) =9 

NEXT I 


PRINT 


PRINT "BERECHNUNG LAEUFT...": PRINT 


B= 1ESAB) = 1 
FORN = 458 TO 1 STEP -1 
HTAB 25: VTAB 9: PRINT Ns" 
FOR I = 8 TOD 2@1 
@ = INT (AT /MN 
R=AMTMD -UHN 
AD =Q 
AI+1) =SAUlI+1) +BH#R 
NEXT I 
AB) = AM +1 
NEXT N 
REM RUNDUNG 
PRINT : PRINT 


PRINT "RUNDUNG AUSFUEHREN. .." 


A(291) 
FOR I = 288 TO 1 STEP -1 

U= INT (AD /B 

A =ATD -B*U 

AII-1I =AUL- I +U 
NEXT I 


REM AUSGABE 
IF PR$ = "J" THEN PR# 1: 


55 


Atzei) + INT (A282) /B + .5) 


GOTO 472 


438 PRINT CHR$ (7)5 CHR$ (MD 

448 PRINT "WEITER MIT AUSGABE <RETURN> "5 
450 INPUT T$ 

468 HOME 

478 FOR I = 8 TO 288 

ABB At = STR$ (AtT)) 

498 IF I = 8 THEN PRINT A$ss".": GODTD 529 
508 PRINT RIGHT$ ("00980" + A$,D)5" "; 
S18 IF1/%&= INT (17/6) THEN PRINT 
5280 NEXT I 

558 PRINT 

349 : 

358 IF PR$ = "J" THEN PR# 9 

369 END 


96 


18284 
77572 
40766 
64274 
35764 
07381 
27531 
21548 
B&688 
34424 
27618 
49338 
870971 
01416 
21112 
34499 
68892 
92295 
31825 
98294 
81978 
32823 
98194 
12589 
85193 
849798 
84445 
70419 
99418 
62531 
642453 
72085 
864196 
30354 


39845 
47093 
32353 
27466 
29043 
32328 
95251 
87149 
82264 
37197 
M1-1-177-) 
25569 
27443 
72936 
32389 
699767 
30098 
76351 
28869 
88793 
68416 
76464 
35815 
961B1 
45807 
92085 
63443 
78623 
49146 
32996 
78149 
18383 
87396 


23536 
697995 
34759 
39193 
37290 
62794 
81921 
75488 
=] 02 8-) 
33997 
26153 
29768 
74704 
819792 
78442 
94686 
79512 
48228 
39849 
32036 
14039 
80429 
30175 
88159 
2738& 
80582 
24296 
29998 
31499 
18569 
39271 
75851 
96552 


02874 
95749 
45713 
20938 
055342 
34987 
15738 
41675 
84774 
77249 
13845 
67371 
72506 
So151 
0569 
44549 
77361 
82498 
64651 
25074 
72198 
35118 
&7173 
Ssg416 
&7385 
37492 
84875 
21699 
34317 
28887 
45635 
81157 
12671 
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71352 
669767 
82178 
39921 
95269 
63233 
34187 
99244 
11853 
72069 
830090 
13208 
96977 
08457 
33696 
25787 
78215 
95193 
25822 
43117 
37679 
02528 
61332 
98351 
87422 
79618 
60233 
70235 
38143 
27016 
AIDbl 
47794 
34688 


9. LINEARE DIOPHANTISCHE GLEICHUNG 


Ein Standardverfahren zur Lösung einer linearen Diophantischen Gleichung 
axtby=c 

gleichbedeutend mit der Kongruenz 
ax=cmody 


ist der erweiterte Euklidsche Algorithmus. Mit seiner Hilfe kann der größte 
gemeinsame Teiler 


ggT (a,b) 
als Linearkombination von a und b bestimmt werden. 


Dies soll an einem Beispiel für a= 286 und b = 121 gezeigt werden: 
262 : 121 = 2 Rest 44 > 286 =2 - 121 + 44 
121: 44 =2 Rest 33> 121 =2- 44 +33 
44: 33=1 Rest 11> 44=1-33+11 
33: 11=3 > 33=3-11 
Somit gilt 
ggT (286, 121) = 11 
Durch Rückwärtsrechnen folgt 
11 = 44 — 33 
= 44 — (121 — 2.44) 
=3.44 — 121 
=3.(286 — 2.121) — 121 
=3.286 — 7. 121 
Dies liefert die gesuchte Linearkombination 
3:286 — 7- 121 = 11 
allgemein 
ax +by =ggT (a,b) 
Sind xo, Yo die Koeffizienten der Linearkombination, so hat die Diophanti- 
sche Gleichung 
axt+tby=c 


die allgemeine Lösung: 


x=CXo t 


+. _. 
ggT (a,b) 


=c b «t; t ganzzahlig 
vYoT ggtlabl 
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Diese Lösung existiert nur, wenn c ein Vielfaches von ggT (a,b) ist; andern- 
falls ist die Gleichung nicht lösbar (vgl. [7]). 


Zum folgenden Programm 


Das Programm berechnet mit Hilfe des erweiterten Euklidschen Algorithmus 
nach dem oben angegebenen Verfahren eine spezielle Lösung der Gleichung 


ax +by =ggT (a,b). 
Damit wird die allgemeine Lösung bestimmt. 


Das im Programm angegebene Verfahren kann prinzipiell auch für Diophanti- 
sche Gleichungen der Form 


ax —by=c 


{b > 0) angewandt werden. In diesem Fall ergibt sich für den ggT ein negati- 
ver Wert, im Gegensatz zur üblichen Definition. Dazu müssen die Zeilen 120 
und 130 gelöscht werden; zwei aufeinander stoßende negative Vorzeichen 
sind als positiv zu lesen. 


Als Beispiel wird die Diophantische Gleichung 
17x +6y=5 

gelöst. Das Programm liefert die Parameterlösung 
x=-5 -6t 
y = 15 +17t 


dabei durchläuft t alle ganzen Zahlen. Für jeden Wert von t ergibt sich wieder 
eine spezielle Lösung, z.B. für 
t=0:x=-5,y=15 
t=1:x=-11,y=22 
=-1:x=1,y=— 
usw. Zahlreiche Diophantische Gleichungen findet man in eingekleideter 
Form in den historischen Aufgaben des Anhanggs. 
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REM LINEARE DIOPHANTISCHE GLEICHUNG 
HOME : INVERSE 
PRINT " " 
PRINT " DIOPHANTISCHE GLEICHUNG " 
PRINT " n 
PRINT " AxX + BiY = C n 
PRINT " " 
NORMAL 
PRINT : PRINT 
INPUT "EINGABE A,B,C ? "5A,B,C: GOTO 149 
IFA>M@ AND B > O0 THEN 170 
PRINT "A UND B MUESSEN POSITIV SEIN": GOTO 119 


PRINT "-------------------- " 
PRINT 
AB = Aral = B 
X = 1ıYi = 1 
YO = 0:1 =9 


REM ERWEITETER EUKLIDSCHER ALBORITHMUS 
Q@ = INT (Ad / AU 

IF AQ / Al = Q THEN 328 
A2=AM -Ux* MA 


2 =2X8 -0Q x Xi 
Y2=Y-o#*Yl 

A@ = Al:Al = A2 
XD = XliXi = X2 

Y@ = Yi:Yl = V2 
GOTO 229 


PRINT "BBT("3A5","5B5")="5A2 

PRINT : PRINT 

PRINT "LOESUNG DER GLEICHUNG A#X+B#Y=GBT (A,B)" 
PRINT 

PRINT Az "x"; 

IF X2 < 8 THEN PRINT "("5X27")"5: GOTO 4990 
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398 
a12]7) 
419 
420 
438 
440 
450 
469 
470 
489 
490 
RB 
319 
320 
3sß 
J4B 
ei] ) 


PRINT X23 
IF B< 8 THEN PRINT Bi: GOTO 428 
PRINT "+"5Bj 


PRINT 


nn; 


IF Y2 < 8 THEN PRINT "("3Y25")"32 GOTO 458 
PRINT Y25 


PRINT 


"= A2 


IFC/A2< > INT {C / A2) THEN GOTD 549 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
END 

PRINT 
END 


: PRINT : PRINT "ALLG. LOESUNG DER DIOPH. 
"X=":C % “2 / A235 "-TR" Br" 2 
"ans # V2 / A2p"+TR"SA; 


" T GSANZZAHLIG" 


. PRINT "KEINE LOESEUNG" 
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GLEICHG, 2" 


DIOPHANTISCHE GBLEICHUNG 


AX + BrY=C 


EINGABE A,B,C ? 17,6,5 


G8T (17,6) =1 


LOESUNG DER BLEICHUNG A#X+B#Y=GGT(A,B) 


17#(-1) +6#3=1 


ALLG. LOESUNG DER DIOPH. GLEICHG.: 


1=-5-T%#6 Y=15+T#17 T GANZZAHLIG 
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a5 0 1 5 5 .. 
fl. 39 3 I>s haupegücnngerwindes erflöjare. 
196875 
413]4375 Anden 
20 6>18>5 
434|1093>5 Diiteen 
21>05468>5 
ass|814843>5 Vierdten 
227907421875 
a>slcossse93>, Jünfften 
23 0392>92068>5 
g02|5358652343>5 Sechften 
25 126>932C1718>75 
&2>|66265840C093>5' Sibendeit 
2638313292489408>% 
554|04579 1420808437 5 Achtete 
2>>022805>10449218>?5 
83 174808090 1943359375 Fleiindee 
209 08>404000959> 1067968>5 
fl.610]83548504 1540527343 > 5 3cherd 
BE 6168>88033232421 8>5000 
+ 20|616409969>2656250000. 
>2 Dierzofferage zjar phanptgürzinsund 
sinszinsı 32 fr 2ßı2 +.Bringesing birzinssing 
12flafr2 +. Darnachdie 250 FR erage 3 jar 
Baupeg: ins vn sinssink 280 Fraß> +1.0nd 
st balger sine des vierdräjare > FR ıf D 15 
Y 


“Aus dem Exempel-Büchlin von Christoph Rudolff, Augsburg 
1530. In diesem Buch wurden erstmals Dezimalbrüche 
systematisch genutzt” 
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10. KETTENBRUCH / BRUCHAPPROXIMATION 


Sucht man zu einem vorgegebenen Dezimalbruch einen Quotienten gleichen 
Werts — z.B. ein Windungszahlverhältnis beim Transformator — so kann dieser 
mit Hilfe eines Kettenbruchs bestimmt werden. 


Ein Kettenbruch hat folgende Form 


Solche Kettenbrüche können für rationale, d.h. für Bruchzahlen durch den 
Euklidschen Algorithmus berechnet werden. 


Für 67/29 ergibt sich z.B. 


67 :29=2 Rest 9 
29: 9=3Rest?2 


9: 2=4Resti 
2: 1=2 
Somit folgt 

67 9 _ 1 

2 2'2°'5 
9 

29 2 1 

2 3tg=3tg 
2 

9 1 

a ttz 


Setzt man diese Gleichungen rückwärts ineinander, so ergibt sich der Ketten- 
bruch 

67 1 

—= + — 

29 2 3+ 1 
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Der Euklidsche Algorithmus zeigt, daß der Kettenbruch einer rationalen Zahl 
abbricht, jedoch nicht der einer irrationalen Zahl. Ein Beweis dazu findet sich 
in [14]. Kettenbrüche können somit durch folgendes Programm berechnet 
werden 

188 INPUT "DEZIMALBRUCH”; X 

118 PRINT INT (X) 

128 IF X=INT(X) THEN 150 

138 X=1/(X-INT(X)) 

148 GOTO 118 

158 END 


jedoch muß dabei noch die beschränkte Genauigkeit des Rechners beachtet 
werden. 


Zum folgenden Programm 


Das folgende Programm berechnet zu einem eingegebenen Dezimalbruch maxi- 
mal 8 Teilnenner nach dem oben angegebenen Programmstück. Dabei wird 
versucht Rechnerungenauigkeiten, so weit wie möglich, abzufangen. 


Der so berechnete Kettenbruch wird auf dem Drucker ausgegeben und die zu- 
gehörige rationale Zahl berechnet. 


Als Programmbeispiel wurde 
na = 3.14159265 

eingegeben. Man erhält die Bruchapproximation 
251690445 
80115557 


Der ausgegebene Kettenbruch kann dem Programmausdruck entnommen wer- 
den. 
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REM KETTENBRUCH/BRUCHAPPROXIMATION 


HOME : INVERSE 
PRINT " " 
PRINT " KETTENBRUCH/BRUCHAPPROXIMATION " 
PRINT " " 
NORMAL 
PRINT : PRINT 

INPUT "DRUCKER EINGESCHALTET ? "1A$ 

IF LEFT$ (A$,1) < > "J" THEN SQ 
PRINT : PRINT 

INPUT "DEZIMALZAHL ? "5Y 

DIM R(12) 

PR# 1 
K=B 

FOR I=1T0OK 
RD = INTV+IE-95 

IF ABS (RD - YV < IE - 5 THEN 242 
Yz1/XY-RD) 

NEXT I 

GOTD 259 
K=1I 

PRINT Rt1);" + 1" 

FORI=2TOK-1 

PRINT TABU I *& - Na" " 

PRINT TAB\ I *& - SJ);R(TD);" + 1" 
NEXT I 

PRINT TABU K #6 - My"! 

PRINT TABUK * & - S)IR(kK) 
z = R(kK) 
N= 1 

FOR J=K- 1 TO 1STEP - 1 


I = 
z=R(U) #*Z +N 
N=] 
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NEXT J 


PRINT " "12 

PRINT "=--------- 
PRINT Z / N 

PRINT " "N 

PRINT : PRINT 

PR# 8 

END 


68 


251690445 


890115557 
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11. PERMUTATIONEN 


Da endliche n—elementige Mengen stets auf die Menge 
11,2,3,...n| 


abgebildet werden können, arbeiten die folgenden Kombinatorik-Programme 
stets mit dieser Anfangsmenge der natürlichen Zahlen. 

Jede Anordnung der Menge 1,2,3,.. .,n, bei der jedes Element genau einmal 
vorkommt, heißt Permutation. 


Permutationen treten vielfach in der Kombinatorik und Wahrscheinlichkeits- 
rechnung auf. Man benötigt sie insbesondere auch in Operations Research, 
wenn man z.B. auf der Suche nach der kürzesten Rundreise alle Anordnungen 
der zu besuchenden Städte durchmustern muß. 


Es gibt 
n! 


Permutationen einer n—elementigen Menge; dabei ist n! die in Programm 6 be- 
handelte Fakultätsfunktion. 


Die Fakultäten der Zahlen wachsen sehr schnell, so gilt 


6! = 720, 71=5049, 8! = 40320 
9! = 362880, 101 = 3628800 
11! = 39916800 


Somit gibt es für 10 Leute bei 10 Sitzplätzen 
3628800 


verschiedene Sitzordnungen. Spielt in einer Fußballmannschaft jeder auf jeder 
Position, so gibt es 


39916809 
verschiedene Mannschaftsaufstellungen. 


Zu folgendem Programm 


Das folgende Programm erzeugt nach einem Algorithmus von Ord-Smith die 
Permütation einer n—elementigen Menge in lexikographischer Anordnung. 


Dabei heißen zwei Zahlenfolgen a; und b; lexikographisch geordnet, wenn für 
die erste Stelle mit 


atbi 


gilt 
ash 


Der Programmausdruck zeigt die Permutation einer 4—elementigen Menge. 
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id REM PERMUTATIONEN 

280 HOME : INVERSE 

3s® PRINT " " 

48 PRINT " PERMUTATIONEN " 

8 PRINT " " 

68 NORMAL © PRINT : PRINT 

70: 

80 INPUT "WIEVIELE ELEMENTE ? "IN 
98 DIM X(N),Y(N 

180 PRINT 

118 INPUT "DRUCKERAUSGABE ? ";PR$ 
128 PRINT 

130 : 

148 IF PR$ = "J" THEN PR& 1 

158 REM STARTWERTE 

160% FORI=1T0ON 

178XD = 1 

188 NEXT I 

IBF=-9P=9 


218 FOR I=1TON 
228 PRINT X(TD}3 

238 NEXT I: PRINT 
2ZAöP=P+|1 

258 GOSUB 338 

268 IF F = 1 G0OTD 218 


288 PRINT : PRINT 

298 PRINT P5" PERMUTATIDNEN" 
Ss22 PR# 9 

Siß END 

320 : 

S5@ REM UNTERPROGRAMM 

Sad NI=N- 1 

s5® IFF = 1 THEN 498 


SB F= 1 
378 FORK=1TO N 
Ss88 Y(K) = N 
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578 
498 
410 
420 
430 
449 
4582 
460 
470 
480 
490 
pa] 717) 
319 
320 
339 
348 
3 
JbR 
378 
382 
378 
[#171 
619 
6209 
63% 
&40 
658 


NEXT K 

IF Y{iNIDD <> N THEN 459 
Y{NDD = NI:L = X(N 

XN = XtNDDEXND = L 

GOTD 658 

FOR J=1 TO NM 
K=-N-J 

IF Y(K) < > K THEN 548 
Y{iK) = N 

NEXT J 


GOTO 589 


ı ı SZ N 
< 
ı uU = 


zZ Rr 
+ 
> 


= X(M:EX(M) = Xi(K) 
ı=L 


KRZIxT" ZER << Ze 
I 
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PERMUTATIONEN 


WIEVIELE ELEMENTE ? 4 
DRUCKERAUSGABE ? J 


1234 
1243 
1324 
1342 
1423 
1432 
2134 


24 PERMUTATIONEN 
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12. ZUFALLSPERMUTATION 


Für viele Zwecke benötigt man zufällige Permutationen, z.B. für die bei Pro- 
gramm B erwähnten Monte-Carlo-Simulation der vertauschten Briefe. 


Als Anfangswert wird die Permutation P(l) der Zahlen 1,.., n auf die natürli- 
che Reihenfolge gesetzt. 


Sodann werden n—1 ganzzahlige Zufallszahlen ausgelost und die Permutation 
mit entsprechendem Index ausgetauscht. Dies wird durch folgendes Pro- 
grammstück durchgeführt. 


268 FOR J=NTO2STEP —1 

278 R= INT (J* RND(N)) +1 

280 H=P(J) : P(J) =P{R) : P(R)=H 
290 NEXT J 


Die so erhaltene Zufallspermutation wird ausgedruckt und bei Bedarf eine 
neue erzeugt. 


Zum folgenden Programm 


Nach Abfrage, wieviele Zahlen die Permutation umfassen soll und wieviele 
Permutationen erwünscht sind, werden nach dem angegebenen Verfahren mit 
Hilfe des eingebauten Zufallszahlengenerators die benötigten zufälligen Permu- 
tationen erzeugt. 

Der Programmausdruck zeigt 20 Zufallspermutationen der Zahlen 1 bis 9. 
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REM ZUFALLSPERMUTATIONEN 

HOME : INVERSE 

PRINT " " 

PRINT " ZUFALLSPERMUTATIONEN " 

PRINT " " 

NORMAL = PRINT = PRINT 

INPUT "WIEVIELE ZAHLEN ? "5N 

DIM P(N) 

PRINT : INPUT "WIEVIELE ZUFALLSPERMUTATIONEN ? "5M 
PRINT 


I= 1: REM ZAEHLER 


REM ANFANGSWERTE 

FOR J=1TDN 
PtJ) =J 

NEXT J 

REM AUSLOSEN DER REIHENFOLGE 
FOR J=NTO2STEP -1 
R= INT (J* RND (1i)) +1 

H = PiW:P(J) = P{R):P{R) = H 
NEXT J 

REM AUSGABE 

FOR J=1TON 

PRINT P(J)5 

NEXT J: PRINT 


+ 1 
“© = M THEN 218 
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ZUFALLSPERMUTATTONEN 


WIEVIELE ZAHLEN ? 9 
WIEVIELE ZUFALLSPERMUTATIONEN ? 28 


6284751539 
813675249 
759421883 
69781545372 
458917326 
362389147 
273158694 
486173527 
217385746 
432571889 
163549782 
765921834 
398521457 
145867923 
243876951 
718542936 
718632457 
864925157 
314596278 
237641895 
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13. KOMBINATIONEN 


Ordnet man nur k Elemente einer n-elementigen Menge ohne Wiederholung 
an, so erhält man eine k-Kombination oder auch Kombination k.-ter Klasse. 


Die Anzahl der k-Kombinationen einer n-elementigen Menge schreibt man 


n 
k 


und bezeichnet dies als Binomialkoeffizient. Diese tragen ihren Namen daher, 
weil sie als Koeffizienten in den binomischen Formeln, z.B. 


= a9 + 3a2b + 3ab2 + b9 


auftreten. 


Die Binomialkoeffizienten | "] i= 0,1,2,...,n finden sich jeweils in der 


(n+1)-ten Zeile des Pascalschen Dreiecks: 


Die Binomialkoeffizienten können auch über die Fakultätsfunktion berechnet 


werden 
n _ n! 
(%) kl (nk)! 


So kann beim Schafkopf (24 Karten, 3 Spieler) jeder Spieler 


2 = 735471 


Blätter haben. Beim Lottospiel gibt es 
2) = 13983816 


verschiedene Lottotips, beim neuen Mittwochslotto entsprechend 


®) = 12620256 
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Zum folgenden Programm 
Das Programm erzeugt nach einem Algorithmus von Kurtzberg alle k-Kombi- 
nationen (oder k-ten Klasse) einer n-elementigen Menge. 
Der Programmausdruck zeigt die 
4-Kombinationen 


einer 6-elementigen Menge. 


| >= >< > 
Eyn 
RHmD tvolgegründre 
Bi | onderweyfüng aller Rauffimanß Rcch- 
gung in Dreyen Büchernymic (chönen Re 
BR geln vn fragffucfen Benriffen. Sunder- 
N ]ilıch was forrl unnd Behendigfaic in der 
| Weliche Practicavii Tolleen geßrauchr | |P 


x WINDE / Des gleychen fürmalß weider in 
Teugfcher noch in Weclfcher fpiach nic 


_ gedeiiche. Durch Petrum Apiand 

a von Zeyinich/d Altronomer 
A x, 30 Ingolftar Ördinas zu 
% AU, verferriger, KT 


"Titelblatt der Arithmetik von Petrus Api 
Pplanus, Ingolstadt 15 
mit der ersten gedruckten Darstellung des Pascalschen Dreiecken 
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REM KOMBINATIONEN 

HOME : INVERSE 

PRINT " " 

PRINT " KOMBINATIONEN " 

PRINT " " 

PRINT = NORMAL 

INPUT "WIEVIELE ELEMENTE ? "5N 

PRINT : INPUT "ZU WELCHER KLASSE ? ";M 
DIM X(M) 


REM STARTWERTE 
L=l 

FOR IT=1TOM 
XD) =N 

NEXT I 

REM M AUS N 
P=|1 

IF M = 9 THEN 26% 
FOR I =1 TOM 
P=P*(N-I+1i) /I 
NEXT I 

PRINT 


GOSUB 378 

FOR I=S1I1 TOM 

PRINT XD; 

NEXT I: PRINT 
L=L+rli 

IFL< = P THEN 269 
PRINT 

PRINT P5" KOMBINATIONEN" 
END 


REM UNTERPROGRAMM 
FORK=1T0OM 
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JeaKk-1l 
A=M-J 
B=N-J 
IF XtA) < B THEN 599 
NEXT K 


FORK=ATOM 
B=B+ti 

X(K) = B 

NEXT K 

RETURN 
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KOMBINATIONEN 


WIEVIELE ELEMENTE ? 6 
ZU WELCHER KLASSE ? 4 


1234 
1235 
1256 
1245 
1246 
1256 
1545 
1346 
1356 
1456 
2545 
2346 
2356 
2456 
3456 


15 KOMBINATIONEN 


ULOADKOMBINATION,D2 
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14. VARIATIONEN 


Ordnet man k Elemente einer n-elementigen Menge unter Berücksichtigung 
der Reihenfolge an, so erhält man eine k-Variation oder Variation zu k-ten 
Klasse. 
Es gibt 
n! 
(n—k)I 


k-Variationen einer n-elementigen Menge. 


Wollen von 20 Personen 8 Leute an einem bestimmten Tisch sitzen, so gibt 
es dort 


201 _ 5979110400 


121 

verschiedene Sitzordnungen. Aus einem Gremium von 30 Personen gibt es 
20 _ 24360 
27! 


verschiedene Möglichkeiten einen Präsidenten, dessen 1. und 2. Stellvertreter 
auszuwählen. 


Zum folgenden Programm 


Das Programm erzeugt alle k-Variationen oder Variationen zur k.-ten Klasse 
einer n-elementigen Menge. 


Der Programmausdruck zeigt alle 3 Variationen einer 4-elementigen Menge. 
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1d REM VARTATIONEN 
28 HOME © INVERSE 


39 PRINT " ö 
48 PRINT " VARIATIONEN " 
59 PRINT " " 
&8 NORMAL : PRINT 

70: 


88 INPUT "ZAHL DER ELEMENTE ? ";N 
978 PRINT : INPUT "ZU WELCHER KLASSE 
188 PRINT 

118. 

128 REM STARTWERTE 

138 FOR I =1 TON 

148 XD) = I 

150 NEXT I 

188 M=K 

170 IFK=NTHENM=M- 1 

188 FOR I1=1 TOM 


248 FOR J=1TOK 
258 PRINT X(W; 
268 NEXT J: PRINT 
27 UV=V+1 


290 IF Y(D > = N THEN 348 
388 Y{ID) = Y{TD) + 1:H = X(I) 
Si XD = XIY{DIEX(YCD) = H 
320 GOTD 239 

338 : 

348 H = X(D:X(D) = X(Y(D)) 

3s58 XY{D) = H!Y(D = Y(T) - { 
360 IF Y(TD) > I THEN 349 
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"K 


Ss978 IF I > 98 THEN 299 

498 : 

418 PRINT 

428 PRINT v5" VARIATIONEN" 
458 END 
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VARIATIONEN 


ZAHL DER ELEMENTE ? 4 
ZU WELCHER KLASSE ? 53 


123 
124 
132 
134 
142 
143 
213 
214 
231 
234 
241 
243 
312 
314 
S21 
324 


24 VARIATIONEN 
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15. 91-TUPEL 


Wählt man k Elemente aus einer n-elementigen Menge aus und läßt dabei be- 
liebige Wiederholung zu, so erhält man die 


k-tupel 
der Menge. Im Fall von k=2 spricht man auch von Paaren, bei k=3 von Tripel. 


Es gibt insgesamt 
nK 


k-Tupel einer n-elementigen Menge. 


Besondere Bedeutung haben die 8-Tupel der Binärzahlen ® und 1, Byte ge- 
nannt, für Computer. Alle ganzen Zahlen werden mit Hilfe zweier Bytes dar- 
gestellt. Da das erste Byte das Vorzeichen angibt, können somit INTEGER- 
Zahlen mit 15 Bit dargestellt werden. Es sind dies die Zahl ® bis 32767. 


Weitere Anwendungen finden die ®91-Tupel, da mit ihrer Hilfe alle Teilmengen 
einer gegebenen Menge erzeugt werden können: Ist die i.-te Stelle des 81- 
Tupel Eins, so gehört das Element i zur entsprechenden Teilmenge, sonst 
nicht. 


Somit hat eine n-elementige Menge soviele Teilmengen, wie der B1-Tupel der 
Länge n gibt, nämlich 


an 


Wichtig sind 91-Tupel auch in der Boole’schen Algebra und Aussagenlogik, da 
mit ihrer Hilfe alle Wahrheitswerte von Aussageverknüpfungen ermittelt wer- 
den können (siehe Logikprogramme im Band 1 der vorliegenden Programm- 
sammlung). Allerdings wird im CBM-BASIC *'wahr” nicht als 1, sondern als 
—1 codiert; dies kann aber im Programm durch eine Vorzeichenänderung er- 
reicht werden. 


Zum folgenden Programm 


Das Programm erzeugt nach einem Algorithmus von Berztiss alle ®1-Tupel der 
gewünschten Länge. 

Läßt man in Zeile 258 statt X(K) das Element A(K) einer Menge — falls 
X(K) = 1 — ausdrucken, so erhält man alle Teilmengen dieser Menge. 


Der Programmausdruck zeigt alle @1-Tupel der Länge 5. 
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8 
9 
30 
33 
T 
3 
4 
Ss 


. Iy IN ee... 


“Entwurf von Leibniz für eine Gedenkmünze, die das von ihm gefundene 
Binärsystem darstellt (1697)'' 
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REM 81-N-TUPEL 
HOME : INVERSE 
PRINT " " 
PRINT " @1-N-TUPEL " 
PRINT " " 
NORMAL : PRINT 
INPUT "LAENGE DER B1-TUPEL ? "IN 
DIM X(N) 
PRINT 
REM STARTWERTE 
J=2”N 
L=1 
FOR K 
X(K) = 
NEXT K 
FORK=1 TON 
PRINT X(K)5 
NEXT Kı PRINT 
L=L+1 
IF L > J THEN 278 
GOSUB 518 
GOTO 1982 


=1T0ON 
Ö 


PRINT 

PRINT Js" "5N3 "-TUPEL" 
END 

REM UNTERPROGRAMM 
I=1 

IF XD < > 98 THEN 356% 
XD = 1 

GOTD 3978 
xD = 9 
I=1Ir+i1l 

IFI< = N THEN 339 
RETURN 
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81-N-TUPEL 


LAENGE DER B1-TUPEL ? 5 


22008 
10829 
819089 
11998 
22108 
10198 
21108 
11108 
229010 
12219 
21018 
11219 
20118 
12119 
21110 
11118 
20201 
19991 
21201 
11291 
28101 
12101 
B11B1 
11181 
DOD11 
10011 
vıali 
11211 
BB11ll 
18111 
Billl 
11111 


32 5-TUPEL 


16. PARTITIONEN 


Unter einer Partition der Zahl n versteht man die Zerlegung von n in Sum- 
manden, wobei es auf die Anordnung nicht ankommt. So hat die Zahl 6 
folgende 11 Partitionen: 


6=6 

5+1 

4+2 

4+1+1 

3+3 

3+2+1 
3+1+1+1 
2+2+2 
2+2+1+1 
2+1+1i+r1+1 
Ii+t1+1+r1+r1+1 


Solche Partitionen werden benötigt bei Geldwechsel- und Frankatur-Proble- 
men. 


Zum folgenden Programm 


Das Programm erzeugt sämtliche Partitionen einer Zahl nach einem Algorith- 
mus von McKay. 


Der Programmausdruck zeigt die 42 Partitionen der Zahl 10. 
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REM PARTITIONEN 

HOME : INVERSE 

PRINT " " 
PRINT " PARTITIONEN " 
PRINT " " 
NORMAL : PRINT 

INPUT "WELCHE ZAHL ? "5N 
DIM P(N) 


REM ANFANGSWERTE 
Pi =N 
K=-1:F=-9272=9 

PRINT 


PRINT "PARTITIONEN VON "5N 
PRINT 

REM UNTERPROGRAMMAUFRUF 
GOSUB 299 
z=2I+]1 

FOR J=1T0K 

PRINT PiJ)5 

NEXT J: PRINT 

IFF = 1 THEN 198 

PRINT 

PRINT Z5" PARTITIONEN" 
END 

REM UNTERPROGRAMM 

IFF = 1 THEN 399 
F=ei 

FOR I=1TOK 

IF P(I) = 1 THEN 519 
NEXT I 


I=K 
GOTO 319 


L=K-]1I 
K=1I 
P{ID = P(I) - i 
IF PiK) > L THEN 488 
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49 L=L- 

AOK=SKHri 
458 PikK) = Pik - 1) 

4608 GOTO 428 

A70 : 

485 K=K +1 

498 PiK) = L + i 

S28 IF PIIDD) < >1LTHENI=K 
18 IF PD) =1ITHENI=1- 
522 IFI=B8THENF=9 

332 RETURN 
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PARTITIONEN 


WELCHE ZAHL ? 18 


PARTITIONEN VON 19 


322111 
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Ss211111 
Ssi111111 
22222 
222211 
2221111 
22111111 
zilliillil 
1111111111 


42 PARTITIONEN 


97 


Hırronrnı CarDanı 
relinquiturprima 6 m:rz 30$,hx autem quantitates proportionales 
Fine, & quadrauum feaında et aquale duplo produdti fecund& in 
primam,cum quadruplo primx,ut proponebatur. 


De cubo&rebusaqualibusnumero. Cap. XI, 


Cipio Ferreus Bononicnfis iam annis ab hincrriginta fer 
mie capısulaın hoc inuenit , tradidit uero Anthonio Mas 
Ey ra Plorıdo Vcneto,qui cu incertamen cu Nicolao Tars 
en ralca Brixclienfe aliguando ucnifler, occafionemdedit, ue 
Nicolaus inuencric & ipfe,qui cum nobis rogantıhus tradıdıffer, fup 
prclla demonitratione, frei hoc auxilio, demonftrarionem quclinie 
mus,eam&g in modos, quod diflicıllımum fuic, redadtam fic fubiccis 
mus. Desonstrartıo, 

Sit igieur cxemplı caufa cubus G H & [cxcuplum lateris « u zqua 
le 20,&% ponam duos cubos a E& c 1,orum dıflerentia lit 20 ‚ita 
cuodproductum A c larerıs, inc x lacuıs, 
fir 2, tertia fcilıcet nunierircrum pars, & 
abfeindam c 5, equalem c x,dıco, quodfi F 3 
ira fucrir,lineam A 5 rcfiduum ‚cfle zqua» 
Icm 6 #,& idco rei altımationem, nam de 
6 Hıam fupponebatur,quod ita cflet, per G 
havam igitur per modum primı fuppofiti K 

.s'capiiulı huius libri, corpora DA,DC,DE A 
‚D p,ur per Dc in.clligamus cubum B c,per 
p Fr cubum A B,per D A triplum c Bın quadrarum A B,per De triplum 
ausm quadratü B c.quiaigitur x AcineKÄfir2,cxAcimcK rer Act 
« nuımzrus rerum, igitur cg AB triplum Acinc K fiunt6 resaB, 
fcu fexcuplum A B,quarc triplum produdti ex AB, B C,A c,eltfexcus 
plum a B,at uero difterenria cubı A c,axubo CK ‚& exiltenti ä aıbo 
3 ccı @äleex fuppofiro,cht 20,& ex fuppofio prımo 6 capirulı ‚el 
aygregatum corporum D A,D E,D F,triaigitur hırce corpora func 20, 
pofita tıero o c m:cubus A B,qualıs eft cubo a c,& trıplo A cinqua 
dıamum c 8,& cubo B c m:& rriplo 8 c in quadratum a c m: per des 
monftrata illıc,diflerencia aurem tripl 8 c inquadracum A c,ä triplo 
ac in quiadracrum a ceft productum A B,B CA C,quarc cum hoc,ut de 
monltratum eft,icquale fir fexcuplo A 5, igitur addıro fexcuplo AB, 
adıd quod Arex A c inquadratum a c ter, Act triplum 8 c ın quadras 
tum A C,cum iQitur a c hir m:iam oftenlum eft,quod produstum ca 
N) 


"Erste Seite aus der Ars Magna von Geronimo Cardano, 
Nürnberg 1545” 
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17. KUBISCHE GLEICHUNG 


Die kubische Gleichung 
axd+bx2+cx+d=0 


kann mit Hilfe der Cardano-Formel gelöst werden. Die Cardano-Formel ist 
für das Rechnen von Hand äußerst umständlich, so daß man meist zu einem 
Iterationsverfahren greift (vgl. Programm 31). Sie liefert jedoch im Gegensatz 
zu den Iterationsverfahren sämtliche Lösungen der Gleichung. 


Dividiert man die kubische Gleichung durch a und führt die Substitution 


b 
= + —. 
yax 3a 


durch, so erhält man die sogenannte reduzierte Gleichung 
y+t3py+2q=0 
mit 
2b3 be, d 3ac — b2 
2q = —__ _ _—_ + _ 3p = 
2723 32 a 3a2 
Aus p und q berechnet sich die kubische Resultante 
D= q2 + p3 


Ihr Vorzeichen gibt Auskunft über die Art der Lösungen: 


Ist D> 6, so gibt es 1 reelle, 2 konjugiert komplexe Lösungen 
Ist D < 0, so gibt es 3 reelle Lösungen 
Ist D = @, so gibt es mehrfache, reelle Wurzeln 


Mit Hilfe der 3. Wurzel 


Y -qatsVD 


können nun die gesuchten Lösungen berechnet werden. Im Fall D<@ aller- 
dings, ist in der Cardano-Formel eine komplexe Wurzel zu ziehen. Deswegen 
wird hier eine trigonometrische Lösung durchgeführt. Die hierbei benötigte 
Umkehrfunktion arccos(x) der Cosinus-Funktion ist nicht in BASIC imple- 
mentiert. Sie muß daher über die Arctangens-Funktion ausgedrückt werden. 
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Es gilt 


vi-x2 fürO<x<1 
X 
Vi,n für_I1sx<öO 


arctan 
arccosx 7 arctan 


a fürx = 0 
2 


Zum folgenden Programm 


Das Programm berechnet die Lösung der kubischen Gleichung nach dem an- 
gegebenen Verfahren von Cardano. 


Eine berühmte kubische Gleichung ist 
x3 + 2x2 + 10x - 20=0 


Sie wurde Leonardo von Pisa (ca. 1170 — 1250), genannt Fibonacci, von 
Johann von Palermo, einem Höfling von Friedrich Il, zur Lösung gestellt. 
Fibonacci fand die Lösung 


x = 1.3688081075 
die auf 8(l) Dezimalen genau ist. 


Die Eingabe der Koeffizienten ins Programm liefert, neben der genannten Lö- 
sung, noch die beiden komplexen Lösungen 


—1.68440405 + 3.43133135i 
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REM KUBISCHE GLEICHUNG 


PI= 3.141592654 


HOME : INVERSE 


PRINT " " 
PRINT " KUBISCHE GLEICHUNG " 
PRINT " " 
PRINT " AX"3+BX24CX+D=B " 


PRINT " n 

NORMAL : PRINT 

PRINT "EINGABE DER KOEFFIZIENTEN" 

PRINT 

INPUT "A, B, GC, D ? "5A,B,C,D 

PRINT 

IF A= @ THEN PRINT "KEINE KUBISCHE GLEICHUNG": END 
A=AA=B/AU:B=C/AULC=D/ A 


REM PARAMETER DER REDUZ.GLEICHUNG 
= -AxA/9I+B/S 
Arash /27-ArB/& +C/2 


P 
Q 
REM KUBISCHE RESULTANTE 

R=QO+Qa+PxPxP 

REM FALLUNTERSCHEIDUNG 

ON SON (R) + 2 GOTO 490,419, 298 

REM i REELLE,2 KOMPLEXE LOESUNGEN 

U= 850R (R) - U = SGN (VD) * ABS (WU) * (1/9 
v= - SOR (RM - Q:V= SGN (V #* ABS (V) * (IL /D 
= U+V-A7/IS 

REM  REAL- UND IMAGINAERTEIL 


2= -(U+rW7/2-A!/S.ıX5=WUW-V) / 2% SOR (3) 
PRINT "Xi1="5X1 

PRINT 

PRINT "X2="5X25 "+1#"5X5 

PRINT 


PRINT "X3="5X2; "-I#"5X3: END 


REM REELLE DOPPELLOESUNG 
U=(-9)*(1/D 
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Yı =2U-A/3 
{2= -U-A/S 

PRINT "X1="5X1 

PRINT 

PRINT "X2=X3="5X2: END 


REM 3 REELLE LOESUNGEN 
= -07 SORt-P*MI 

REM DEFINITION DES ARCCOS 

DEF FN AC(X) = ATN ( SOR (1 -XX  / MD 
IF X > 8 THEN PHI = FN AC(X): BOTO 568 

IF X < 8 THEN PHI = FN AC(X) + PI: GOTD 568 
PHI=PI/2 

Xi =2%* SOR(-P) * COS (PH /IM -A/3 


2= -2%* S50R(-P) * COS ((PH+PD /ID) -A/S 
S= -2# SOR(t-P) * COS ((PH-PND /I -A7/3 


PRINT "X1=e"5X1 
PRINT 

PRINT "X2="jX2 
PRINT 

PRINT "X3="5X3 
PRINT 

END 
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KUBISCHE GLEICHUNG 


AX"S+BXT2+CX+D=9 


EINGABE DER KOEFFIZIENTEN 
A, B, C, D ? 1,2,19,-28 
X1=1.36880811 
X2=-1.68449405+1#3.43133135 


X3=-1.68440405-1%#3.43133135 
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18. GLEICHUNG VIERTEN GRADES 


Ebenso wie die kubische Gleichung kann auch die Gleichung 4. Grades 
axd +bx3+cx2+dx+te=d 


mit Hilfe von Wurzeltermen gelöst werden. Nach Ferrari kann die Lösung der 
Gleichung 4. Grades auf die kubische Gleichung 


8y3 — Acy2 + (2bd-Be)y + el4c-b2) - d2=0 
zurückgeführt werden. Setzt man eine Lösung y in die quadratische Gleichung 
1 by-d 
+ — (b+A)x + (y+t ——) = 
xt, (b+A)x + (y r\ )=0 
mit 


A= t,V8y+b2 -4c 


ein, so liefert deren Lösung schließlich die gesuchten Lösungen der ursprüng- 
lichen Gleichung 4. Grades. 

Es ist klar, daß die Lösung nach Ferrari für das Rechnen von Hand zu kompli- 
ziert ist. Daher werden Gleichungen 4. Grades meist mit Iterationsverfahren 
angegegangen. Da jedoch Gleichungen 4. Grades nicht notwendig reelle Lösun- 
gen haben, kann es hier zu Konvergenzschwierigkeiten kommen. Da aber die 
Gleichungen vom 3. und 4. Grad häufig vorkommen — insbesondere im Zu- 
sammenhang mit dem charakteristischen Polynom einer Matrix — ist das vor- 
liegende Programm sehr nützlich. 


Nachdem Ferrari 1545 seine Formel veröffentlicht hatte, bemühten sich die 
Mathematiker — insbesondere Euler und Lagrange — über 250 Jahre lang eine 
entsprechende Formel für die Gleichung 5. Grades zu finden. Erst 1824 zeigte 
Abel, daß Gleichungen von höherem als dem 4. Grad prinzipiell nicht mit 
Radikalen, d.h. mit Wurzeltermen zu lösen sind. Hier muß man stets zu Itera- 
tionsverfahren greifen. 


Zum folgenden Programm 


Das Programm löst die allgemeine Gleichung 4. Grades nach dem Verfahren 
von Ferrari. Die dabei auftretende kubische Gleichung wird nach der Cardano- 
Formel von Programm 17 gelöst. 


Als Beispiel wird die Gleichung 
x4 — 24x3 + 150x2 — 200x - 375 = 0 


behandelt; sie ergibt sich in’Programm 36 als charakteristisches Polynom der 
dort angegebenen Matrix. 


Das Programm liefert die reellen Lösungen 
—], +5, B, 15 
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REM "GLEICHUNG VIERTEN GRADES" 


PI = 3.141592654 


4. GRADES”: END 


HOME : INVERSE 
PRINT " 
PRINT " GLEICHUNG VIERTEN GRADES 
PRINT " 
PRINT " AxX%4 + BRXYI + CaX*2 + DaX + E " 
PRINT " 

NORMAL : PRINT 

INPUT "A, B, C, D, E ? "}A,B,C,D,E 
PRINT 

PRINT "LOESUNGEN: ": PRINT 

IF A= 8 THEN PRINT "KEINE GLEICHUNG 

IF B= 8 AND D = 8 THEN 829 

REM KUBISCHE GLEICHUNG 


GDSUB 339 

REM GUADRATISCHE GLEICHUNG 
J=8*Y+BxsB-4+*«C 

IF ABS (J) <SE-8THENJ = 9 
Hl = S0R (J):G = Hi 
H (BE #* Hi + )) / 2 

I Y+H+Bx*VY-D 

GOSUB 558 


HIi= -HIEH= (B*Hi + J) / 2:6 = H 
I=sYsHi+B*Y-D 

GOSUB 559 

END 


REM KUBISCHE GLEICHUNG 
All= - C/2:B1=B#D/4-Eili = 
8 

= -Al#*A1 /9+B1 7/35 
=Al*37/27-Al* Bl /6+C1/2 
REM KUBISCHE RESULTANTE 
R=Q0*Q+P*3 


pP 
Q 


ON SGN ({R) + 2 GOTO 3359, 480,429 


E*(4A*C-B*B) -D*D) / 


428 REM 1 REELLE, 2 KOMPLEXE LOESUNGEN 


48 U= SUR (R - HU= SGN (VD) * ABS (UV) * (1 /D 
49 V= - SOR (RD - B:EV= SGN (VW #* ABS (V) * (1 /D 
450 YwU+r+V-pAL/3 

45% RETURN 

470 : 


488 REM REELLE DOPPELLOESUNG 
WB U=-D) *L/MD 

SB Y=2#U- AL 7/3 

318 RETURN 


352 REM 3 REELLE LOESUNGEN 

SBBX= -090/ SORI-P* MD 

ssß REM DEFINITION DES ARCCOS 

9 DEF FNACKK) = ATN (SOR (L-X XD) / X 
780 IFX >98 THEN PHI = FN ACtX): GOTO 590 

s88 IFX < 68 THEN PHI = FN AC(K) + PI: GOTO 6%9 
398 FHI = PI/2 

BE Y=2* SOR (-P) * COS {PHI /3 - Al /3 
612 RETURN 

628 : 

638 REM GQUADRATISCHE GLEICHUNG 

6409 DII=SH*H- 46 #1 

658 IF DI < 8 THEN 748 

b6B : 

678 REM REELLE LOESUNG 

688 Xi = (-H + SOR (DI) / (2 * 6 

698 X2 = -H- SOR (DI) / (260 

788 PRINT X1: PRINT : PRINT X2 

718 PRINT 

728 RETURN 

730 5 

748 REM KOMPLEXE LOESUNG 

758DI= -DI | 

7668 PRINT -H 7/7 (2 *G);5" + I#*"5 SOR (DD) / (2 * © 
778 PRINT 

788 PRINT -H 7 (2 *G)5" + I"; - SQR (DI) / (2 * © 
778 PRINT 

898 RETURN 

gißı: " 

828 REM BIQUADRATISCHE GLEICHUNG 

SSOK=-C HC -AHrAHE 

849 IF K< 2 THEN 1999 

859 YyYı =(-C+ SIR IC CC -AHAxE)) / (2A 
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868 Y2 = ( 


879 
21-77) 
879 
708 
919 
929 
750 
748 
9582 
768 
970 
988 
970 
1098 
1818 
1028 
1030 
1040 
1950 
10590 
1278 
1082 
1999 
1108 
1119 
1128 
1138 


IF Yı 
PRINT 
PRINT 
GOTO 9 
PRINT 
PRINT 
IF Y2 
PRINT 
PRINT 


U 
ee 
[ug | Bauen | u | u || 


R2 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
END 


-C- SORIE*HE -ArAHE)) / (2A 
< 8 THEN 918 

SOR (Yidz PRINT 

- SAR (Yı)z PRINT 

39 

"I1x"5 SOR t - Yi)z PRINT 

"I#"3 - SOR t - YDDE PRINT 
< 8 THEN 969 

SQaR (Y2)s PRINT 

- SQR (Y2): PRINT : GOTO 988 
"I#"s3 SOR ( - Y2): PRINT 

"I#"5 - SOR t - Y2)z PRINT 


-C/ (2 *A):IM= SR (-K / 2 HM 

SER {RE * RE + IM « IM) 

SaR ({BT + RE) / 2:11 = SQOR ((BT - RE) / 2 
SOR (Ri # Ri + I 1 * II) 

SOR (BE + Ri) / 23:12 = S5SOR (BE - Ri) / 2 
R2; "#1"5 - 12 


R2; "#1"312 
I R2: "I": = 12 


- R23 "1#"512 
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GLEICHUNG VIERTEN GRADES 


ARX“A + BRrX"3 + C#X”2 + DAX + E 


A, B, C, D, E ? 1,-24, 158, -208, -375 
LOESUNGEN: 

4.999994 

-,999999993 

5.00899599 


15 
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19. POLYNOMBERECHNUNG 


Das folgende Programm bestimmt zu vorgegebenen Nullstellen das zugehörige 
Polynom. 


Die Koeffizienten a; (1 Si sn) eines normierten Polynoms n-ten Grades 
p(x) =xNn + an _ 1x0 I +2 xn-2+an_1x0-1+...+00 


sind nach dem Satz von Vieta elementarsymmetrische Funktionen der Null- 
stellen x;: 


xıtx2+x3+....txXn”= -an—-]1 
x1x2 + x1x3 +x1x4t...tXn-1Xn” an-?2 


x1X2X3 + X1X2x4 + X1X2X5 +... + Xn-2 Xn—1 Xn= —an-3 


xX1'X2°X3-X4°...."Xn= (10.00 


Zum folgenden Programm 


Das Programm berechnet die Koeffizienten des gesuchten Polynoms nach 
oben genannten Formeln iterativ innerhalb zweier verschachtelter Schleifen 

220 A(d)=9:A(N)=1 

238 FOR I=N TO 1STEP —1 

248 FOR J=H TO 1STEP —1 

258 AIN-J) = AIN-J) — AIN-J+1) * X{N—I+1) 

260 NEXT J 

2760 NEXTI 


Als Programmbeispiel wird das Polynom mit den Nullstellen 
xı=1, x2=2, x3=3, x4=4, x5=5undxg=6 
gesucht. Es ergibt sich das Polynom 
P(x) =x6 — 21x5 + 175x4 — 735x3 + 1624x2 — 1764x + 720. 


Das Programm funktioniert nur für reelle Nullstellen. Für komplexe Nullstellen 
müssen die entsprechenden quadratischen Polynome multipliziert werden, 
dies kann mittels Programm 21 erfolgen. 
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REM POLYNOMBERECHNUNG 


HOME : INVERSE 
PRINT " n 

PRINT " POLYNOMBERECHNUNG " 

PRINT " 
NORMAL : PRINT 

INPUT "WIEVIELE NULLSTELLEN ? "}N 
DIM AN) ,X(N) 


PRINT 
FOR I=1 TON 
PRINT I15".NULLSTELLE "5: INPUT X(D 


AI) =9 
NEXT I 
PRINT " " 
NORMAL : PRINT 
AB) = DEAN = 1 
FORI=SNTOISTEP - 1 
FÜR J=NTOASTEP -1 
AN -W =EAN- I) -AN-J+1) «XN-I+N) 
NEXT J 
NEXT 1 


PRINT "DIE KOEFFIZIENTEN SIND: " 
PRINT 

FORI=SNTOBSTEP -1 

PRINT "AUS = "5 

IFA > =®THEN PRINT ® "; 
PRINT A(T) 

NEXT I 

END 
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POLYNOMBERECHNUNG 


WIEVIELE NULLSTELLEN ? 6 


1.NULLSTELLE 
2.NULLSTELLE 


S.NULLSTELLE 


4. NULLSTELLE 


3. NULLSTELLE 


6. NULLSTELLE 


DIE KOEFFIZIENTEN SIND: 


A(6) 
AS) 
A(4) 
A(3) 
A(2) 
AD) 
A(9) 


1 
-21 
175 
-7535 
1524 
-1764 
720 


21 
?2 
73 
?4 
?5 
?6 
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20. KOMPLEXES HORNERSCHEMA 


Mit Hilfe des Hornerschemas kann der Wert des Polynoms 
Fix) = an tan _ıxn-1 Han _a2xNn-2+...+ 00. 


an der Stelle x ohne Potenzieren ausgewertet werden. Dies kann mit folgen- 
dem Programmstück geschehen 


100 F=A(N) 

116 FOR I=N-1TO BSTEP —1 
126 F=F*X+All) 

130 NEXTI 


Dabei ist N der Polynomgrad und F der gesuchte Funktionswert. Das 
Hornerschema ist auch im Betriebssystem von CBM-Rechnern enthalten; mit 
ihm werden alle transzendenten Funktionen wie sin, arctan usw. über eine 
Polynomapproximation berechnet. Das angegebene Hornerschema kann auch 
für komplexe Zahlen angewendet werden. Dazu muß die komplexe Addition 
und Multiplikation geeignet sein. Sind die komplexen Zahlen als Zahlenpaare 


(RE(1), IM(1)) (RE(2), IM(2)) 
gegeben, so kann die Addition wie folgt definiert werden 
500 REM KOMPLEXE ADDITION 


51d RE = RE(1) + RE(2) 
520 IM = IM(1) + IM(2) 


Entsprechend die Multiplikation 


608 REM KOMPLEXE MULTIPLIKATION 
618 RE = RE(1) * RE(2) — IM(1) * IM(2) 
620 IM = RE(1) * IM(2) + RE(2) * IM(1) 


Diese Rechenoperationen werden als Unterprogramm definiert und bei Bedarf 
angesprungen. Da jedoch in BASIC alle Variablen global sind, muß zuvor eine 
geeignete Parameterübergabe durchgeführt werden. 


Da hier die komplexen Zahlen als Zahlenpaare behandelt werden, müssen ent- 
sprechend auch die Koeffizienten des Polynoms und der Argumentwert als 
Zahlenpaare eingelesen werden. 
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Zum folgenden Programm 
Als Beispiel wird der Funktionswert des komplexen Polynoms 
ptz) = (1 + 1.5129 — 22 +32 + (1 + 2i) 
an der Stelle 
z=-1-i 
berechnet. Die Koeffizienten werden im Programm mittels DATA-Werten ein- 


gelesen. 
Es ergibt sich der Polynomwert 


5-4i. 
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REM KOMPLEXES HORNERSCHEMA 


HOME : INVERSE 

PRINT " n 
PRINT " KOMPLEXES HORNERSCHEMA " 
PRINT " " 
NORMAL : PRINT : PRINT 

READ N: REM POLYNOMGRAD 

DIM R{N),ItN),RE(2),IM(2) 


REM EINLESEN DER KOEFFIZIENTEN 
FORK= N TO 8STEP - 1 
READ RK) ,I(K) 

NEXT K 

READ X,Y: REM ARGUMENT 

REM HORNERSCHEMA 

RE(1) = R{NSIM(I) = IN 
FORK=N-1 TO 8STEP - |] 
RE(2) = X: 1IM(2) = Y 

REM MULTIPLIKATION 

GOSUB 459 
RE(1) = RE:IMil) = IM 


RE(2) = R{K)zIM(2) = IR) 
REM ADDITION 


GOSUB 38% 
RE(1) = RE:IM(I) = IM 
NEXT K 


REM AUSGABE 
PRINT "POLYNOMWERT = "5 


IF IM) > = 9 THEN PRINT RE(1)5" + I#"51M(1): GOTO 368 


PRINT RE(1)3" -I#"; ABS (IM(iD) 
END 


REM UNTERPROGRAMM KOMLEXADDITION 
RE = RE(1) + RE(2) 
IM = IMt1) + IM\2) 

RETURN 
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438 REM UNTERPROGRAMM KOMPLEXMULTIPLIKATION 
448 RE = RE(1) * RE(2) - IM(1) * IM(2) 

458 IM = RE(1) * IM(2) + RE(2) * IM) 

468 RETURN 

470 : 

488 DATA 3 

498 DATA 1,1.5,-1,9,3,0,1,2 

598 DATA 1,-1 


KOMPLEXES HORNERSCHEMA 


POLYNOMWERT = 5 -I#4 
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21. POLYNOM-MULTIPLIKATION 


Zwei Polynome vom Grad n und m 
an +tan_ 1x1 +an_2xn-2+...+20 
bmxM + bm—_1xm-1+bm_2xm-2+...+bo 


werden multipliziert, indem man jeden Term des ersten Polynoms mit jedem 
des zweiten multipliziert. 


Dies kann mit Hilfe zweier verschachtelter Schleifen geschehen 


2968 FOR I=8 TON 

300 FOR J=@9TOM 

318 K=1+J 

328 C(K)=CIK) + Al) * BU) 
3380 NEXT J 

349 NEXT I 


Das Verfahren läßt sich direkt auch auf Reihen übertragen ; d.h. es liefert auch 
das Cauchy-Produkt zweier Reihenentwicklungen. 


Zum folgenden Programm 


Das Programm berechnet nach dem angegebenen Verfahren das Produkt 
zweier Polynome und gibt die entsprechenden Koeffizienten aus. 


Als Beispiel werden die Polynome 
p(x) =x4 _5x3+3x2+6x— 2 -%0 
ax)=x2-1 


multipliziert. Eingabe der Koeffizienten von links nach rechts ins Programm 
liefert die Koeffizienten des Produktpolynoms 


1,—-5,2, 11, -5, -6,2 


Das Ergebnis ist somit 
x6 _ 5x5 + 2x4 + 11x93 — 5x2 - 6x +2 
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Ijerlil 


"Darstellung des Polynoms 
2x3 + 15x2 + 166x — 4460 = 0 


in Neue Schritte der Berechnung von Li Ye, 
China 1259 
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REM POLYNOM-MULTIPLIKATION 
HOME : INVERSE 
PRINT " " 
PRINT " POLYNOM-MULTIPLIKATION " 
PRINT " " 
NORMAL : PRINT 

INPUT "GRAD DES 1.POLINOMS ? ";N 
PRINT 

DIM AN) 

FOR I=N TO BSTEP - 1 

PRINT "KOEFFIZIENT VON X("5 15)”; 
INPUT A(D) 

NEXT I 

PRINT 

INPUT "GRAD DES 2.POLYNOMS ? "5M 
PRINT 

DIM B(M),C(N + WM 

FOR I=MTOBSTERr - 1 

PRINT "KOEFFIZIENT VON X("515'")";5 
INPUT B(D) 

NEXT I: PRINT 


K=1I+) 

C{K) = ECiK) +AID) * Bi) 
NEXT J 

NEXT I 

REM AUSGABE 

PRINT "PRODUKTPOLYNDM: " 


PRINT 

FOR I=N+MTOBSSTER -1 
PRINT C(DD5" "I 

NEXT I 

END 
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POLYNOM-MULTIPLIKATION 


GRAD DES 1.PODLINOMS ? 4 


KOEFFIZIENT VON X(4)?1 
KOEFFIZIENT VON X(3)?-5 
KOEFFIZIENT VON X(2)73 
KOEFFIZIENT VON X(1)76 
KOEFFIZIENT VON X(9)7-2 


GRAD DES 2.POLYNOMS ? 2 
KOEFFIZIENT VON X(2)?1 
KOEFFIZIENT VON X(1)78 
KOEFFIZIENT VON X) ?-1 
PRODUKTPOLYNDM: 


1-5 211-5 -6 2 
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22. POLYNOM-DIVISION 


Die Polynom-Division findet vielfache Anwendung beim Abdividieren von ge- 
fundenen Nullstellen, bei der Partialbruchzerlegung u.s.w. 


Die Polynom-Division kann in völliger Analogie zur gewöhnlichen Division 
durchgeführt werden: 

Es wird geprüft, wie oft der entsprechende Koeffizient des Divisors in den 
Dividenden hineingeht. Das entsprechende Vielfache des Divisors wird dann 
vom Koeffizienten des Dividenden subtrahiert. 


Zum folgenden Programm 
Das Programm dividiert zwei Polynome nach dem gewöhnlichen Divisionsver- 
Es wurde versucht Rundungsfehler aufzufangen, indem sehr kleine Reste Null 
gesetzt werden. Die Polynomgrade und Koeffizienten werden in Form von 
DATA-Werten eingelesen. 
Als Beispiel werden die Polynome 

2x7 + 9x6 — 5x5 + 42x4 — 24x3 + 10x2 + 21x — 8 

x3+5x2 - 3x +4 
dividiert. Das Programm liefert den Quotienten 

2x4 _x3+6x2+x—7 
mit dem Rest 

24x2 — 4x + 20 
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REM POLYNOMDIVISIDN 


HOME : INVERSE 


PRINT " " 
PRINT " POLYNOMDIVISION " 
PRINT " " 
NORMAL : PRINT 

READ N 


DIM A{N + 1),B{N + 1),Q0(N + 1),R(N + 1) 


FOR I=N+1TO1STEP - 1 
READ A(I): REM DIVIDEND 
NEXT I 

READ M 

FOR I=M+1 TO 1STEF - 1 
READ B(I): REM DIVISOR 
NEXT I 


PRINT "DIVIDEND: ": 

FOR I=N+1 TO 1STEP - 1 
PRINT A(D5" "5 

NEXT I: PRINT 

PRINT 

PRINT "DIVISOR: ": 
FORT=M+1TO1STEP -1 
PRINT B{DD;" "5; 

NEXT I: PRINT 


R=N-M 

FORI=N+ITOM+1STEP -i 

Qg(I-M =A(I /BiM +1) 

FOR J=8TOM 

AT-J) =AUL- ID -QUI-M *BIM+I- UT 


IF ABS (AL -J)> = ABS (BII-M *BiM+I1 - 
J) * IE - 8 THEN 370 
AI-J)=9 

NEXT J 
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428 NEXT I 

439 : 

448 PRINT 

459 PRINT "QUOTIENT: "3 

48 FORI=SR+ILTOLSTE - 1 
478 PRINT Q(D)5" "5 

488 NEXT I: PRINT 

498 : 

598 PRINT 

518 PRINT "REST: IF 

528 IF R(R) = Ö THENR = R - 1: GDTO 528 
S38 FORI=RTOISTE -1 

548 PRINT R(DI5" "; 

SSQ NEXT I 

568 PRINT 

578 END 

589: 

590 DATA 7 

600 DATA 2,9,-5 „42, -24,18,21,-8 
10 DATA 3 

628 DATA 1,5,-3,4 


POLYNOMDIVISION 


DIVIDEND: 29-5 42 -24 18 21 -B 
DIVISOR: 15-54 
QUOTIENT: 2-1651-7 


REST: 24 -4 28 
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23. MATRIZEN-MULTIPLIKATION 


Matrizen sind rechteckige Anordnungen von reellen Zahlen. Eine Matrix hat 
die Ordnung (m, n), wenn sie m Zeilen und n Spalten besitzt. 


Zwei Matrizen A und B werden multipliziert, indem man jede Zeile von A mit 
jeder Spalte von B multipliziert. Dies zeigt, daß das Produkt von A und B nur 
existiert, wenn die Spaltenzahl von A mit der Zeilenzahl von Bübereinstimmt. 
Ist A von der Ordnung {n, m) und B der Ordnung (p, q) mit m=p, so ist die 
Produktmatrix von der Ordnung {n, q). 


Die Matrizenmultiplikation kann mit Hilfe von folgendem Programmstück 
durchgeführt werden 


320 FOR I=1 TON 

330 FORJ=-1TOQ 

340 S=0 

350 FORK=1 TOM 

366 S=S + All,K) * B(K,J) 
376 NEXTK 

380 C(1,J)=S 

398 NEXT J 

480 NEXT I 


Zum folgenden Programm 


Das folgende Programm berechnet das Produkt zweier Matrizen, die in Form 
von DATA-Werten eingelesen werden. Dabei wird geprüft, ob die Spaltenzahl 
der ersten Matrix mit der Zeilenzahl der zweiten übereinstimmt, andernfalls 
wird eine Fehlermeldung ausgegeben. 


Im Programmbeispiel werden die Matrizen 


8 5 2 3 1 — 

6 4 1 1 —2 3 

1 2 2 4 1 1 

1 2 

miteinander multipliziert. Die Produktmatrix ist 

3 — 
—1 1 
3 14 


(vgl. Programmausdruck). 
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REM MATRIZENMULTIPLIKATION 


HOME : INVERSE 
PRINT " " 
PRINT " MATRIZEN-MULTIPLIKATION " 
PRINT " " 
NORMAL : PRINT 
READ N,M: REM ORDNUNG DER 1.MATRIX 
DIM AtN,M) 

FOR I = 1 TON 

FOR J=1 TOM 

READ A(1,J) 

NEXT J 

NEXT I 

READ P,@: REM ORDNUNG DER 2.MATRIX 

DIM B(P,O) 

IF M = P THEN 228 


PRINT "SPALTENZAHL DER 1.MATRIX MUSS MIT DER " 
PRINT "ZEILENZAHL DER 2.MATRIX UEBEREINSTIMMEN": END 
FOR I=1T0P 

FORJ=17T00 

READ B(1,0) 

NEXT J 

NEXT I 

FORI=1 TON 

FORJ=1700 

s=9 

FORK=1 TOM 
5= 5 + A(I,K) * B(K,J) 

NEXT K 
CL) =8S 

NEXT J 

NEXT I 

PRINT "PRODUKTMATRIX: " 

PRINT 

FOR I = 1 TON 

FORJ=1T0Q 
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420 
430 
449 
450 
469 
470 
480 
490 
302 
310 
328 
330 
549 
JR 
1-17) 
378 
89 


IF C(1,D > 


= @ THEN PRINT " "5 


PRINT C(I,J)3" "3 


NEXT J 
PRINT 
NEXT I 


END 


DATA 
DATA 
DATA 
DATA 


DATA 
DATA 
DATA 
DATA 
DATA 


MATRIZEN-MULTIPLIKATION 


PRODUKTMATRIX: 
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24. ZWEIERPOTENZEN VON MATRIZEN 


Da für quadratische Matrizen die Zeilenzahl mit der Spaltenzahl überein- 
stimmt, können quadratische Matrizen mit sich selbst multipliziert werden. 


Solche Potenzen von Matrizen haben zahlreiche Anwendungen in der numeri- 
schen und Wirtschaftsmathematik, Wahrscheinlichkeitsrechnung und Graphen- 
theorie. 


Als Beispiel werden hier die (homogenen) Markowketten herausgegriffen: 

Hat ein System endlich viele Zustände und geht ein Zustand i, unabhängig 
von seinem vorhergehenden Zustand, in den Zustand j mit der Wahrschein- 
lichkeit 


Pij 
über, so stellt das System eine Markowkette (auch Markoff geschrieben) dar. 
Die Matrix P mit den Elementen 


Pij 
heißt eine stochastische Matrix. Solche Matrizen sind dadurch gekennzeichnet, 
daß ihre Zeilensumme den Wert 1 haben, da die Summe aller Übergangswahr- 
scheinlichkeiten natürlich den Wert 1 hat. 
Hat dieses System ein stationäres Verhalten, d.h. strebt es einer Grenzvertei- 
lung zu, so kann diese Grenzverteilung durch fortgesetztes Potentieren berech- 
net werden [10]. 
Multipliziert man die Matrix mit sich selbst und quadriert das jeweilige Ergeb- 
nis, so erhält man die Zweierpotenzen der Matrix. 


Zum folgenden Programm 


In einem bestimmten Land gebe es 3 Parteien A, Bund C. Folgendes Wahlver- 
halten sei gegeben: 

80% aller A-Wähler wählen bei der nächsten Wahl wieder A, 10% B und 
10% C. 

70% aller B-Wähler wählen bei der nächsten Wahl wieder B, 15% A und 
15% C. 

60% aller C-Wähler wählen bei der nächsten Wahl wieder C, 20% A und 
20% B. 
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Damit ergibt sich folgende stochastische Matrix 


0.8 9.1 9.1 
0.15 0.7 0.15 
0.2 0.2 0.6 


Wie man dem Programmausdruck entnimmt, nähern sich die Zeilen der Matri- 
zenpotenzen dem Vektor 


(0.461912194 0.307419108 0.230668699) 


Dieser stellt die Grenzverteilung dar; d.h. auf lange Sicht wählen 


46.2 % aller Wähler A 
30.7% B 
231% C 
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REM ZWEIERPOTENZ EINER MATRIX 


20 HOME : INVERSE 

39 PRINT " " 
48 PRINT " ZWEIERPOTENZ EINER MATRIX " 
58 PRINT " n 
60 NORMAL : PRINT 

78: 

80 INPUT "WELCHE ZWEIERPOTENZ ? "}P 
96 P= LOG (P) / LOB (2) 

188 : 

11dö PRINT 

120 INPUT "ORDNUNG DER MATRIX ? "IN 
1308 PRINT 

140 PRINT "GIB MATRIX ELEMENTWEISE EIN..." 
150 PRINT 

168 FOR I=1 TON 

170 FOR J=1 TON 

188 INPUT A(I,.) 

198 NEXT J 

208 NEXT I 

218 : 

228 FORL = 1 TOP 

238 PRINT 

248 PRINT 2 * L$".POTENZ:" 

259 FOR I=1 TON 

260 FOR J = 1 TON 

2795=9 

268 FORK = 1 TON 

290 8 = 58 + A(I,K) * AlK,J) 

388 NEXT K 

318 B(I,J) = 8 

328 IFS> = THEN PRINT " "3 

339 PRINT 55" "} 

340 NEXT J: PRINT 

359 NEXT I 

360 : 

376 REM UMSPEICHERN 

388 FOR I =1 TON 

398 FOR J = 1 TON 

488 Al1,J) = B(I,)) 

410 NEXT J 

420 NEXT I 

438 NEXTL 

449 END 
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ZWEIERPOTENZ EINER MATRIX 


WELCHE ZWEIERPOTENZ ? 16 


ORDNUNG DER MATRIX 7? 


3 


GIB MATRIX ELEMENTWEISE EIN... 


en 


u 


IV II ID NV 
rNDNDeNe mr (0 


2. POTENZ: 
.675 .17 
‚255 
.Si .28 

4.POTENZ: 
‚247025 
.37385 
40775 


8. POTENZ: 
„475442597 
448546516 
453719452 


16. POTENZ: 
„461912194 
„461128662 
461337397 


935 


„2491 
„388375 
„3173 


„155 


„2i 


„Ai 


„27495 


2976977678 
„319420823 
„312043548 


„327419108 
. 307973985 
«387836579 


„203875 
257975 


2268597726 
„234032661 
234237091 


250668599 
230877434 
. 23588260825 
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25. MITTELWERTE 


Neben dem arithmetischen Mittel A der Zahlen x; (1 Si<.n) 
=. +x2+x3+X4...+Xn) 


sind auch andere Mittelwerte von Bedeutung: 
das geometrische Mittel 


G=Y xX1'X2°’X3'X4°....*Xn 


das harmonische Mittel 
I-n(L+ Lı _L; 1,1 
H x X2 X3 X4 Xn 


das quadratische Mittel 
_,/1 1.2..2..2,.2 2 
Q= ara rRatRg rt... tx) 


Das meistgenutzte arithmetische Mittel ist sehr empfindlich gegen Ausreißer; 
d.h. bei kleinen Zahlenreihen wird das arithmetische Mittel durch einen zufäl- 
lig großen Wert verfälscht. In diesem Fall sollte man lieber den Median (vgl. 
Programm 26) nehmen. 


Das geometrische Mittel ist stets bei prozentualen Größen zu nehmen. Wird 
z.B. eine Gehaltserhöhung von 3% erst zum 1. Juni ausgezahlt, so ist die effek- 
tive Gehaltserhöhung nur 


(19 1.837 - 1): 100 %=1.7% 


Hier findet man häufig falsche Angaben, auch in angesehenen Tageszeitungen. 


Fährt ein Auto eine Strecke bergauf mit 30 km/h und bergab mit 90 km/h, so 
ist die Durchschnittgeschwindigkeit keineswegs 60 km/h, sondern nur das har- 
monische Mittel von 45 km/h. 


Für das arithmetische, geometrische und harmonische Mittel gilt die Unglei- 
chung von Cauchy 


A>G>H 
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Das quadratische Mittel hat große Bedeutung in der Statistik, da die Standard- 
abweichung einer Meßreihe gleich dem quadratischen Mittel der einzelnen Ab- 
weichungen ist. 


Zum folgenden Programm 


Das Programm berechnet für beliebige Zahlenreihen das arithmetische, geome- 
trische, harmonische und quadratische Mittel nach den oben angegebenen 
Formeln. 


Als Beispiel wurde die Zahlenreihe: 
19, 22,20, 18, 19, 23, 17, 21,25, 17 
eingegeben. Es ergeben sich folgende gerundete Mittel: 


arithmetisches M.: 20.10 
geometrisches M.: 19.95 
harmonisches M.: 19.80 
quadratisches M.: 20.26 
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REM MITTELWERTE 


HOME : 
PRINT " 
PRINT " MITTELWERTE " 
PRINT " 
NORMAL : PRINT 


INVERSE 


HEN PRINT =: PRINT "HARMONISCHES MITTEL 


INPUT "WIEVIELE ZAHLEN ? "}N 
PRINT 
A=BG6= 1 
H=8:0=9 
FOR I =1 TON 
PRINT I5".TE ZAHL "32 INPUT X 
A=SA+X 
G=-6#+X 
IFX=9T 
NICHT DEFINIERT"5 CHR$ (7): END 
H=H+1/X 
a=OrX#X 
NEXT I 
PRINT 
A=A/N 
G=-6* (1/N 
H=N/H 
Q= SIR M/N 
PRINT "ARITHMETISCHES MITTEL= "JA 
PRINT "GEOMETRISCHES MITTEL = "376 
PRINT "HARMONISCHES MITTEL = ";H 
PRINT "QUADRATISCHES MITTEL = ";Q 
END 
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MITTELWERTE 


WIEVIELE ZAHLEN ? 1% 


ZAHL 719 
ZAHL ?22 
ZAHL 720 
ZAHL 718 
ZAHL 717 
ZAHL 723 
ZAHL 717 
ZAHL ?21 
ZAHL 725 
E ZAHL ?17 


- ION UPUN 
+4 4-4 444-4 +4 
mmmmmmmmm 


nn ® 


ARITHMETISCHES MITTEL 
GEOMETRISCHES MITTEL 
HARMONISCHES MITTEL 

QUADRATISCHES MITTEL 


20.1 

19.97481892 
19.8812891 
29, 2558634 
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26. STATISTISCHE MITTELWERTE 


Neben dem arithmetischen Mittel werden in der Statistik auch noch andere 
Mittelwerte betrachtet. 


Der Modus ist der häufigste Wert einer Stichprobe oder einer Meßwertreihe. 
Für Durchschnittsiohnangaben ist sicher der Modus aussagekräftiger als der 
Mittelwert, der durch Bezieher sehr großer Einkommen verfälscht wird. 


Denkt man sich alle Stichprobenwerte der Größe nach geordnet, so stellt das 
in der Mitte liegende Element den Median dar. Ist die Zahl der Werte gerade, 
so wählt man den Mittelwert der beiden mittleren Elemente zum Median; 
z.B. beim Stichprobenumfang 100 ist der Median das Mittel aus dem 50. und 
51. Stichprobenwert. Der Median hat sehr große Bedeutung beim Testen auf 
vorgegebene Verteilungen (siehe z.B. [16]). Bei einer geringen Zahl von Stich- 
probenwerten ist das arithmetische Mittel besonders empfindlich gegen 
“Ausreißer”’. So verwendet die Stiftung Warentest als mittleren Preis eines 
Gerätes den Median, wenn weniger als 5 Geräte gekauft werden. 


Weitere Parameter, die die Verteilung von Zahlen kennzeichnen, sind die 
Varianz, die Schiefe und der Exzeß. 


Sind x; (1 <=i<=n) die Stichprobenwerte, so stellt 
(Zxj)2 
V- Sen 
n—1 
die Varianz dar. Sie ist ein Maß für die Streuung der Werte um ihren Mittel- 
wert. Die Quadratwurzel aus der Varianz wird Standardabweichung 0 ge- 
nannt. Im Falle einer Normalverteilung gibt letztere Zahl den Bereich um 
den Mittelwert an, in dem sich 65% aller Werte befinden. 


Die Schiefe wird berechnet aus 
. ‚2 . .\3 
.. (>: 22 Di, 2) )r 


n n 
Sie ist ein Maß für die Symmetrie der Verteilung. 


Der Exzeß 
.4 3 2 2 :\4 
en (> 4 Ex? Zxi ,, Exi (2) 4 2=) ). 12a 


n n n n n 
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(auch Wölbung genannt) ist ein Maß für die Steilheit der Kurve. Positiver 
Exzeß zeigt, daß die Verteilungskurve steiler als die Normalverteilung verläuft, 
negative Werte zeigen einen flacheren Verlauf. Für die Normalverteilung ver- 
schwindet sowohl Schiefe wie Exzeß. 


Zum folgenden Programm 


Gibt man die Stichprobenwerte in Form von DATA-Werten ein, so werden 
die oben genannten statistischen Kennwerte nach den angegebenen Formeln 
berechnet. Existieren zwei Werte mit gleich großer Häufigkeit, so wird der 
kleinere als Modus gewählt. 


Für die Stichprobe 
(18, 20, 21, 19, 22,23, 17,21, 18, 24) 


können die sich ergebenden Maßzahlen dem Programmausdruck entnommen 


werden. 
Die Schiefe ist klein, da genau 5 Werte unter und 5 Werte über dem Mittel von 


20.3 liegen. 
Der Exzeß ist negativ, dies zeigt eine flache Verteilung an. 
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REM STATISTISCHE MITTELWERTE 


HOME : INVERSE 


PRINT " " 
PRINT " STATISTISCHE MITTELWERTE " 
PRINT " " 
NORMAL : PRINT 

READ N 


DIM X{N) ,H(N) 


M=9:5=-:0:T=-9U=9 
FOR I=1 TON 


REM BUBBLESORT 
F=ei1l 

FOR I=1T0N-i 

IF XD < = XtI + 1) THEN 288 
H= XID:XID = X (I + D 
X(ıI+1)=HiF=98 

NEXT I 

IF F = 8 THEN 238 

PRINT "MITTELWERT=";j: HTAB 22: PRINT A 
v=S8/N-A*2 

PRINT "VARIANZ=";j: HTAB 22: PRINT V 
SB = SOR (V) 

PRINT "STANDARDABREICHUNG="}: HTAB 22: PRINT SG 
SF=I(T/N-S#AxS/N+2x (Ar SD) / (80 * 
3) 

PRINT "SCHIEFE=")5 HTAB 22: PRINT SF 


E= W/N-A#AHT/N+E*SH(A*D /N- 
ZH (AA) /ISE*HM -3 
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428 PRINT "EXZESS(WDELBUNG)="j: HTAB 22: PRINT E 


498 IFN/2< > INT (N/ 2) THENME = XUN+ 1 / 
2): GOTO 458 

a8 ME= (XKIN/ 2) +X(N/2+1)) /2 

468 PRINT "MEDIAN="5: HTAB 22: PRINT ME 

478: 

488 REM FREQUENZZAEHLUNG 

4985 FÜR I =2 TON 

22 IF X < > XI - 12) THEN 528 

Sig HD = HD +1 

522 NEXT I 

SöK=9 

42 FOR I =2TON 

8 IFH{D > HI - DD THENK=] 

5682 NEXT I 

378 PRINT "MODUS="5: HTAB 22 

s89 IFK< > BTHEN PRINT X(K): END 

378 PRINT "NICHT DEFINIERT" 

&0 END 

&iR DATA 19 

620 DATA 18,29,21,19,22,23,17,21,18,24 


STATISTISCHE MITTELWERTE 


MITTELWERT= 20.3 
VARIANZ= 4.88999997 
STANDARDABWEICHUNG= 2.19317121 
SCHIEFE= „127493993 
EXZESS (WDELBUNG) = -1.15515941 
MEDIAN= 20.5 
MODUS= 21 
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27. PYTHAGOREISCHE ZAHLEN 


Ganzzahlige Lösungen # @ der Diophantischen Gleichung 
x2 + y2 = z2 


heißen pythagoreische Zahlen. Ein Dreieck ist genau dann rechtwinklig, wenn 
seine Seitenlängen pythagoreische Zahlen sind. 


Setze man zwei solcher rechtwinkeligen Dreiecke mit gleicher Kathete zusam- 
men, so erhält man ein Dreieck, dessen Höhe ganzzahlig ist. Die Seitenlängen 
eines solchen Dreiecks nennt man auch Heronsche Zahlen, da sie rationale 
Werte der Heronschen Flächenformel liefern. 


Bekannte pythagoreische Tripel sind 
(3, 4, 5) 
(8, 6, 10) 
(5, 12, 13) 
(9, 12, 15) 


Pythagoreische Zahlen können durch die Formeln 

y = 2uv 

z = u2 + v2 
erzeugt werden; dabei durchlaufen u und v alle natürlichen Zahlen. Die 
obengenannten Tripel erhält man für 


u=2,v=1l 
u=3, v=l 
u=3,v=2 
u=4 v=2 


Wie man an den Beispielen sieht, sind auch stets alle Vielfachen eines solchen 
Tripels wieder pythagoreisch. Im folgenden Programm werden daher nur tei- 
lerfremde Tripel erzeugt. 


Zum folgenden Programm 
Das folgende Programm berechnet nach dem oben angegebenen Verfahren be- 
liebig viele pythagoreische Zahlen, wobei Vielfache anderer Tripel ausgeschlos- 


sen werden. 
Die ersten 20 vom Programm erzeugten Zahlentripel können dem Programm- 


ausdruck entnommen werden. 
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18 REM PYTHAGDREISCHE ZAHLEN 

2B : 

38 HOME : INVERSE 

405 

328 PRINT " " 
68 PRINT " PYTHABOREISCHE ZAHLEN " 
78 PRINT " u 
88 NORMAL : PRINT 

98 INPUT "WEITER <RETURN>"IT$ 

188 HOME 

19 U= 2:7 = 1 

128 GOSUB 4298 


2498 IF(Z-1) /18< > INT ((Z - 1) / 190) THEN 328 


2508 PRINT 

268 PRINT "WEITER <RETURN> ODER <E>NDE ?"3 
278 BET A$: IF A$ = "" THEN 270 

2588 IF At = "E" THEN END 


298 HOME 

SßR GOSUB 428 
sie: 
S2BEV=eV-2 


Sss@ IFV< =98THENU=U + 1: GOTO 149 

Sad vi vMil2 = U 

S5@ N= INT (U2 / v1) 

Ss68 Ui =U2-N * Vi 

s’@ IFUlI< >98 THEN U2 = Vi:Vi = Ul: G0OTO 358 
s8@2 IF Vi< > 1 THEN 329 

35978 GOTO 159 

408 END 
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420 PRINT "PYTHAGOREISCHE ZAHLEN 
458 PRINT "--------------------- 
448 PRINT 

458 PRINT "NR. A B 

468 PRINT 

478 RETURN 


PYTHAGDREISCHE ZAHLEN 


WEITER <RETURN> 
PYTHAGDREISCHE ZAHLEN 


= 
a 
2 
oJ 
OD 


l 3 4 J 

2 = 12 13 
3 7 24 25 
4 15 8 17 
=) 9 AD 41 
& 21 29 29 
7 11 &R 61 
8 35 12 37 
9 13 84 85 
12 33 36 65 


WEITER <RETURN> ODER <E>NDE ? 
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28. KOORDINATEN BESONDERER DREIECKSPUNKTE 


Sind p1, p2, p3 die Ortsvektoren der Dreiecks-Eckpunkte, so ergeben sich die 
Ortsvektoren der besonderen Dreieckspunkte aus folgenden Formeln: 


Schwerpunkt 
_ Pitp2+p3 
3 


Inkreismittelpunkt 
_ a1Bl + a2p2 + a3p3 
a ta2 +a3 


Höhenschnittpunkt 
_ t1p1 +t2p2 + t3p3 
t1t2t3 


Umkreismittelpunkt 
_ °151B1 + c2s2p2 + c3s3 p3 
2 s1s2s3 


dabei stellen aj, a2, a3 die Seitenlängen, 
s1, 52, s3 die Sinuswerte, 
c1, c2, c3 die Cosinuswerte, 
t1, t2,t3 die Tangenswerte 


des Dreiecks dar. 


Zum folgenden Programm 


Das Programm berechnet die 4 besonderen Punkte eines Dreiecks nach den 
oben angegebenen Vektorgleichungen. 
Als Programmbeispiel wird das räumliche Dreieck 


(5101-4), (-21613), (61-31 --5) 
behandelt. Es ergeben sich folgende Werte, 


Schwerpunkt: 
(3111-2) 
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Inkreismittelpunkt: 
(4.60337933 | -0.50521699 | —3.60337933) 


Höhenschnittpunkt: 
(14.6 | 17.9666654 | —13.56) 


Umkreismittelpunkt: 
(—2.8 | -7.03333264 | 3.8) 


Zur Kontrolle kann man den Satz anwenden, daß der Schwerpunkt.die Strecke 
Umkreismittelpunkt — Höhenschnittpunkt im Verhältnis 1:2 teilt 


s= 1(2.u+h) 
3 


Dies ist, wie man leicht nachrechnet, erfüllt. 
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10 REM KOORDINATEN V.BES. DREIECKSPUNKTEN 

30 HOME : INVERSE 

48 PRINT " " 

58 PRINT " KOORDINATEN BESONDERER " 

&8 PRINT " n 

78 PRINT " DREIECKSPUNKTE " 

80 PRINT " " 

90 NORMAL : PRINT 

100 PRINT "GIB KOORD.DES 1.PUNKTS EIN "3: INPUT X(1), 
X(2),X(3) 

110 PRINT "GIB KOORD.DES 2.PUNKTS EIN "$: INPUT Y(1), 
Y(2),Y(3) 

120 PRINT "GIB KOORD.DES 3.PUNKTS EIN "j: INPUT Z(1), 
2(2),2(3) 

138 PRINT 

140 : 

150 REM SEITENVEKTOREN 

168 Pi1,1) = Yil) - Z(1):PA1,2) 


Ya) - ZIZI:P(Ll,S) 


YiS) - ZI) 

179 Pt2,1) = Z(1) - Xt1L):P(2,2) = Z(2) - Xt2):P12,5) = 
Z(3I - XI) 

188 P(3,1) = XL) - YEODEPIS,ZI = Xt2) - YIZJSPIS,S) = 
(3) - YO) 

1980 : 


2290 REM SEITENLAENGEN 
218 FOR I1=1T1T035 
2285=9 

238 FOR J=1T03 
24985=5 + P(1,J * 2 

258 NEXT J 

2608 AI) = SOR (9) 

278 NEXT I 

289 : 

298 REM COSINUSWERTE 
320C=d 

318 FORI=1T03S 
328C=cC + Pi2,I) *# P(1,D 
339 NEXT I 

348 C(3) = -C / (AU) #* ADII:C= 9 
3598 FRRI1=-1T03 
38C=cC + P(3,D * P(1,D 
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CD) = -C/ AU) «ADI:C= 9 
FOR I = 1703 

C=C +P(3,I) * P(2,1 
NEXT I 

cl = -C/ AMD *AB)) 


REM SINUS- U. TANGENSWERTE 

FOR 1=1T7T03 
S{I) = SOR (1 - Ci“ 2 

IF C(I) = 2 THEN T{DD = 1E3@: GOTD 478 
TID =Sım /CMD 

NEXT I 


PRINT "SCHWERPUNKT: " 

FOR I=1T03 

PRINT (X{D + Y{D + ZIDI / 3" "5 
NEXT I: PRINT © PRINT 


PRINT "INKREISMITTELPUNKT: " 

U = Al) + AD) + A) 

FOR I = 1703 

PRINT (All) * XD + AD) # VID + AD * ZUM) / 
Us N ng 

NEXT I: PRINT : PRINT 


T=- TU TD * TOD 

IF T= 8 THEN 898 

PRINT "HOEHENSCHNITTFUNKTE " 

FOR IT=1T0 5 

PRINT (T(1) # XD + TI2) * Y{D + T{D) * ZUI)) / 
T3 n "; 


NEXT I: PRINT 2 PRINT 


5 = 2 # S(1) # S(2) * S(3) 
IF S = 8 THEN 759 
PRINT "UMKREISMITTELPUNKT: " 
FOR I = 1 T03 
PRINT (Ct1) * St1) *# XI) + C(2) # S(2) * Y(I) + 
stgH2 #cC(I) * ZI) 7 Sg" 9 
NEXT I: PRINT : PRINT 
END 
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KOORDINATEN BESONDERER 


DREIECKSPUNKTE 


GIB KOORD.DES 1.PUNKTS EIN ?5,2,-4 
GIB KOORD.DES 2.PUNKTS EIN ?-2,6,3 
GIB KOORD.DES 3.PUNKTS EIN ?&,-3,-5 


SCHWERPUNKT: 
31-2 


INKREISMITTELPUNKT: 
4.608337933 -.508521699 -3.605337933 


HOEHENSCHNITTPUNKT: 
14.5999994 17.0666654 -13.59999796 


UMKREISMITTELPUNKT: 
-2. 79999992 -7.03333264 3.79999977 
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29. POLYGON-BERECHNUNG 


Das folgende Programm berechnet Flächeninhalt und Umfang eines Vielecks 
bei bekannten Koordinaten der Eckpunkte. Der Flächeninhalt wird nach der 
Formel, vgl. [4] 


A=-[xı tra) (yı-vatlXatxgllya-v3)t..-. 
+ (Xn +x1) (Yn - yı)]/2 
berechnet. 
Da der Abstand des Punktes (xj+1 | yj+1) vom Punkt (x; |y;) nach Pythagoras 
gleich 
d=vkxirı - xi)2 + (yirı - vi)? 
ist, ergibt sich der Umfang des Vielecks durch Aufsummieren der Abstände 


der einzelnen Eckpunkte. 


Der Schwerpunkt ergibt sich aus 


ın 

x = —— 2%, (xiyviri — xXiriyi) (Xi + xi+1) mit Xn+1 = %1 
6A ie] 
1 n 

Ys = —— 2%. Ixiyirl - xir1vi) (vi + vi+1) Yn+1=y1 
6A i=1 


Zum folgenden Programm 
Das Programm berechnet für ein beliebiges Vieleck Umfang und Flächeninhalt. 
Bei der Eingabe der Koordinaten muß darauf geachtet werden, daß der Um- 
laufsinn eingehalten wird, da die Formel nicht für "überschlagene”’ Vielecke 
gilt. 
Als Programmbeispiel wird das Sechseck 

(4 12), (10 16), (1219), (8111), (315), (013) 
gewählt. Es ergibt sich 


Flächeninhalt = 67 
Umfang = 35.5551423 


Schwerpunkt = (7.35323384 | 5.20895523) 


153 


"Visieren von Winkeln mit Hilfe des Jakobstabs. Holzschnitt aus dem 16. Jahrhundert‘ 


154 


REM POLYBON-BERECHNUNG 


HOME : INVERSE 
PRINT " u 
PRINT " POLYGON-BERECHNUNG " 
PRINT " " 
NORMAL : PRINT 
INPUT "WIEVIELE ECKPUNKTE ? "}N 
DIM X(N + DI,Y(N + D 
PRINT 
PRINT "GIB KOORDINATEN DER ECKPUNKTE" 
PRINT "IM GLEICHEN UMLAUFSSINN EIN !" 
PRINT 
FOR I=1i TON 
PRINT "XtHs rg) Yang) 
INPUT X(DD,Y(D 
NEXT I 
XN+1) = X) 
Y{iN+ 1) = Y(l) 


REM FLAECHENBERECHNUNG 
F=-9 
FOR I=1 TON 
F=F +XtI *Y{iI +1) - XI +1) * Y(D 
NEXT I 
F= ABS (F) /2 


REM UMFANGSBERECHNUNG 
U=98 
FOR I=1T0ON 
U=U+r SaR (XI +1) - XD) r* 2 + (YıI+ 1) - 
mr 2 
NEXT I 


REM SCHWERPUNKTSKOORDINATEN 
XS = 8:Y5 = 9 

FOR I =1 TON 
5 =XS5 + (XII) # YII+ UV) - X I+ DVD) # (X 
(ID +X(I+ VD) 
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MO YS=YS + (KDD A YIL+ 1) - XL + DD * VID) * (Y 
(ID +YiI+D) 

420 NEXT I 

430 X3 = XS / (6 + F} 

449 YS = 5 / (6 * FF) 

450 : 

468 REM AUSGABE 

470 PRINT 

480 PRINT "FLAECHENINHALT="j: HTAB 16: PRINT F 

498 PRINT "UMFANG="}: HTAB 16: PRINT U 

5088 PRINT "SCHWERPUNKT:"}: HTAB 16 

510 PRINT "(H3xG5" , "pyag)" 

526 END 


POLYGON-BERECHNUNG 


WIEVIELE ECKPUNKTE ? & 


GIB KOORDINATEN DER ECKPUNKTE 
IM GLEICHEN UMLAUFSSINN EIN ! 


Xt1),Yt1) 74,2 
X(2),Y(2) ?19,8 
X{3),Y(D) 712,9 
X(4),Y(4) 78,11 
Xt5),Y(5) 73,5 
X(6),Yib) 79,3 


FLAECHENINHALT=67 


UMFANG= 35.5551423 
SCHWERPUNKT: (7.35323384 „ 5. 20895523) 
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30. DETERMINANTE 


Determinanten haben zahlreiche Anwendungen: Mit ihrer Hilfe kann geprüft 
werden, ob 


Vektoren linear abhängig sind 

ein lineares Gleichungssystem lösbar ist 

zwei Polynome gemeinsame Nullstellen haben 

— ein Extremum einer Funktion mehrerer Variablen vorliegt 


USW. 


Die Resultante der Polynome 
fix) =x3+2x2 + 3x +5 
g(x) = 6x2 +8x +9 


12 35 dd 
01 2 35 
Def = 6 8 9 06 © 
06 8 9 9 
06066 8 9 


z.B. ist genau dann Null, wenn f und g mindestens eine gemeinsame Null- 
stelle haben. Hier gilt Det (f,g) = 2873, somit haben f und g keine gemeinsa- 
men Nullstellen. Ersetzt man g durch die Ableitung f’ von f, so gibt die De- 
terminante Auskunft, ob f mehrfache Nullstellen hat. 


Erwähnenswert ist auch die Cramersche Regel, die erlaubt einzelne Unbekann- 
te eines Gleichungssystems zu finden, ohne das ganze System lösen zu müssen. 
Ein Programm dazu findet sich in [15]. 


Zum folgenden Programm 


Das Programm berechnet die Determinante durch fortgesetztes Entwickeln 


nach der 1. Spalte. 
Dazu werden die Zeilen 2 bis n mit dem Element a1] # @ multipliziert. Da- 
durch vervielfacht sich zunächst der Wert der Determinante um 


n—1 
a1 
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Der entsprechende Faktor wird sodann abdividiert. Subtrahiert man von jeder 
der Zeilen 2 bis n das 

aj1-fache von Zeile 1 
so wird in der 1. Spalte der Vektor 

(11818 |.... 0) 


erzeugt. 


Entwickelt man nun nach der 1. Spalte, so verbleibt eine Determinante der 
Ordnung n—1. 

Setzt man das Verfahren fort, so verringert sich die Ordnung der Determinan- 
te bei jedem Schritt, bis sie nur noch zweireihig ist. Im letzteren Fall kann der 
Wert leicht ermittelt werden. 

Die Determinante wird im Programm in Form von DATA-Werten eingelesen. 


Als Beispiel wird die Determinante 


12 345% 
2 293 45% 
3393 45% 
a4 4.45 6 
555 55% 
666666 6 


berechnet. Das Programm liefert den Wert —6. 
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REM DETERMINANTE 


HOME : INVERSE 
PRINT " n 
PRINT " DETERMINANTE " 
PRINT " u 
NORMAL : PRINT 
READ N: REM ORDNUNG DER DETERMINANTE 
DIM A(N,N),B{N,N) 
PRINT "DETERMINANTE" 
PRINT 
FORI=1TON 
FOR J=1TON 
READ AI,J): PRINT AUI,Dg" "; 
NEXT J: PRINT 
NEXT I 
D = 1: REM DETERMINANTE 
M = N: REM ZAEHLER 
REM PIVOTSUCHE 
P=9 
IF ABS (Al1,1)) > IE - 5 THEN 488 
FORK=2 TOM 
IF ABS (Atl,Kd) < = ABS (Ali,K - 1)) THEN 288 
P= Atll,Kd!L=K 
NEXT K 
IF P= 69 THEN D = 8: GOTO 418 
REM SPALTENTAUSCH 
FOR J=1 TOM 
H=Al=L) 
AlJ,L) = AlJ,1) 
Ald,i) =H 
NEXT J 
D= -D 
REM ENTWICKLUNG NACH DER 1.SPALTE 
FOR J=2 TOM 
FORK=2 TOM 
AtJ,K) = Atl,1) * Ald,K) - Ald,1) * Ali,K) 
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439 NEXT K 

440 NEXT J 

458 T= SEN (All,1)) * ABS (Atl,1) * IM - 2) 

48D=D/T 

470 FOR J = 2 T0 

488 FORK=2T0 

49 AU -1L,K-1 

598 NEXT K 

518 NEXT J 

520 : 

SSoM=M-1 

540 IF M > 2 THEN 229 

559g : 

569 REM ZWEIREIHIGE DETERMINANTE 

578D=D* (All,i) * Al2,2) - All,2) * A(2,1)) 

580 : 

598 REM AUSGABE 

6ö0 PRINT 

618 PRINT "DETERMINANTE = "}D 

620 END 

638 : 

6Ad DATA & 

6559 DATA 1,2 

668 DATA 2,2 

670 DATA 3,3 
4,4 
5,5 
6,6 


x = 


) = AlJ,K) 


680 DATA 
698% DATA 
80 DATA 


DETERMINANTE 


DETERMINANTE 


Oo UP GN 
rOABUND 
on (A Ca A 
rUülbbn»5 
ramnad u U 
eo cc von 


DETERMINANTE = -& 
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31. MÜLLER-ITERATION 


Die Iteration ist eine Programmiermethode, bei der durch zahlreiche Iterations- 
schritte ein gesuchter Punkt immer näher eingekreist und schließlich mit der 
vorgegebenen Genauigkeit bestimmt wird. 


Iterationsverfahren werden hauptsächlich zur Lösung von nichtlinearen Glei- 
chungen genutzt, da z.B. Nullstellen von Polynomen vom Grad >4 nicht mehr 
formelmäßig berechnet werden können. Iterationen finden jedoch auch An- 
wendung bei der Lösung von großen linearen Gleichungssystemen (z.B. beim 
Gauß-Seidel-Verfahren). 


Ein einfaches Verfahren ist die sog. Fixpunkt-Iteration 
Xn+1 ” F (Xn) 


Dabei wird der Graph der Funktion F{x} mit der Geraden y = x zum Schnitt 
gebracht. 


Bei der Regula falsi 
Xn — Xn—1 


Xn+1 = Xn — Fix) ———— — 
n+1 n n Fxn) — fin) 


wird die Nullstelle der Verbindungsgerade zweier Graphenpunkte von f(x) be- 
stimmt. 


Bessere Konvergenz erzielt das Newton-Verfahren 
f(xn) 
Pen) 
bei dem jeweils die Nullstelle der Tangente bestimmt wird. Die Newton-Itera- 
tion läßt sich auf nichtlineare Systeme und somit auch auf komplexe Glei- 


chungen ausdehnen (siehe Programm 33). Programme zu den erwähnten Itera- 
tionsverfahren finden sich z.B. in [15] und [17]. 


Die Müller-Iteration (1956) (in der englischsprachigen Literatur Muller ge- 
nannt) kann als Verallgemeinerung der Regula falsi aufgefaßt werden. Im Ge- 
gensatz zu dieser wird durch 3 Kurvenpunkte eine Parabel gelegt und deren 
Nullstelle bestimmt. Die Bestimmung der Parabelgleichung erfolgt mittels 
Interpolation. 


Xn+1 ” Xn 
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Sind x1, x2 und x3 die 3 Kurvenpunkte bzw. die Startwerte, SO berechnet 
man den nächsten Iterationswert aus 


x=x3+1(x3 —x2) 


Der Faktor | berechnet sich aus dem Quotienten 
— 2Df(x3) 


G + /62--4DCfix3) 


Die Konstanten ergeben sich aus 
G=L2f(x1) — D2 f(x2) + (L+D) flx3) 


C=L(Lflx1) - Dflx2) + fixg)) 


mit 
= x3—-X2 
2X] 
_x%3-%1 
2 %1 


Das Vorzeichen im Nenner von | wird aus numerischen Gründen so gewählt, 
daß der Nenner möglichst groß wird. 


Wegen der im Nenner auftretenden Wurzel könnte das Verfahren auch ins 
Komplexe gehen. Da die komplexe Newton-Iteration effektiver ist, ist das 
hier angegebene Müller-Verfahren auf das Reelle beschränkt. 


Zum folgenden Programm 


Die zu lösende nichtlineare Gleichung ist im Programm in Zeile 150 zu defi- 
nieren. Der Startwert, die Abbruchgenauigkeit und die max. Zahl der Iteratio- 
nen werden über INPUT eingegeben. 

Liegt der Startwert in der Nähe der gesuchten Nullstelle, so ist die Konver- 
genz gut. Bei Funktionen, die keine reellen Nullstellen haben, kann die Itera- 
tion manchmal zwischen zwei Werten hin- und herpendeln und so Konvergenz 
vortäuschen. Zur Kontrolle wird daher im Programm noch der Funktionswert 
des letzten Iterationsschritts mit ausgedruckt. 


Als Programmbeispiel wird das kubische Polynom 
x3-x—-1 

behandelt. Beim Startwert 1 ergibt sich in 5 Iterationsschritten die Nullstelle 
x = 1.32471796 

auf Maschinengenauigkeit. Der zugehörige Funktionswert ist 
f(x) = -4.65661287 . 19-10 
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REM MUELLER-INTERATION 


HOME : INVERSE 


PRINT " " 
PRINT " MUELLER-INTERATION " 
PRINT " " 
PRINT" FiX) = Kr 
PRINT " " 


NORMAL =: PRINT 


REM DEFINIEREN DER FUNKTION IN 158 

DEF FNFIX) =S X 3 - X - 1 

INPUT "STARTWERT 2 I,X0 

PRINT : INPUT "ABBRUCHGENAUIGKEIT ? SE 
PRINT : INPUT "ZAHL DER INTERATIONEN ? ";N 
PRINT 


REM STARTWERTE 


K=1:D= .5 

X53 = X: PRINT X8 
1 =X3-D 

2 =X5 +D 


REM MUELLER-INTERPOLATION 

IF X2 = Xi THEN X2 = X2 + IE - 8 
Li = (X3 - X2) / X2- XD 
Di = (X3S-X1) / X2- XD 

IF K > 1 THEN 329 


Ei= FNF(D 

E2 = FNF(X2) 

ES = FNF(XS 

Ai =Li”’2#+E1 -Di*2#E2 + (Li + Di) * ES 
Ci =Li* (LL*EI - DI*+E2 +ES) 


B=Al*2-4#*DI CI *ES 
REM INTERATION 
IFB< 8 THENB= 9 
IF AI < 8 THEN Al = Al - SOR (BD) 
IF AL > 8 THEN Al = Al + SOR (B) 
IF ABS (Al) + ABS (B) = BD THENAl =4# Di * E3 
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420 IF Al =9 THEN Al = IE - B 
4328 L= -2 DI #*E3 / A 
AN = XI + L*# (XI - X2 

458 PRINT X 

460 : 


472 REM PRUEFEN AUF KONVERGENZ 

486 IF ABS (X-XD <E* ABS (X5) THEN 5978 

499 IF K< N THEN 520 

S0808 PRINT 

s18 PRINT "KEINE KONVERGENZ NACH'"SNS"INTERATIONEN": END 


S2DK=SKrI 

338 X1 = X2 

40 X2 = X3 

Sa XS = X 

368 E1 = EZ!E2 = E3 

578 GOTO 268 

S8B : 

398 PRINT 

689 PRINT "FUNKTIONSWERT = "5 FN F(X3) 
610 END 


MUELLER-INTERATION 


FIX = X"S-X-1 


STARTWERT ?»1 
ABBRUCHGENAUIGKEIT ? 1E-8 
ZAHL DER INTERATIONEN ? 59 

l 

1.51344432 

1.32456943 

1.352471888 

1.52471796 

1.324717986 


FUNKTIONSWERT = -4. 55641287E-10 


164 


32. MÜLLER-ITERATION (ZWEIDIMENSIONAL) 


Die Müller-Iteration (Programm 31) läßt sich auch auf Funktionen 
fix, y) 
zweier Variablen ausdehnen. 


Dabei wird die x-y-Ebene in einfachen Müller-Iterationsschritten rechteckig 
in x- und y-Richtung durchlaufen 


x=x3+11 (x3 —x2) 
y=sy3+tla(y3—y2) 


Das Verfahren bricht ab, wenn die Änderung der x- bzw. y-Werte die vorgege- 
bene Abbruchgenauigkeit unterschreitet. 


Zum folgenden Programm 


Startwerte und Abbruchgenauigkeit werden mittels INPUT eingegeben. Das 
Programm ist gegen das Verschwinden des Nenners gesichert, jedoch nicht das 
Überschreiten des Wertebereichs (Overflow). 


Als Beispiel wird die Funktion 
fix,y)=eY-x2-2 
behandelt. 


Für die Startwerte x=1, y=1 erhält man in 2 Schritten die Lösung 
x = 0.847515091,y = 1 


Analog liefern die Startwerte x=1, y=® die Lösung 
x = -1.15802104, y = 1.20627398 


Da die Funktion unendlich viele Nullstellen hat, erhält man für verschiedene 
Startwerte wechselnde Lösungen. Sie erfüllen die Gleichung 


y=In(x2 +2) 
die man erhält, wenn man die gegebene Funktionsgleichung nach y auflöst. 
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298 
308 
s1ß 
329 
332 
348 
35g 
Ki-], 
378 
Ss8ß 
398 
429 
419 
420 


REM MUELLER-ITERATION ZWEIDIMENSIONAL 


HOME =: INVERSE 


PRINT " 

PRINT " MUELLER-ITERATION ZWEIDIMENSIONAL 
PRINT " 

PRINT " F(X,Y) = EXP(Y) - X*2 - 2 
PRINT " 


NORMAL = PRINT 


REM FUNKTION IN ZEILE 708 EINGEBEN 


INPUT "STARTWERTE X,Y 7 "3x9, VQ 


INPUT "ABBRUCHBENAUIGKEIT ? IE 
INPUT "ZAHL DER ITERATIONEN ? ";N 
PRINT 


REM ANFANBSWERTE 
K=1:Bli = „5:B2 = Bi 


REM MUELLER-INTERPOLATION 

IF X2 = X1 THEN X2 = X2 + IE - B 
L1i=(X3-X2 7 X - X) 
Di m (X3 - X1) / (X2 - Xi) 
Y=YaX=X 

GOSUB #70:E1 = WU 
X=X2 

GOSUB 678:E2 = W 
X=X3 

GOSUB 578: E3 
GOSUB 729 
Bl=L# (XI -X2) 
X = X3 + Bi 


u 
= 


REM PRUEFEN AUF KONVERGENZ 
IF ABS (Bi) + ABS (B2) < E THEN 838 
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SONG = X 


440 Y3 = YO 
450 Y1 =Y3 - B2 
468 Y2 = YI + B2 


478 IFY2=Y1 THENY2=Y2 +1E-8 


4808 Li = (Y3 - Y2) / (v2 - Mi) 
498 Di = (Y3 - Yı) / (Y2 - YUV 
Sa Y = Yı 

S12 GOSUB 47M:E1 = W 

S20 Y=Y2 

352 GOSUB 57@:E2 = W 

Sa0 Y = Y3 

sd GOSUB #70:E3 = W 


360 GDSUB 729 

378 B2 =L * (Y3 - Y2) 

86 Y = Y3 + B2 

78 IF ABS (Bi) + ABS (B2) < E THEN 838 

688 IF K < N THEN 639 

818 PRINT 

628 PRINT "KEINE KONVERGENZ NACH "3N5" TITERATIONEN": END 


&SBK=K ri 
Ba YB=NY 
658 GOTD 219 


678 REM FUNKTIONSAUSWERTUNG 

&B8 W= EXP (Y)- X 2-2 

&98 RETURN 

788 : 

718 REM HILFSPROGRAMM 

20 Ali =Li*2 Ei - Di*2#E2 + (Li +Di) *# ES 
738C1 =Li* (LI*Ei - DIL +E2 + EI 

74585B=A1i * 2-4 Di * Ci *E3 

50 IFB<9THENB = 9 

7628 IF Al < 98 THEN AI = Al - SOR (BD 

778 IF AL >®@ THEN Al = Al + SOR (B) 

788 IF ABS {Al) + ABS (B) = B THENAL =4 # DI * EI 


778 IF Al = 8 THEN Al = IE - 8 
BG L= -2 DI *E3 / Al 

818 RETURN 

828 ı 


838 PRINT : PRINT "FUNKTIONSWERT = "; 
848 Y = Y3:X = X3: GOSUB 480 
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852 PRINT W 
B68 END 


MUELLER-ITERATION ZWEIDIMENSIONAL 
F(X,Y) = EXP(V) - X*2 - 2 

STARTWERTE X,Y 24,1 

ABBRUCHGENAUIGKEIT  ? 1E-6 

ZAHL DER ITERATIONEN ? 3@ 


1 1 
847515891 1 


FUNKTIONSWERT = -9.31322575E-18 


MUELLER-ITERATION ZWEIDIMENSIONAL 
F(X,VY) = EXP(V) - X*2 - 2 
STARTWERTE X,Y 21,8 


ABBRUCHGENAUIGKEIT ? 1E-6 
ZAHL DER ITERATIONEN ? 3@ 


1 9 
-. 279797997999 1.29627398 
-1.15802104 1.28627398 


FUNKTIONSWERT = @ 
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33. NEWTON-ITERATION IM KOMPLEXEN 


Die Newton-Iteration im Reellen 


f(Xn) 
Xn+] =Xn —- 
f’(xn) 
kann auch auf die komplexe Ebene ausgedehnt werden 
f(zn) 
Zn+1 = Zn = 
f (Zn) 


Um eine komplexe Arithmetik zu vermeiden, wird die Funktion f und ihre 
Ableitung f’ in Real- und Imaginärteil zerlegt 


flz) =u+iv 
f’(z) = ux + iuy 


Damit erhält man 
f utiv Uux +Yuy _ VUx — Uuy 
f’ Ux + iuy ux2 + uy? ux2 + uy? 

durch Erweitern mit dem Konjugiert-Komplexen des Nenners. Dabei stellen 
Ux, uUy 


die partiellen Ableitungen von u nach x bzw. y dar. Trennt man Real- und 
Imaginärteil, so erhält man die komplexe Newton-Iteration 


uux + Vuy 


Xn+1 =Xn — 
n+1 n uy2 +uy2 


Vux — UUy 


+1 = Yn - 
Yn+1 n ux2 + uy2 


Polynome können mit Hilfe der binomischen Formel in Real- und Imaginär- 
teil zerlegt werden. Z.B. gilt 


233+1 =(x+iy)3+1 
= x3 + 3ix2y — I3xy2 — iy3 +1 


Dies liefert den Realteil u = x3 — 3xy2 +1 und Imaginärteil v = 3x2 y—y3. 
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Andere Funktionen können mit Hilfe eines geeigneten Tafelwerks zerlegt wer- 
den. Nach [4] gilt z.B. 


sin z = sinx coshy + i cosx sinh y 
Dabei treten die hyperbolischen Funktionen sinh bzw. cosh auf. 


Zum folgenden Programm 


Der Realteil u bzw. Imaginärteil w der Funktion muß im Programm als Unter- 
programm in den Zeilen 360 - 370 vereinbart werden. Die partiellen Ablei- 
tungen von u nach x und y folgen in den Programmzeilen 390 - 400. 


Die Startwerte, die Abbruchgenauigkeit und die maximale Zahl der Iteratio- 
nen wird über INPUT-Anweisungen eingegeben. 


Als Programmbeispiel wird die Gleichung 
23 +1 

gelöst. Die Zerlegung in Real- und Imaginärteil ist 
u=x3—-3xy2+1, v=3x2y—y3 

die partiellen Ableitungen von u nach x bzw. y sind 
ux = 3x2 — 3y2, uy = —Öxy 


Eingabe der Startwerte x=1, y=1 liefert die Nullstelle 
0.5 + i * 0.866025403 

Entsprechend liefern die Startwerte x=1, y=—1 die Nullstelle 
0.5 — i + 0.866025493 

bzw. x=1, y=® die Nullstelle 
—0.999999999 + i = Q, 


Die exakten Nullstellen sind hier 
+ 
—1 und | v2 
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id REM NEWTON-VERFAHREN IM KOMPLEXEN 


20 : 

39 HOME : INVERSE 

48 PRINT " " 
5 PRINT " NEWTON-VERFAHREN KOMPLEX " 
68 PRINT " " 
78 NORMAL : PRINT 

88 INPUT "STARTWERTE ? "3X0, Yo 
90 INPUT "GENAUIGKEIT ? "}E 
120 INPUT "ZAHL DER ITERATIONEN ? ";N 
118 PRINT 

128 : 

130 REM ANFANGSWERTE 

140 K=9 

150 PRINT X9,Y@ 

160 : 

170 REM ITERATION 

IBK=K+i 


198 X = X0:Y = 9 

200 GOSUB 359 

2196=-U1*2 + U2*%2 

228 IF 6 = 8 THEN PRINT "SRADIENT IST NULLVEKTOR": END 


236 X=X0 + W#U2 -U* U) 7/8 

28Y=YO -(V#UL+U*U2) /6 

250 PRINT X,Y 

260 : | 

270 REM KONVERGENZTEST 

2 IF X -D72+ (WdE -YVr2<CE* 2 THEN END 


299 IFK< N 1THEN 3182 
388 PRINT "KEINE KONVERGENZ NACH ";NJ" ITERATIONEN": END 


SIE XD = X 
328 Yo = Y 
338 GODTO 189 
348 : 


358 REM FUNKTIONSUNTERPROGRAMM 
ZBAU=-XI-IHKHYNT 22H 1 
370V=3#HX*2#Y-Y”53 
38@ REM ABLEITUNGEN 

379 11 =-3+X*2-3#Y°%2 
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19 0U2= - 6X H“Y 
418 RETURN 


NEWTON-VERFAHREN KOMPLEX 


STARTWERTE ? 1,1 

GENAUIGKEIT ? 1E-5 
ZAHL DER ITERATIONEN ? 19 

1 1 

AR-1-1-1-1-1-1-1- 7 835333333 

9398691918 .8410997874 

49793299796 8662697171 

„479799711 . 8660249053 

.9 . 866825403 

URUN 


NEWTON-VERFAHREN KOMPLEX 


STARTWERTE 1,1 
GENAUIGKEIT ? 1E-5 
ZAHL DER ITERATIONEN ? 29 

1 -1 

BR -1-1-1-1-1-1-1- 94 -. BI3IIIII3 
3088691916 -.841899874 
.47935299976 -. 866269171 
„49997979911 -. 866024995 
.9 -. 8660825493 
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NEWTON-VERFAHREN KOMPLEX 


STARTWERTE 
GENAUIGKEIT 


1 
333333334 
-2. 77777777 
-1.89595185 
-1.35618883 
-1. 88555956 
-1.00653788 
-1. 00004256 
-1 
-.999979979977 


gssaoaıs aaa Ss 5 & 


? 1,9 
? 1E-5 
ZAHL DER ITERATIONEN ? 29 
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34. NEVILLE-INTERPOLATION 


Unter Interpolation versteht man Aufsuchen von Zwischenwerten einer tabel- 
lierten Funktion. Im Gegensatz dazu wird bei der Extrapolation ein Wert 
außerhalb des tabellierten Bereichs gesucht. 


Sind die Funktionswerte Z(l) an den Stellen X(l) gegeben und soll an der 
Stelle X interpoliert werden, so läßt sich die Neville-Interpolation wie folgt 
beschreiben 


268 FORK=1 TON 

278 FOR I=K+1 TON 

280 IF I=N+1 THEN 330 

298 ZU) = (ZU) » Y{K) — ZIK) » YLDI /LYIK) - YIN) 
3060 NEXTI 

318 NEXTK 


Dabei gilt 
Y)=-X-x() 


N ist die Anzahl der vorgegebenen Funktionswerte, Z(N) der gesuchte, inter- 
polierte Funktionswert. 


Die in Programmzeile 290 auftretenden Differenzen heißen dividierte Diffe- 
renzen. Mit ihrer Hilfe kann auch das Polynom vom Grad höchstens N—1 be- 
stimmen, das durch die vorgegebenen Punkte geht. Ein Programm zur Newton- 
Interpolation findet sich in [17]. 


Zum folgenden Programm 


Im Programm wird nach dem Neville-Verfahren im angegebenen Punkt inter- 
poliert. 


Folgende Funktionsstellen (xIy) seien gegeben: 
(-4 | -26), (-1 14), (812), (2 116) 
Gesucht werde der Funktionswert im Punkt x=1. 


Das Programm liefert den Funktionswert 
f1)=4 

Das interpolierende Polynom ist hier 
px)=x3+2x2-xt+2 


wie man leicht nachrechnet. 
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REM NEVILLE-INTERPOLATION 


HOME : INVERSE 


PRINT " " 
PRINT " NEVILLE-INTERPOLATION " 
PRINT " " 


NORMAL © PRINT 
INPUT "WIEVIELE STUETZSTELLEN ? ";N 
DIM X{N),Y(N),Z(N) 


PRINT 

INPUT "X-WERT 2 ."5X 
PRINT 

FOR I=1T0ON 

PRINT "yemsjz"m),yYerglg®) us 


INPUT " ? "3X(D,Y(D 
Zi = YiDEYOD = X - XD 

NEXT I 

FORK= 1 TON 

FOR I=K+1TON 

IFI=N + 1 THEN 299 
ZID = (YK #* ZID - VID * ZIK)) / (YiK) - Y(D)) 
NEXT I 

NEXT K 

PRINT 

PRINT "INTERPOL. WERT Y = "5Z(N) 
END 
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NEVILLE-INTERPOLATION 


WIEVIELE STUETZSTELLEN ? 4 


X-WERT 21 
XtD),Y(1) ? -4,-26 
X(2),Y(2) 72 -1,4 
X(3),Y(3) 20,2 
X{4),Y(4) ? 2,16 


INTERPOL. WERT Y = 4 
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Pars immobilis 


gg aD dB GB GB GG BGB BA BGB 


tr mare ae Mihusculd - Tr Aara- 

| 0 ! 3 , 2 

K IE EIN TERH y: 
Magna Rota (1 aluscula 


“Skizze der Rechenmaschine von Gottfried Leibniz (Entwurf 1673)” 
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35. GAUSS-JORDAN-VERFAHREN 


Löst man im Gleichungssystem 
a11 X1 +2412x2+....+31nXn=bi 
a21 xX1 +ta22x2+....+232n Xn”=b2 
aniX Fan2X2+....tTannXn” PBn 
die k-te Gleichung nach der Unbekannten xx auf, so erhält man 
xk = (-ak1 X1 —aK2X2 —.... —akn Xn + br) /akk 


Dies setzt voraus, daß das Diagonalelement agk # @ ist. Ist ak = @ oder sehr 
klein, so führt man eine entsprechende Spalten- oder Zeilenvertauschung 
durch. Findet sich keine Zeile oder Spalte, die ein Element # @ liefert, so ist 
die entsprechende Matrix singulär, d.h. die Determinante der Matrix ver- 
schwindet und das Gleichungssystem ist nicht eindeutig lösbar. 


Setzt man diese Unbekannte xx in alle anderen Gleichungen ein, so verschwin- 
det sie aus diesen Gleichungen, d.h. sie wird eliminiert. Führt man dies für 
alle Unbekannten durch, so bleibt in der 1. Zeile nur noch die x1, in der 2. 
noch x2 und entsprechend in der K-ten Zeile xx übrig. 


Man hat somit ein Gleichungssystem erhalten, in dem die Unbekannten nur 
noch in der Diagonale der Matrix stehen 


d11xX1 = ci 
d22x2 =Q 
d33 x3 =cz3 


dnn*n "En 


Dieses Gleichungssystem kann sehr einfach gelöst werden 


Ck 
x = ——; dek #0 
KT dkk 
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Lineare Gleichungssysteme finden zahlreiche Anwendungen bei 
— Dimensionierung elektrischer Schaltkreise 
— der Lösung von partiellen Differentialgleichungen 


— Vektorrechnung 
— Input- und Outputmodellen der Wirtschaft 


— Ausgleichsrechnung 
— elastische Verformungen von Fachwerken, Karrosserien 


USW. 
Zur Konstruktion eines Flugzeugflügels treten z.B. Gleichungssysteme mit 
mehr als 20.000 Unbekannten auf. 


Zum folgenden Programm 


Das Programm löst beliebige lineare Gleichungssysteme, bei denen die Zahl 
der Unbekannten mit der Zahl der Gleichungen übereinstimmt. 
Gibt es mehr Unbekannte als Gleichungen, so heißt das System unterbe- 
stimmt, im umgekehrten Fall überbestimmt. Zur Lösung von überbestimmten 
Gleichungssystemen, wie sie bei der linearen Regression auftreten, gibt es spe- 
zielle Verfahren wie z.B. die Methode der kleinsten Quadrate. 


Als Programmbeispiel wird das Gleichungssystem mit 6 Unbekannten behan- 
delt. 
6x1 3x2 +2x3+ XA— xX5+ x6 = 11 
3x1 - 7x2 +t4x4—2x5+ x7 = —5 
4x1 - 3x2 +6x3— x4 +2x5 + x7 = 28 
2x1 + 4x2 + 5x3 — 7x4 — 3x5 + 2x7 = -6 
—x1 + 5x2 — 4x3 +8x5 — 2x7 = 25 
3x1 + 4x3 — 2x4 + 5x5 —6x7 = —4 


Die Anzahl der Unbekannten und die Gleichungsmatrix mit rechter Seite wird 
im Programm in Form von DATA-Werten eingelesen. 


Das Programm liefert hier die exakten Lösungen 
xı=1,x2=2, x3=3, x4=4, x5=5 und x6=6. 
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REM GAUSS-JORDAN-VERFAHREN 


HOME : INVERSE 
PRINT " 


PRINT " GRAUSS-JORDAN-VERFAHREN " 


PRINT " 
NORMAL : PRINT 
PRINT 


PRINT "GLEICHUNGSSYSTEM: " 


PRINT 
READ N 
DIM A{N,N + 1),X{N) 


REM EINLESEN DES GLEICHUNGSSYTEMS 


FOR I=1T0 
FOR J=1T0 
READ A(I,J) 


a = 


+ 1 


IF AUI,JD) > =98THEN PRINT " "; 


PRINT ACI,J)5" "3 
NEXT J 

PRINT 

NEXT I 

FORK = 1 TON 

IF K= N THEN 399 

5 = Alk,K) 

M=K 
FORI=K+1TON 
IF ABS (AlI,K)) € 

G= ABS (Atı,K)) 

MI 
NEXT I 
FORI=KTON+I 

H = AtK,J) 

AlK,J) = AM, J) 

AM, = H 
NEXT J 
IF ABS (AlK,K)) < 
ER": END 
FORI=1TON 
IF I =K THEN 468 


= 5 THEN 338 


= IE - 8 THEN PRINT "SINGULA 
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F= Atll,K) / AlK,K) 

FOR J=K+IiTON+I 
AlI,J) = AllL,J) - AK, # F 
NEXT J 

NEXT I 

NEXT K 


REM AUFLOESUNG 
FOR I=1TON 
XD) = AtI,k,N + 1) / A(I,D 


NEXT I 

PRINT 

PRINT "LOESUNG: " 

PRINT 

FOR I=1 TON 

PRINT "X; 15") = "3X 
NEXT I 

END 

DATA 6 

DATA 6, -3,2,1,-1,1,11 
DATA -3,-7,0,4,-2,1,-5 
DATA 4,-3,6,-1,2,1.28 
DATA 2,4,5,-7,-3,2,-6 
DATA -1,5,-4,9,8,-2,25 
DATA 3,0,4,-2,5,-6,-4 
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GRAUSS-JORDAN-VERFAHREN 


GLEICHUNGSSYSTEM: 


6-32 1-11 1 
3-79 4-2 1-5 
4-3 6-1 2 1 28 
2445-7-3 2-6 
-1 5-4 89 8-2 25 
s 84-2 5-64 


LOESUNG: 


(1) 
X(2) 
(5) 
X(4) 
X(5) 
X(6) 


ap GN 
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36. MATRIZEN-INVERSION NACH FADDEJEW 


Eine einfache Methode zur Invertierung einer quadratischen Matrix, die zu- 
gleich noch das charakteristische Polynom liefert, ist das Verfahren von 
Faddejew [9]. 


Ist die zu invertierende Matrix Aı von der Ordnung n, so wird iterativ die 
Matrix An in n-1 Schritten über die Gleichung 


Ai+rı = Aı (Ah -ciE 


berechnet; dabei ist E die Einheitsmatrix und c; = Spur (A;)/i. Die Spur einer 
Matrix ist die Summe der Diagonalelemente. 


Die gesuchte inverse Matrix ergibt sich dann aus 


AT! = (An-1--1E) . 
n 


Der Wert cn liefert außerdem noch den Wert der Determinante 
Det(A) = (-1)nt1 cn 


Programme zu anderen Inversionsverfahren finden sich z.B. in [17]. 


Zum folgenden Programm 


Das Programm berechnet die inverse Matrix nach dem angegebenen Verfahren. 
Die zu invertierende Matrix wird in Form von DATA-Werten eingelesen. 


Als Beispiel soll die Inverse von 


5 7 6 5 

7 10 8 7 

6 8 10 9 

5 7 9 10 
berechnet werden. Das Programm liefert hier die exakte Inverse 

68 -41 —17 10 

a1 = 41 25 19 -6 

-17 19 5-3 

1 -6 -3 2 


A= 


Durch Matrizen-Multiplikation (Programm 23) kann geprüft werden, daß gilt 
AA—1 = A—1 A=E. 
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REM MATRIZENINVERSION NACH FADDEJENW 


HOME : INVERSE 
PRINT " " 
PRINT " MATRIZENINVERSION NACH FADDEJEW " 
PRINT " n 
NORMAL : PRINT 
READ N: REM ORDNUNG DER MATRIX 
DIM A(N,N) ,AL{N,N),B{N,N) ,G{N,N) ,H{N,N) 
PRINT "MATRIX: " 

FORI=1 TON 
HT = 1 

FOR J=1 TON 

READ A(I,J:AL(I,TD) = All,J) 

HTAB HT 

IF A(I,DD > = @ THEN PRINT" 9; 

PRINT A(1,J)3 
HT =HT +6 

NEXT J 

PRINT 

NEXT I 

FOR I=1 TON 

5S=9 

FOR J=1 TON 

S=8+AU,D 

NEXT J 


EN H(J,K) = AtJ,K): GOTO 349 


IFI<X >N- 1 THEN 449 
FOR J = 1 TON 

FORK = 1 TON 
G(J,K) = H(J,K) 

NEXT K 

NEXT J 
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4US= GH 

4786 FORL= 1 TON 

488 5 = 5 + Al(J,L) * HiL,K) 
490 NEXT L 

509 B(J,K) = $ 

510 NEXT K 
520 NEXT J 
539 : 

548 FOR J 
559 FOR K 
568 A(J,K) 
5786 NEXT K 
589 NEXT J 
598 NEXT I 


I TON 
I TON 
B(J,K) 


o 
a 
SQ 


619 IFC = 8 THEN PRINT : PRINT "MATRIX SINGULAER": END 


630 PRINT "INVERSE MATRIX: " 
648 FÜR I=1 TON 


=1T0N 
678 G(1,J) = GlI1,J) /C 


698 IF G(I,J) > =98THEN PRINT " "5 

788 PRINT G{1,93;" "5 

710 HT=HT +6 

728 NEXT J: PRINT 

738 NEXT I 

7D=-i(-D*N+D«+C 

75@ PRINT 

768 PRINT "DETERMINANTE 

778 END 

788 : 

7998 DATA 4 

8008 DATA 5 

818 DATA 7 
6 
3 


"sD 


828 DATA 
838 DATA 
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MATRIZENINVERSION NACH FADDEJEW 


MATRIX: 

I 7 & I 
7 18 8 7 
=) 8 18 9 
= 7 9 19 


INVERSE MATRIX: 

68 -4 -17 19 
Al 25 id -6 
-17 12 = -3 

18-6 -3 2 


DETERMINANTE = 1 
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37. CHARAKTERISTISCHES POLYNOM NACH FADDEJEW 


Das charakteristische Polynom hat genau die Eigenwerte der Matrix als Null- 
stellen. 


Eigenwerte A und Eigenvektoren x sind Lösung der Gleichung 
Ax=\x 


Eigenwerte von Matrizen geben an 

— Klassifikation von Flächen n-ter Ordnung 
— Schwingungsfrequenzen eines Systems 

— Kreiseleigenschaften von festen Körpern 
— Stationäres Verhalten von Markow-Ketten 
usw. 


Wie beim vorhergehenden Programm entstehen beim Faddejew-Verfahren [9] 
Ai+j = A1 (Ai -ciE) i=1,2,...,n—1 

die Koeffizienten ci; aus 
ci = Spur (A;)/i 


Versieht man die Werte c; mit alternierenden Vorzeichen, so erhält man die 
Koeffizienten des charakteristischen Polynoms der Matrix. Bestimmt man die 
Nullstellen des Polynoms mit Hilfe der Programme 17, 18, 31, 32 oder 33, so 
erhält man die gesuchten Eigenwerte der Matrix. 

Programme zu anderen Eigenwertverfahren, die auch gleichzeitig die Eigen- 
vektoren bestimmen, finden sich in [17]. 


Zum folgenden Programm 


Analog zum vorhergehenden Programm werden die Koeffizienten c; bestimmt 
und mit den richtigen Vorzeichen versehen. Die Matrix wird wieder in Form 


von DATA-Werten eingelesen. 
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Als Beispiel wird die Matrix 


644 1 
ja6ı [4 
Al ı 64 
1446 


behandelt. Das Programm liefert das charakteristische Polynom 
p(x) = x4 — 24x3 + 150x2 — 200x — 375 


dessen Nullstellen bereits in Programm 18 bestimmt worden sind. Die Eigen- 
werte von A sind somit 


A1=-1; %9=5; X3=-5 ; Aa = 15. 
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REM  CHARAKT.POLYNOM NACH FADDEJEN 


HOME : INVERSE 
PRINT " " 
PRINT " CHARAKT.POLYNOM NACH FADDEJEW " 
PRINT " " 
NORMAL : PRINT 
READ N: REM ORDNUNG DER MATRIX 
DIM A(N,N),AL(N,N),B{N,N) ‚G{N,N),H{N,N),C(N) 
PRINT "MATRIX: " 
FOR I = 1 TON 
HT = 1 
FOR J = 1 TON 
READ AtI,D:AltI,D = All, 
HTAB HT 
IF AU,D > = 8 THEN PRINT" "35 
PRINT A{I,JV$ 
HT =HT +6 
NEXT J 
PRINT 
NEXT I 
PRINT 
FORI=1 TON 
g=9 
FORI=S1 TON 
8 =85 + AlJ,)) 
NEXT J 
CD =5/1 


EN H(J,K) = AtlJ,K): 60TO 359 
- C(D 


IFI< >N- 1 THEN 459 
FOR J=1TON 

FOR K = | TOD N 
GtJ,K) = HUJ,K) 

NEXT K 
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430 NEXT J 

448 : 

458 FORJ = 1 TO 

468 FORK= 1 TOD 

405=6 

4880 FORL = 1 TON 

4998 5 = 5 + Al(J,L) * HiL,K) 

Ss NEXT L 

519 BiJ,k) =85 

520 NEXT K 

538 NEXT J 

549 : 

S59 FOR J 

560 FOR K 

570 AlJ,K) 

588 NEXT K 

598 NEXT J 

&Öd NEXT I 

iR: 

628 PRINT 

638 PRINT "KOEFFIZIENTEN DES CHARAKT.POLYNOMS: " 

6495 = 1 

6598 IFN/2 

668 PRINT - 

678 FORI = 

688 Lil) =5 

698 PRINT C 

708 NEXT I 

718 END 

728 : 

730 DATA 4 

748 DATA &,4 

758 DATA 4,6, 
4,1 
1,4 


N 
N 


1 TON 
1 TON 
B(J,K) 


= INT (N/ 2) TENS= -1 
5 

1 TON 
* C(I) 
(1) 


768 DATA 
778 DATA 
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CHARAKT.POLYNOM NACH FADDEJEW 


MATRIX: 

& 4 4 l 
4 6 1 4 
4 1 & 4 
1 4 4 & 


KOEFFIZIENTEN DES CHARAKT.POLYNOMS: 
l 

-24 

150 

-2028 

-375 
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38. NUMERISCHE DIFFERENTIATION 


Obwohl Funktionen —im Gegensatz zur Integration — analytisch differenziert 
werden können, ist dennoch das folgende Programm von Interesse, da es einen 


Überblick -über Funktions- und Ableitungswerte im ganzen Intervall ermög- 
licht, 


Das Programm benützt die 3-Punkte-Formeln zur Ableitung: 
—3f(x) + 4f(lx+th) — f(x+2h) 


f(x) = ah 


(f(x) — 2flx+h) + f(x+2h) 


f(x) - 2 


Da die Formeln nach links "'übergreifen‘’, ist zu beachten, daß die Funktion 
auch links vom betrachteten Intervall definiert sein muß. 


Da in BASIC die Verkettung von Funktionen implementiert ist, können die 
beiden Formeln direkt als Funktionen mittels 


DEFFN.. 
vereinbart werden. 


Dazu muß nur in Zeile 150 die zu differenzierende Funktion als 

DEF FNFix) =... 
definiert werden. 
Das Intervall und die Schrittweite h werden über eine INPUT-Anweisung ein- 
gegeben. Wegen der unvermeidlichen Rundungsfehler werden alle Funktions- 
werte im Programm auf 4 Dezimale gerundet. Dies kann, wenn gewünscht, 


leicht geändert werden. Die Schrittweite zur Ausgabe ist im Programm mit 
@.1 vorgegeben (Zeile 280), kann aber ebenfalls variiert werden. 


Zum folgenden Programm 
Als Programmbeispiel wird die Funktion 
f(x) = x2e—% + sinx 

gewählt. Das Intervall soll [1,2], die Schrittweite 8.1 sein. Wie man dem Pro- 
grammausdruck entnimmt, wächst der Funktionswert bis x = 1.6 an und 
nimmt rechts davon wieder ab. Die Funktion hat somit zwischen 1.6 und 1.7 
ein Maximum, wie es auch am Vorzeichenwechsel der 1. Ableitung erkennbar 
ist. Die 2. Ableitung ist im ganzen Intervall negativ, d.h. der Graph ist rechts- 


gekrümmt. 
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REM NUMERISCHE DIFFERENTIATION 


HOME : INVERSE 


PRINT " 


PRINT " NUMERISCHE DIFFERENTIATION " 


PRINT " 
NORMAL ; 


REM DEFINIERE FUNKTION IN 159 


PRINT 


DEF FNF(X) =X*2%* EXP(-X) + SIN (X 
PRINT "FUNKTION: "5 


PRINT "F(X)=X"2#EXP(-X) +5INIX)" 


PRINT 


REM ERSTE ABLEITUNG 


DEF FNAX = 


REM ZWEITE ABLEITUNG 
DEF FNB(X) = (FNF(X) -2* FNFiX+H + FN 
FX+2s*H) /H*2 


INPUT "INTERVALL 


INPUT "SCHRITTWEITE 


PRINT " 
x)" 


X 


PRINT %;: 


IFF > 


PRINT F}3: 


IFA > 


PRINT A;: 


IFB> 

PRINT B 
NEXT X 

END 


FIX) 


-(-3%# FNF(O +4* FNF( + 
H) - FNF(X+2xH)) / (2x H 


? "I3A,B: PRINT 


TOB+IiE-%& 


FN F(X) 
FN A(X) 
FN B(X) 
® THEN 
HTAB 7 
® THEN 
HTAB 17 
8 THEN 
HTAB 27 
8 THEN 


* 1E4 
* 1E4 
* 1E4 
PRINT 
PRINT 
PRINT 


PRINT 
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? "3HE PRINT 


F’(X) 


STEP „i 


F"3 CHRS (34) 1" 


+.5) / 164 
+.5) / 1E4 
+ .9) / 1E4 


NUMERISCHE DIFFERENTIATION 


FUNKTION: FOX)=X"2#EXP (-X) +SIN (X) 


INTERVALL 


?1,2 


SCHRITTWEITE ? „1 


SONDGTUOIPWN 


N me mu me pa pa pa me me ie me DC 
® 


FIX) 
1.2094 
1.294 
1.5658 
1.4241 
1.4688 
1.4995 
1.5164 
1.5176 
1.5094 
1.4862 
1.4596 


FIX 
-. 9196 
-,. 7848 
-,6525 
-,5157 
-,3769 
-,2374 
-,. 2989 
0372 
„1696 
„2971 
„4187 


F'M 

-1. 2854 
-1.3399 
-1.3755 
-1.3885 
-1.3867 
-1.3699 
-1.3396 
-1.29697 
-1.2431 
-1.1792 
-1.1862 
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33. NUMERISCHE INTEGRATION (TRAPEZREGEL) 


Da viele Funktionen wie 


sin (x2}, ex2, SImX 


nicht analytisch integriert werden können, ist man hier auf numerische Inte- 
gration angewiesen. Es gibt zahlreiche Integrationsmethoden. Ein bekanntes 
Verfahren ist die Trapezregel (auch Sehnentrapez-Regel) genannt: 


b 
/ fixldx =D [ffa) + 2flarh) + 2ffar2h) +. + 2b-h) + fib)] 
a 
mit der Schrittweite h = >, n Zahl der Intervallunterteilungen. 


Die Trapezregel eignet sich auch zur Extrapolation (siehe [33], [35], [36)). 
Weitere Programme zur numerischen Integration finden sich in [15] und [17]. 


Zum folgenden Programm 


Das Programm berechnet das bestimmte Integral einer Funktion nach der an- 
gegebenen Trapezformel. 


Die zu integrierende Funktion ist in Zeile 90 mittels 


DEF FNF(X) >=... 


zu definieren. Die Intervallgrenzen a, b und Anzahl n der Unterteilungen wer- 
den über INPUT eingegeben. 


Als Beispiel wird das Integral 
1 


4 
dx 
J 1+x2 
() 
mit 5@ Unterteilungen berechnet. Das Programm liefert 
3.14152599 


Vergleicht man dies mit dem exakten Wert 
a = 3.14159265 

so liefert dies den Fehler 
7:10. 
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18 REM NUMERISCHE INTEGRATION{TRAPEZREGEL) 
28 HOME : INVERSE 


38 PRINT " " 
48 PRINT " TRAPEZREGEL " 
38 PRINT " " 
62 NDRMAL © PRINT 

ww: 


88 REM DEFINIEREN DER INTERGRANDENFUNKTION 
9 DEF FNFIX =4/ (L+X«N 

188 PRINT "FUNKTION: "5 

118 PRINT "FiX)=4/(1+rX#X)" 

128 PRINT 

159 : 

148 INPUT "INTERVALLGRENZEN ? "5A,B 
150 PRINT 

168 IFA>BTHENH = MA =BEB=H 
178 INPUT "ZAHL DER INTERVALLE ? "}N 
188 PRINT 


198 : 

200 1 FN FiA) + FN FiB 

21ö H= (B-A) /N: REM SCHRITTWEITE 
220 : 

25385=98 


248 FORK=1T0O N-1 
2385=5+ FNFA+K*+H 
260 NEXT K 
I + 
I* 


2%*5 
H/2 


rn 

00 

& 
—. 4 


298 
588 PRINT "INTERGRAL = "51 
sid END 


TRAPEZREGEL 


FUNKTION: F(X)=4/{1+X#X) 
INTERVALLGRENZEN ?R,1 
ZAHL DER INTERVALLE ? 59 


INTERGRAL = 3, 14152599 
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40. HAMMING-VERFAHREN 


Bekannte Verfahren zur Integration von gewöhnlichen Differentialgleichun- 
gen sind die Verfahren von Runge-Kutta und Adams-Bashforth (siehe [33], 
[35], [37]. Programme dazu finden sich z.B. in [17]. 

Ein neueres Verfahren stammt von R.W. Hamming (1959) [12]. 


Die Differentialgleichung 
y'’=flx,y) 
mit dem Anfangswert 
y(a}) = yg mit xe [a,b] 
habe die Lösungsfunktion u(x). Im folgenden wird 
ulx;) = u, fx, u) = fi 
gesetzt. Das Hamming-Verfahren läuft wie folgt in 5 Schritten ab: 
(1) Zunächst wird mit dem Prädiktor der Funktionswert 


vi+1=ui-3 + a h (2fi — fi-1 + 2fji2) 
vorhergesagt. 
(2) Dieser Prädiktorwert wird mittels des vorhergehenden Korrektorwerts zj 
modifiziert 


112 
wi+1= ir = 777 Zi Vi) 


(3) Für diesen modifizierten Wert wird die rechte Seite der Differentialglei- 
chung bestimmt 


w'i+1 = flxi+1, wi+1) 


(4) Der neue Korrektor lautet 


zu41 = 5 (9u; — 2 + 3h (wirt +2 - A) 


(5) Der gesuchte Funktionswert uj+1 wird schließlich mit Hilfe des Korrek- 
tors verbessert 
9 


av 
121 (zi+1 i+1) 


uj+1 = zirl + 
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Ist der rechte Rand des Integrationsintervalls [a, b] noch nicht erreicht, so 
wird x um die Schrittweite h erhöht und das Hamming-Verfahren erneut 


angewandt. 


Wie man der Prädiktorformel entnimmt, benötigt man bei jedem Schritt die 
Kenntnis der Funktionswerte fi, fi_1 und fj_2. Man sagt, das Verfahren ist 
dreifach zurückgreifend. Da aber nur der Startwert ug und somit auch fg be- 
kannt ist, müssen die fehlenden Funktionswerte f1, fa mit Hilfe eines anderen 
Verfahrens berechnet werden. Man wählt hier meist die Methode von Runge- 
Kutta. 


Zum folgenden Programm 

Als Beispiel soll die gewöhnliche Differentialgleichung 
y'=y-2x/y 

mit dem Anfangswert y(B) = 1 im Intervall [@,1] integriert werden. 


Die Differentialgleichung wird in Zeile 710 in der Form 
10 F=Y—2xxX/Y 


vereinbart. Das Intervall [a,b], der Anfangswert yg und die Schrittweite h 
werden in From von DATA-Werten eingelesen. Die numerische Lösung kann 
dem Programm-Ausdruck entnommen werden. 


Zum Vergleich ist hier zusätzlich die exakte Lösung der Differentialgleichung 


yvayıx+1 


ausgedruckt worden. Der relative Fehler am rechten Intervallrand ist somit 
3.6. 106, 
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REM HAMMING-VERFAHREN 


HOME : INVERSE 

PRINT " " 
PRINT " HAMMING-VERFAHREN " 
PRINT " " 


NORMAL : PRINT 
REM EINGABE 
PRINT "Y'=Y-2#X/Y"ı PRINT 
READ A,B: REM INTERVALLGRENZEN 
READ Y: REM ANFANGSWERT DER FUNKTION 
READ H: REM SCHRITTWEITE 
x = AıYll) =Y 
PRINT " X"5: HTAB 9 
PRINT "Y" 
REM STARTPROZEDUR NACH RUNGE-KUTTA 
FOR I = 1 T03 
T=Y 
GOSUB 718 
Ki) =F 
H=H/2 
x=X+ Hi 
FOR J=1T0 3 
Y=X+ INT(J/ID) *H 
Y=T+Hi*Kt)) * INT (J/5S+D 
GOSUB 718 
KU+i1)=F 
NEXT J 
REM SCHRITTFUNKTION 
Y=T+H/6%* (KU) +2 #Ki2) +2 #KII + Ki 
) 
YyYı+1)=Y 
NEXT I 
FOR I=1T03 
Y=A+IrH 
Y=Yiıl+D 
GOSUB 718 
FII+1) = F 
GOSUB 789 
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NEXT I 


REM HAMMING-VERFÄAHREN 
XKmfA+4% Hik = diPi = 9 


REM PRAEDIKTOR 
P=Yi)+4/3#H* (2 »#F(4) - FI +2 # F(2) 
) 
REM MODIFIKATOR 
M=P- 112/121 x (k- PÜ 
YaeıM 
GOSUB 719 
M=F 
REM 
K= (9 #Yt4) - Yi2) + 34H (M +2 %*# F(A) — F(3 


I) /8B 

REM SCHRITTFUNKTION 
YS5)=K+9/121* (K-P 
GOSUB 78% 
Y=Yi(5) 

GOSUB 710 
F(5) 
FOR 
YıD 
FiD 
NEXT I 
X{=X+H 

IFX<B + H THEN 478 


END 


H 
Tr 
U 
de 


1 1 
< 
- 

+ 
_ 


REM UNTERPROGRAMM ZUR AUSWERTUNG DER DIFF.GLEICH 
UNG 
F=eY-2#X/Y 

RETURN 

DATA 9,1 

DATA 1 

DATA „1 


REM UNTERPROGRAMM ZUR AUSGABE 
IFX> = 8 THEN PRINT " "ı 
PRINT X3 

HTAB 8 

IFY> =38THEN PRINT " "; 
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830 PRINT Y 
849 RETURN 


HAMMING-VERFAHREN 


Y’aY-2#X/Y 


> 


Song >» WN 


mi 08 


Y 
1.09544553 
1.183216875 
1.26491223 
1.54156395 
1.41410298 
1.4831788 
1.5491585 
1.61243786 
1.67331976 
1.73205797 
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ANHANG A 


HISTORISCHE RECHENAUFGABEN 


Altbabylonische Keilschrifttafeln (ca. 2000 v. Chr.) 


1) Länge und Breite habe ich multipliziert und so die Fläche gemacht. Wie- 
derum was die Länge über die Breite hinausgeht zur Fläche, habe ich addiert 
und (es ergibt) 183. Wiederum Länge und Breite addiert (gibt) 27. Länge, 
Breite und Fläche (ist) was?. Lösung: (15, 12, 180) (14, 13, 182) 


2) 1/4 Breite und Länge sind 7 Handbreiten, Länge und Breite zusammen 
sind 10 Handbreiten. Lösung (6,4) 


3) Die Fläche und die Seite des Quadrates habe ich addiert und 0,75 ist es. 
Lösung: 1/2 


4) Ein Balken, 0,5 lang (steht senkrecht) an einer Wand. Oben ist er um 0,1 
herabgekommen. Von unten (wie weit hat er sich entfernt)? Lösung: 3/10 


5) Finde, wie lange es dauert, bis eine Geldsumme sich bei 1/5 Jahreszins 
verdoppelt. Lösung: 3.801784 


Papyrus Rhind (ca. 1700 v. Chr.) 


6) Haufen: sein Ganzes, 2/3 hinzu, {von der Summe) 1/3 hinweg, bleibt 30. 
Lösung: 27 


7) 7 Häuser, in jedem 7 Katzen, jede frißt 7 Mäuse, jede von diesen 7 Ähren, 
von jeder Ähre könnte man 7 Scheffel Korn ernten, wieviele Dinge sind es ins- 
gesamt)? Lösung: 19607 


8) Zu verteilen 100 Brote unter 5 Personen in arithmetischer Folge, daß 1/7 
der Summe der drei ersten Anteile gleich der Summe der letzten drei Anteile 
ist. Lösung: (38 1/3, 29 1/6, 20, 10 5/6, 1 2/3) 


9) Ein Quadrat und ein weiteres, dessen Seitenlänge 3/4 des ersten ist, bil- 
den die Fläche 100. Wie lang ist die Seite des ersten? Lösung: 8 


10) Haufen, 2/3 davon, 1/2 davon, 1/7 davon, zusammen genommen ergeben 
33. Lösung: 14 28/97 


Anthologia Graeca (ca. 500 n. Chr.) 


11) Wieviele Äpfel werden gebraucht, wenn 4 von 6 Personen 1/3, 1/8, 1/4 
und 1/5 jeweils vom Ganzen erhalten, während die 5. Person 10 Äpfel erhält 
und genau ein Apfel für die 6. Person übrig bleibt. Lösung: 120 
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12) Demarchos verlebte 1/4 seines Lebens als Junge, 1/5 als Heranwachsen- 
der, 1/3 als Mann und verbrachte 13 Jahre im Alter. Wie alt ist er geworden? 
Lösung: 60 


13) Ziegelmacher, ich muß dringend an meinem Haus weiterbauen. Es ist 
schön heute und ich brauche nicht mehr als die üblichen 300 Ziegel. Du allein 
kannst soviele am Tage machen, dein Sohn 200 und dein Schwiegerson 250. 
Wenn ihr alle drei miteinander arbeitet, in welcher Zeit könnt ihr die ge- 
ten Ziegel herstellen? Lösung: 2/5 


14) Ich bin ein Löwe aus Erz; aus meinen Augen, aus meinem Mund und aus 
meiner rechten Fußsohle springen Fontänen hervor. Mein rechtes Auge füllt 
das Becken in 2 Tagen, das linke in 3 und mein Fuß in 4 Tagen. Wie lange 
dauerts, wenn alles vereint? Lösung: 3 33/37 


Algebra des al-Hwarizmi (ca. 825 n. Chr.) 
15) Ein Quadrat und zehn seiner Wurzeln sind gleich 39. Lösung 3, —13 


16) Ein Quadrat und 21 sind gleich zehn seiner Wurzeln. Lösung 3,7 


17) Ich habe 10 in zwei Teile geteilt. Dann habe ich den einen Teil durch den 
anderen geteilt und der Quotient war 4. Lösung (8,2) 


18) Ich habe 10 in zwei Teile geteilt. Dann habe ich jeden Teil mit sich selbst 
multipliziert, das kleinere vom größeren abgezogen und der Rest war 40. 
Lösung (7,3) 


“Aufgaben zur Schärfung des Verstandes’’ von Alkuin (ca. 790) 


19) Wenn 100 Buschel Korn so auf 100 Leute verteilt werden, daß jeder 
Mann 3 Buschel, jede Frau 2 Buschel und jedes Kind 1/2 Buschel erhält, ist 
die Frage, wieviel Männer, Frauen und Kinder waren es? 

Lösung: (17,5,78) (14, 10,76) (11, 15,14) (8, 20,72) (5, 25, 70) 


20) 30 Körbe, 10 davon voll, 10 halbvoll und 10 leer, sollen so unter drei 
Söhne verteilt werden, daß die Körbe und Inhalte gleich verteilt sind. Wie 
kann dies geschehen? Lösung: z.B. 10 halbvolle 


21) Ein Hund jagt ein Kaninchen, das 150 Fuß Vorsprung hat. Der Hund 
springt immer 9 Fuß weit, wenn das Kaninchen 7 Fuß springt. Nach wieviel 
Sprüngen hat der Hund das Kaninchen eingeholt? Lösung: 75 


22) Ein Wolf, eine Ziege und ein Kohlkopf müssen in einer Fähre einen Fluß 
überqueren, die außer dem Fährmann noch einen weiteren Gegenstand beför- 
dern kann. Wie muß der Fährmann sie über den Fluß bringen, damit weder 
der Wolf die Ziege noch die Ziege den Kohl frißt? 


208 


Lilavati des Bhaskara Il (1150) 


23) Das Quadrat aus 1/8 einer Herde Affen tummelt sich im Walde, während 
die übrigen zwölf auf dem Gipfel eines Hügels brüllen. Lösung: 16, 48 


24) Jemand hat 300 Rupien und 6 Pferde. Ein anderer hat 10 Pferde vom 
gleichen Wert, aber eine Schuld von 100 Rupien. Beide haben dasselbe Ver- 
mögen. Was ist der Preis eines Pferdes? Lösung: 57 1/7 


25) Schönes Mädchen mit den glitzernden Augen, sage mir, so du die richtige 
Kunst der Umkehrung verstehst: Welches ist die Zahl, die mit 3 vervielfacht, 
sodann um 3/4 des Produkts vermehrt, durch 7 geteilt, um ein Drittel des 
Quotienten vermindert, durch Ausziehung der Quadratwurzel, Addition von 
8 und Division durch 10 die Zahl 2 hervorbringt? Lösung: 28 


26) Von einem Schwarm Bienen läßt sich 1/5 auf einer Kadambablüte, 1/3 
auf der Silindhablume nieder. Der dreifache Unterschied der beiden Zahlen 
flog nach den Blüten einer Kutuja, eine Biene blieb übrig, die in der Luft hin 
und her schwebt, gleichzeitig angezogen von dem lieblichen Duft einer Jasmine 
und eines Pandamus. Sage mir, reizendes Weib, die Anzahl der Bienen! 
Lösung: 15 


Liber abaci des Leonardo von Pisa (1202), 


27) Jemand reist in Geschäften nach Lucca, von da nach Florenz und von da 
heim nach Pisa. In jeder Stadt verdient er so viel, daß sich sein Geld jeweils 
verdoppelt; in jeder Stadt gibt er aber auch 12 Denare aus. Nach der Heim- 
kehr hat er nichts mehr. Lösung: 10 1/2 


28) Erhält ein Mann vom zweiten 7 Denare, so hat er fünfmal so viel wie der 
zweite; erhält der zweite Mann vom ersten fünf Denare, so hat er siebenmal so 
viel wie der erste. Wieviel hatte jeder ursprünglich? 

Lösung: (12 7/17, 46 2/17) 


29) Ein König sandte 30 Mann in seinen Obstgarten um Bäume zu pflanzen. 
Wenn sie 1000 Bäume in 9 Tagen pflanzen können, wie lange würden 36 
Mann für 4400 Bäume benötigen? Lösung: 33 


30) Ein Mann hinterläßt seinem ältesten Sohn 1 Bezant und 1/7 des Verblei- 
benden. Von dem Rest gab er dem zweitältesten 2 Bezant und wieder 1/7 des 
verbleibenden Rests. Diese Teilung setzte er fort, er gab dem nächstjüngeren 
1 Bezant mehr als dem vorhergehenden und je 1/7 des verbleibenden Rests. 
Wieviele Söhne waren es, wenn der letzte Sohn den Rest erhält, aber alle Söhne 
gleichviel? Lösung: 6 
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Algorismus Ratisbonensis (1450) 


31) Item ain fraw hat veigen vnd hat auch kinder vnd sy gibt iglichen kind 12 
veigen, so pleibt ir 37 feigen. Nu nympt sy dij veigen widerumb von den 
kindern vund gibt ander wais hyn vnd gibt iglichen kind 15 feigen, so zw rint 
(fehlen) ir 44 veigen. Nu frag ich, wije vil sind der feigen vnd der kinder gewe- 
sen? Lösung: (361, 27) 


32) Item: diuidatur 8 in tales 2 partes, ut veniant in diuisione 3. (Teile 8 so 
in zwei Summanden, daß ihr Quotient 4 ist). Lösung: (6 2/5, 1 3/5) 


33) Si tunc vixisses quantum vixisti et iterum tantum et dimidium et dimi- 
dium dimidij tanti, 100 annos compleuisses. (Wenn du noch solange leben 
würdest, wie du schon gelebt hast und nocheinmal so viel, und noch halb mal 
soviel und ein viertelmal soviel, so würdest du 100 Jahre alt werden). 

Lösung: 26 2/3 


34) Item ein thurnn (Turm) ist 1/4 im ertrich (Erdreich) vnd 10 schuch im 
wasser vnd 3/5 ym luft. Nu frag ich, wye lang der thurnn seij vnd uil schuch 
im wasser sey vnd vieuil im luft. Lösung: 66 2/3 


Bamberger Rechenbuch (1483) 


35) Item ist ein Faß, das hat drei Zapfen. Wenn man den ersten zieht, geht es 
aus in 2 Tagen. Mit dem anderen Zapfen geht es aus in 3 Tagen. Mit dem drit- 
ten Zapfen geht es aus in 4 Tagen. Und wenn man sie alle drei zieht, wielange 
muß es ausrinnen? Lösung: 12/13 


36) Es lauf ein Has gegen ein Holz und ein Hund lauf ihm hinten nach. Und 
wenn der Has 12 Sprünge tut, so tut der Hund 15, und der Has ist vor dem 
Hund 100 Schritt. Nun ist die Frag, wann der Hund den Hasen erläuft, in wie- 
viel Schritten? Lösung: 33 1/3 


37) Es sind zwei Gesellen, diegehen nach Rom. Einer geht alle Tage 6 Meilen. 
Der ander geht am ersten Tag 1 Meile, an dem anderen Tag 2, etc. Um alle 
Tag eine Meile mehr denn vorher. Nun willst du wissen, in wieviel Tagen einer 
soviel gegangen ist als der andere? Lösung: 71 


38) Item ein Turm gebaut nach solchen Sitten: 1/4 des Turmes ist im Erd- 
reich, 1/5 im Wasser und 100 Schuh in der Luft. Nun frag ich, wieviel Schuh 
sind im Wasser und Erdreich und wieviel Schuh sind an dem ganzen Turm? 
Lösung: 181 9/11 
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“Titelblatt von Adam Rieses Rechenbuch (2. Aufl, 1529). Es zeigt einen 
Wettbewerb zwischen einem Zahlen- und einem Linienrechner’' 


Rechenbuch des Adam Riese (1526) 


39) Item ein kauffmann zeucht (verleiht) hinweg mit gelt /gewinnt ein drittheil 
seines Hauptguts (Kapital)/vnd 4 fl. mehr/legt an hauptgut vnd gewinn ’/ge- 
winnt den vierdten theil /bringt zusammen 40. Lösung: 21 


40) Item einer spricht zum andern /wann ich noch so vil/ein drittheil/vnd ein 
vierdtheil so viel hett /so wer meines gelts über 100 fl. so viel als jetzt darun- 
der. Lösung: 55 35/43 


41) Itemein fuhrmann fehrt von Leiptzig ghen Nürenberg in 6 tagen /vnd ein 
ander fuhrman fehrt desselbigen tags aus von Nürenberg /kompt in 8 tagen 
ghen Leiptzig / in wie vile tagen komen sie zusamen? Lösung: 33/7 


42) Item 21 Personen / Männer und Frawen / haben vertrunken 81 d /ein 
Mann sol geben 5 d vnd eine Fraw 3 d. Nu frag ich wie viel jeglicher in sonder- 
heit gewesen sind? Lösung: (9, 12) 
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43) Item einer hat zween silberne Becher / vnd ein oberlied (Deckel) / so das- 
selbig auff edn ersten gesetzt wirdt / helt er viermal des andern gewicht. Wirdt 
es aber auff den andern gesetzt / so ist er dreymal schwerer dann der erste / 
und das oberlied wigt 16 loth / wieviel wigt ein jeglicher Becher in sonderheit? 
Lösung: (7 3/11, 5 9/11) 


Algebra von Christopher Clavius (1608) 


44) Um seinen Sohn zum Rechnen zu ermuntern, zahlt ein Vater seinem 
Sohn 8 Pfennige für die Lösung einer Aufgabe, wenn sie richtig gelöst wurde 
und zieht ihm 5 Pfennige ab, wenn sie falsch war. Nach 26 Aufgaben schuldet 
weder der Vater dem Sohn etwas noch umgekehrt. Wieviele Aufgaben hat der 
Sohn richtig gelöst? Lösung: 10 


45) Wenn ich jedem Bettler vor der Tür 7 Pfennige geben würde, würden mir 
24 Pfennige verbleiben. 32 Pfennige würden fehlen um jedem 9 Pfennige ge- 
ben zu können. Wieviele Bettler sind es und wieviel Geld habe ich? 

Lösung: (28, 220) 


46) Einem Diener werden 100 Denare und ein Mantel als Jahresgehalt ver- 
sprochen. Nach 7 Monaten verläßt der Diener das Haus und erhält 20 Denare 
und den Mantel. Wieviel ist der Mantel wert? Lösung: 92 


Vollständige Anleitung zur Algebra von Euler (1770) 


47) Ein Vater hinterläßt seinen 3 Söhnen ein Vermögen von 1600 Reichs- 
talern. Nach seinem Testament soll der älteste Sohn 200 Rtl. mehr haben als 
der zweite, der zweite aber 100 Rtl. mehr als der dritte; wieviel bekommt 
jeder? Lösung: (700, 500, 400) 


48) Suche ein Zahl von der Beschaffenheit, daß, wenn man sie mit 5 multipli- 
ziert, das Produkt so viel unter 40 liegt, wie die Zahl selbst unter 12. 
Lösung: 7 


49) Zwey Bäuerinnen haben zusammen 100 Eyer, die erste spricht: wann ich 
die meinigen zu je 8 überzähle, so bleiben 7 übrig; die andere spricht: wann 
ich die meinigen zu je 10 überzähle, so bleiben mir auch 7 übrig; wieviel hat 
jede Eyer gehabt? Lösung: (63, 37) (23, 77) 


50) Eine Gesellschaft von Männern und Weibern sind in einem Wirtshaus: ein 
Mann verzehrt 25 Cop., ein Weib aber 16 Cop. und es findet sich, daß die Wei- 
ber insgesamt einen Cop. mehr verzehrt haben als die Männer: wieviel Männer 
und Weiber sind es gewesen? Lösung: (7 + 16t, 11 +25 t) tganzzahlig 
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ANHANG B 
ZITATE 


Wenn jemand die Mathematiker fragt: Ihr Wunderlichen — von was für Zahlen 
sprecht Ihr eigentlich — Was werden sie wohl antworten? daß sie von Zahlen 


reden, die man nur denken kann, unmöglich aber auf irgendeine andere Art 
handhaben. 
Plato 


Geometrie ist die Erkenntnis von Gegenständen ewigen Seines. 
Plato 


Wer die Geometrie begreift, vermag alles auf der Welt zu verstehen. 
Galilei 


Ceux qui ne sont pas mathematiciens sont portes a considerer les mathema- 
tiques comme une science inhumaine. 


E. Borel 


Gott existiert, wie die Mathematik widerspruchsfrei ist, der Teufel, weil wir es 


nicht beweisen können. 
A. Weil 


Die Zahl ist das Wesen aller Dinge. 
Pythagoras 


Ich hörte mich anklagen, als sei ich ein Feind der Mathematik überhaupt, die 
doch niemand höher schätzen kann als ich, da sie gerade das leistet, was mir 


zu bewirken völlig versagt blieb. 
Goethe 


Die Mathematiker sind eine Art Franzosen: Redet man mit ihnen, so überset- 
zen sie es in ihre Sprache, und dann ist es etwas ganz anderes. 
Goethe 


Ein Mathematiker, der nicht irgendwie ein Dichter ist, wird nie ein vollkom- 
mener Mathematiker sein. 
Weierstraß 


Die Wissenschaft der reinen Mathematik in ihrer modernen Entwicklung darf 
von sich behaupten, die ureigenste Schöpfung des menschlichen Geistes zu 


sein. Whitehead 


Es gibt jedoch noch einen anderen Grund für die hohe Wertschätzung der 
Mathematik: sie allein bietet den exakten Naturwissenschaften ein gewisses 
Maß an Sicherheit, das ohne Mathematik nicht denkbar wäre. 

Einstein 
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Zur Mathematik führt kein Königsweg. 
Menachmos 


Es gibt nichts Unpopuläreres als die moderne Mathematik, und auch darin 
steckt ein Stück Symbolik der unendlichen Ferne, der Distanz. Alle großen 
Werke des Abendlandes von Dante bis Parsifal sind unpopulär. 

Spengler 


Nicht etwa, daß bei größerer Verbreitung des Einblicks in die Methode der 
Mathematik notwendigerweise viel mehr Kluges gesagt werden würde als heu- 
te, aber es würde sicher viel weniger Unkluges gesagt. 

K. Menger 


Die ganzen Zahlen hat der liebe Gott gemacht — alles andere ist Menschen- 

werk. Kronecker 

Die Zahlen sind freie Schöpfungen des menschlichen Geistes, sie dienen als 

Mittel, um die Verschiedenheit der Dinge leichter und schärfer zu fassen. 
Dedekind 


Der Begriff der Mathematik ist der Begriff der Wissenschaft überhaupt. Alle 
Wissenschaften sollen daher Mathematik werden. 
Novalis 


Wer ein mathematisches Buch nicht mit Andacht ergreift und es wie Gottes 


Wort liest, versteht es nicht. . 
Novalis 


Die Mathematik ist wie die Gottseligkeit, zu allen Dingen nütze, aber wie die- 


se nicht jedermanns Sache. 
J. Kraus 


Ich glaube ..., daß es, im strengsten Verstand, für den Menschen nur eine ein- 
zige Wissenschaft gibt, und diese ist die reine Mathematik. 
Lichtenberg 


Die sogenannten Mathematiker von Profession haben sich, auf die Unmündig- 

keit der übrigen Menschen gestützt, einen Kredit von Tiefsinn erworben, der 

viel Ähnlichkeit mit dem vom Heiligen hat, den die Theologen für sich haben. 
Lichtenberg 


Die Mathematik ist eine Art Spielzeug, welches die Natur uns zuwarf zum 
Troste und zur Unterhaltung in der Finsternis. 
d’Alembert 


So kann also die Mathematik definiert werden als diejenige Wissenschaft, in 
der wir niemalsdas kennen, worüber wir sprechen, und niemals wissen, ob das, 


was wir Sagen, wahr ist. 
B. Russell 
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Stets sind sie eilig, nur zu messen und zu rechnen, halten es für die Hauptsa- 
che, und le calcul! le calcul! ist ihr Feldgeschrei. Aber ich sage: ou le calcul 
commence, l’intelligence des phenomenes cesse: während Einer bloße Zahlen 


und Zeichen im Kopf hat, kann er nicht dem Kausalzusammenhang auf die 
Spur kommen. 
Schopenhauer 


Das Wesen der Mathematik liegt gerade in ihrer Freiheit. 
G. Cantor 


Die Geometrie ist einzig und ewig, ein Widerschein aus dem Geiste Gottes, 


daß die Menschen an ihr teilhaben, ist eine Ursache dafür, daß der Mensch 
Ebenbild Gottes ist. 
Kepler 


Die Mathematik ist eine gar herrliche Wissenschaft, aber die Mathematiker 


taugen oft den Henker nicht. , 
Lichtenberg 


Wer die Sicherheit der Mathematik verachtet, stürzt sich in das Chaos der Ge- 


danken. >. 
da Vinci 


Befaßt Euch zuerst mit Mathematik, denn die Mathematik ist die Zahl, und 
die Zahl ist Gott. 
J.H. Fabre 


Die Natur ist mathematisch, die Meisterwerke der Kunst sind im Einklang mit 
der Natur; sie drücken die Naturgesetze aus und bedienen sich ihrer. Folglich 
ist das Kunstwerk mathematisch, und der Wissenschaftler kann das strengste 
Urteil anlegen und unfehlbare Formeln anwenden. 

Le Corbusier 
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gramme ohne Wortezeit ausgeführt 


werden. 
1983.186Seiten.Spiralh.DM42,-/Fr.42.- 
S 378.-. ISBN 3.88322-023-X 


G.Daubach D.Herrmann 


Pr 


Dieses Buch enthält 40 mathemati- 
sche ProgrammeausdenBereichen: 
Mehrregister-Arithmetik - Zahlen- 
theorie - Kombınatorik - Algebra - 
Geometrie - numerische Mathemao- 
tik Neu ıst die Langzaonl-Arıthmelik. 
Sie gestattet die Grundrechenarlen 
für Zahlen bis 255 Stellen. 

In Vorb. Juni 1984. 250 Seiten. Spiralh. 


Ca.DM 42,-/ca. Fr. 42.-/ca.S 378,- 
ISBN 3-88322-036-1 


Der Apple Software-Wegweiser 
wendet sich in erster Linie an die 
Besitzer von Apple-Computern, die 
noch geeigneten Programmen fürih- 
re Anwendung suchen. Aber auch 
Händler, die ihre Kunden gezielt be- 
roten wollen, werden auf ihn nicht 
verzichten. 


1983. 304 Seiten. Kart. DM 42,-/Fr. 42.-/ 
5 378,-. ISBN 3-88322-028-0 


Diese Sammlung enthält 40 BASIC- 
Programme für den Apple Il/lle aus 
den Bereichen Finanzmathematik, 
Unternenmensforschung und Be- 
Iriebswirtschaft.. Wirtschaftliche 
Fragestellungen treffen jeden von 
uns, entweder als Steuerzahler, Spo- 
rer oder Kreditnehmer. Zahlreiche 
Beispiele erläutern die vollständigen 
Programme. 


nn 200 Seiten. Kart. DM 38,-/Fr. 38.-/ 
S 342,- 
ISBN 3-88322-035-3 


Dies ist die umfossendste Samm- 
lung deutscher und internationaler 
Hardworezusätze und Peripherie-Er- 
weiterungen für Apple-Computer. 
Sie wendet sich an Anwender von 
Apple-Computern, die nach Erweite- 
rungen ihrer Systeme suchen. Uner- 
setzlich ist sie auch für alle Händler, 
die ihre Kunden gezielt beroten 
wollen. 

Herausgeber Armin P. Landsberg. 1983. 


272 S.Kart. DM 42,-/Fr. 42.-/8 378. - 
ISBN 3-88322-025-6 . 


vl’: Pi 
Eine Einführung 
Vonder Idee 
zur Konzeption 


Dieses Buch beschreibt an Beispie- 
len ous der Praxis das Arbeiten mit 
dem wohl meistbenuizten Doten- 
banksystem der Walt. Eswerdenkei- 
ne Vorkenntnisse vorausgesetzt, 
sondern es wird Schritt für Schritt 
dos komplette Wissen vermittelt, 
dos man zum Arbeiten mil eınar 
Dotenbank braucht. 

1984. 226 Seiten 


Geb. DM 56,-/Fr. 56.-/S 498,- 
ISBN 3-88322-038-8 


Basic 
aus der Praxis 
Von passe 


Dieses Buch enthält 30 Programme 
aus den Arbeitsbereichen: Suchen - 
Schreiben - Rechnen - Sortieren. Sie 
sind on keinen Rechner gebunden, 
da sie unter CP/M und MBasic ge- 
schrieben sind. Sie entstonden aus 
der praktischen Arbeit und haben 
sich bereits bewährt. Auch der New- 
comer« kann sie ohne Schwierigkei- 
ten einsetzen. 


1983.168Seiten.Spiralh.DM40,-/Fr.40.-/ 
5 360,-. ISBN 3-88322-031-0 


Vonder Konzepton 
zum Programm 


Programme aul Deskette erhählich 


Das Buch, auch als Einstieg in das 
Programmieren geeignet, zeigt dBo- 
se ll ols mächtige Programmierspro- 
che für die Datenbankanwendung. 
Inhalt: Was ıst eine Datenbank? 
Worum dBase II? Wie wird program- 
miert? Einführung: Autokostenver- 
waltung, fürFortgeschrittene:Lilero- 
turverwoltung, kommeniierte Li- 
stings, Tips und Tricks. 

In Vorb. Juli 1984. Ca. 240 Seiten Geb. 


Ca. DM 56,-/ca. Fr. 56,-/ca. S 498,- 
ISBN 3-88322-039-6 


RB KW.Hillerkus ia 


alls der Praxis 


Programm:Beispiele fü 
kaufmannisch orientierte Berufe 


Dieses Buch umfaßt 34 Programme 
und Routinen aus den Arbeitsberei- 
chen: Suchen - Schreiben - Rech- 
nen - Sortieren - Statıstik. Die Pro- 
gramme sind unter CP/M in MBasic 
geschrieben und können somit 
leicht auf andere Rechner umge- 
setzt werden. 


a 192Seiten.Spiralh.DM40,-/Fr.40.-/ 
ISBN 3.88322-042-6 


Günther Daubach 


fün16 Bir Personalcom 
Een 86 et 


Teil 1 gıbt eine Einführung ın die 
Grundlagen: Sprochinterpreter und 
Compiler werden erläutert. Teil2 ent- 
hölt eine umfassende Beschreibung 
oller Anweisungen, Befehle und 
Funktionen, abgerundet mit zohlrei- 
chen Beispielen Teil 3 enthält eine 
Zusammenstellung alier Befehle mit 
Ihrer Syntax, weitere Beispiele und 
Tabellen 

In Vorb. Aprill984.Ca.300Seiten Spiralh. 


Ca. DM 58.-/ca. Fr. 5B.-/ca. $ 522,-. 
ISBN 3-88322-073-6 


Mit dem von der Beratungstirma In- 
terplomar entwickelten Buch kann 
der Manager den EDV-Ablauf des Be- 
triebes schnell beherrschen Das 
Monuolistdurchdie Arbeitsblätterin 
ollen Bereichen einsetzbar. 

1981. 222 S. Zahlr. Abb. u. Organisations- 
schemota. Geb. DM 128,-/Fr. 98.-/ 
S 1152,-. ISBN 3-88322-019-1. RingO. 
DM 158,-/Fr.128.-/5 1422, -. 

ISBN 3-88322-020-5. Zusötzl. Arbeitsbiät- 
ter u. Organisotionsformulore. DM 48,-/ 
Fr. 35.-/$ 432,-. Best.-Nr. IWT 205-X 


Peter PVolzing 


Dieses Buch beschreibt den Lei- 
stungsumfong von FORTRAN-8O für 
Personolcomputer, wie es im Relea- 
se 3.44 von MicroSofl festgelegt 
wurde. Alle Sprachelemente sind in 
detaillierter Form abgehandelt und 
ın Beispielen dargestellt. Das Buch 
kann ın dieser Form als Program. 
mierhandbuch und auch als Nach- 
Schlogewerk genutzt werden. 

1984 360 Seiten. 


Spiralh, DM 56,-/Fr. 56.-/S 498,- 
ISBN 3-88322-079-5 


Dieses Buch ermöglicht die schnelle 
Einorbeitung in die Programmier- 
sprache CBasic. Alle Eigenschaften 
werden mit Beispielen erklärt. Syn- 
toxdiagramme zeigen grafisch die 
Einsotzmöglichkeiten. Die Funktio- 
nen des Compilers, des Linkers und 
des Bibliotheksprogrammes LIB 
werden erläutert. 


en Sprialh.DM56,-/Fr.56.-/ 
498,- 
ISBN 3-88322-034-5 


fur 16BitPersona { 
mit CP/M B60derMSDOS 


Teil 1 gibt eine Einführung in die 
Grundlagen: Sprachinterpreter und 
Compiler werden erläutert. Teil2 ent- 
hält eine umfassende Beschreibung 
aller Anweisungen, Befehle und 
Funktionen, abgerundet mit zahlrei- 
chen Beispielen Teil 3 enthält eine 
Zusammenstellung aller Befehle mit 
ihrer Syntax, weitere Beispiele und 
Tobellen 

InVorb Aprıl1984.C0.300Seiten. Splralh. 


Ca. DM 58,-/ca. Fr. 58.-/ca. S 522,-. 
ISBN 3-88322-073-6 


Dos Buch ist in drei Teile gegliedert: 
Teil 1 führt in die Eigenschoften von 
Computern und CP/M ım besonde- 
renein. AufbauendaufdiesenKennt- 
nissen werden im 2. Teil die zentra- 
len CP/M-Hilfsprogramme vorge- 
stellt.Der3.Teilgeht-nocheinerEin- 
führung in die Funktionsweise des 
8080-Prozessors - auf die CP/M-Sy- 
stembesonderheiten ein. 

1982. 386 S. Mit zahlr. prakt. Beispielen 


Geb. DM 48,-/Fr. 48.-/5.432,- 
ISBN 3-88322-004-3 


Die Programmiersprache :C: ist be- 
sonders zur Erstellung schneller, mo- 
schinennaher Programme geeignet, 
weist jedoch gegenüber der Assem- 
blerprogrammierung wesentliche 
Vorteile auf, wie z.B. lokale und glo- 
bole Variable, Prozeduren, Funktio- 
nen usw. Diese Einführung setzt 
beim Leser keine umfassenden 
Kenntnisse voraus. 

In Vorb 1984.Ca. 250 Seiten. 


Geb. ca. DM 56,-/ca. Fr. 56.-/ca.$ 498,- 
ISBN 3-88322-041-8 


Buch von Anwendern für Laien, auch 
als Einstieg in die Computerei geeig- 
net. Inhalt: Geschichte der PC, 
Grundlagen, Arbeiten mit Disketten 
und Festplatten, Dateien, Überblick 
a. eo. Ar das Arbeiten mit 
m Standard-Hilfspr i 
und Tricks. en 


In Vorb. Junı 1984. Ca. 26 i 
Ca.DM58,-/ca. Fr. TR 522 een. 
ISBN 3-88322-072-8 “ 


Dieses Buch führt Sie zur sicheren 
Hondhabung der Tabellenkolkula- 
tion. Zahlreiche Anwendungen mit 
ousführlicher Beschreibung schlie- 
Ben sich on. Unter anderem: Werbe- 
planung, Reisekosten, Kapozitäts- 
auslostung, Baufinanzierung, Reak- 
tionszeiten.- Auch für die englische 
Version nutzbar. 

InVorb. Juni1984.Ca.250Seiten. Spiralh. 


Ca. DM 56,-/ca. Fr.56.-/ca. 5 498,- 
ISBN 3-88322-074-4 


G Daubach D Harman 


MICRD 
Er 
BASIC 


une Bar) em 


MATHE 
MAIT|IIK| 


Fertige Progratıme, 
een 
wi, 


Dieses Buch enthält 40 mathemati- 
sche Programme.ausdenBereichen: 
Mehrregister-Arithmetik - Zahlen- 
theorie - Kombinatorik — Algebra - 
Geometrie - numerische Mathemo- 
tik. Neu ist die Langzahl-Arithmetik. 
Sie gestatlet die Grundrechenarten 
für Zahlen bis 255 Stellen. 

In Vorb. Juni 1984. 250 Seiten. sproh, 


Ca.DM 48,-/ca. Fr. 48.-/ca. $S 432,- 
ISBN 3-88322-077-9 


Dieses Buch wendet sich an Super- 
Calc-Anwender- Anfänger führt es 
schnell zur sicheren Handhabung, 
Fortgeschrittene finden ausgereifte 
Anwendungen. Beispiele mit aus- 
führlicher Beschreibung machen 
das Buch auch für VisıColc-Anwen- 
der interessant. Die Unterschiede 
zwischen VisiColce und SuperColc 
sind in einem Kopitel dargestellt 

1984. u. „en Spiralh. DM 56.-/ 


Fr. 56.-/5 4 
ISBN Sr X 
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G Daubuch D Hemmann 
ae PTOgramme in zum 


MICRD 
sSOoFT_ 
BASIC 


zummmmme Barıd 2 memmmuemn 


WIRT _ 
SCHAFT 


BASKC Progumme fur un Amiendkı 
ıntgrafscher Dirsteiung 


Programme zf Dsxetteerhaitich 


Eine Hilfestellung für wirtschoftliche 
Entscheidungen sind Programm- 
sammlungen, die die guten Grafik- 
und Farbmöglichkeiten des Compu- 
ters nutzen. Diagramme, Sprites, op- 
tische Darstellungen von Simulatio- 
nen werden eingesetzt, die die Er- 
gebnisse verdeutlichen. Die finonz- 
malhemalischen Grundlagen sind 
zu jedem Programm beschrieben. 

In Vorb. Juni1984.Ca.200 Seiten. Spiralh. 


Ca. DM 48,-/ca. Fr. 48.-/ca. $S 432,- 
ISBN 3-88322-076-0 


ia % 


f L.Dodt 
] u. nm! 
ILOUISLZ3, 


u —a > 
mit Anwendungen fürdie Praxis. 
Teil 1:Kalkulation und Grafik 
Programme auf Diskette erhäktlich. 


Die Verbindung von Kalkulation, Gro- 
fik und Dotenbank verhalf Lotus zu 
schnellem Erfolg In diesem Buch er- 
lernen Sie zunächst dıe sichere Nut- 
zung von Kalkulation und Gralik. 
Zahlreiche Beispiele zeigen die Leı- 
stungstöhigkeit von 1-2-3: Umsatz- 
prognose, Inveslitionsrechnung, Til- 
gungsplon, Private Ausgaben u.a. 

1984. Ca. 250 Seiten. Geb. Ca. DM 58,-/ 


ca. Fr. 58.-/ca. S522,- 
ISBN 3-88322-085-X 


[®) 


GDauccn HDurkopf HSterner 
usa PTOgramme in mem 


MICRD 
sort 
BASIC 


sum Barıd 3 mem 


CHEMIE 


Anscegungan uno Fogemm- 
Bag an urliedenen Borechen 


Programme auf Diskette ernaltich 


Programme reichen von stöchiome- 
trischen Berechnungen bis zur Ele- 
mentedotei.Besonders genutzt wer- 
den die grafischen Möglichkeiten. 
Die Anwendungsbeispiele gehen 
von Strukturformeln, grafischen Dar- 
stellungen von Versuchsanordnun- 
gen und Funktionen bis zur 3-D Gro- 
tik im Detail. 

In Vorb. Aug.1984.Ca. 220 Seiten. Splralh. 


Co. DM 56,-/ca. Fr. 56.-/ca. 5 498,- 
ISBN 3-88322-078-7 


G. Daubach, D. Herrmann 


Mathematik 
auf dem Applell, lle,2c 


Dieses Buch enthält 40 mathematische Programme aus den Bereichen: 
— Mehr-Register-Arithmetik — Zahlentheorie — Kombinatorik — Alge- 
bra — Geometrie — Numerische Mathematik. 


Die hier neu vorgelegte Langzahl-Arithmetik gestattet die Grundrechen- 
arten für Zahlen bis 255 Stellen. Zahlreiche Anwendungen finden auch 
die angegebenen kombinatorischen Prozeduren. 


Sie finden auch das Rechnen mit Polynomen, Matrizen und komplexen 
Zahlen: Neben dem komplexen Hornerschema werden insbesonder Algo- 
rithmen zur Polynomdivision und Matrizeninversion gegeben. 


Anwendungsbezogen sind die Programme der Numerischen Mathematik. 
Hier werden Verfahren zur Lösung von linearen und nichtlinearen Glei- 
chungen gegeben, zum Eigenwertproblem von Matrizen und zur numeri- 
schen Integration und Differentiation. 

Programme auf Diskette erhältlich. 


ISBN 3-88322-.036-1 


